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FOREWORD 


The  significance  of  this  effort  to  the  Air  Force  Is  that 
It  provides  military  operations  research  analysts  with  a 
comprehensive  reference  on  computer  simulation  techniques.  This 
report  summarizes  the  concepts  and  methods  used  In  building 
simulation  models  for  operations  research  analysts  who  nave  had 
limited  experience  In  simulation  and  for  veteran  model  builders 
who  desire  an  easily  available  reference.  Additionally,  It 
provides  background  Information  for  users  of  simulation  analyses 
so  that  they  may  better  understand  the  source  of  the  data 
generated  for  their  reports. 

This  report  Is  an  In-house  effort  conducted  by  the  Air  Force 
Avionics  Laboratory,  Electronic  Warfare  Division,  Analysis  and 
Evaluation  Branch.  The  project  engineer  was  Capt  W.  K.  McQuay 
(AFAL/WRA) . 

This  report  has  been  reviewed  and  Is  approved  For  publication. 
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SECTION  I 


AN  INTRODUCTION  TO  OPERATIONS  RESEARCH 

Operations  research  (OR)  is  the  application  of  mathe- 
matics  and  scientific  methods  to  aid  the  decision  making 
otf  management.  The  implications  of  this  definition  mv/ 
tie  more  apparent  to  the  reader  through  a  closer  e  ',u»  .  :n 
of  the  specifics  of  this  statement.  Aiding  dec! sic. 
makers  is  an  inherent  characteristic  of  OR.  No  model  or 
analysis  study  is  ever  so  sufficient  unto  itself  as  to 
become  independent  of  the  judgment  supplied  by  knowledgeble 
managers.  Scientific  suggests  emphasis  on  objective 
methods  of  assessing  a  situation.  Applied  mathematics  in 
a  broad  sense  is  a  predominant  part  of  most  studies. 

First,  mathematics  is  used  as  a  tool  to  solve  problems 
and,  secondly,  it  plays  a  role  in  formulation  of  the 
problem  typically  as  a  mathematical  model. 

During  World  Wa:  ii  operational  problems  associated 
witn  radars,  aircraft,  submarines,  and  weapons  allocation 
gave  rise  to  highly  skilled  teams  o;'r  mathematicians,  physi* 
cists,  and  engineers  which  attempt  vi  ic  sui  te  the  problems. 
These  quantitative  methods  evoive.t  operations  research 

and  later  variations  such  as  syst:.  *  engines:;  ing,  manage- 
ment  science,  cost-effective  less  analysis,  and  systems 
analysis.  After  the  war  thsse  early  operations  researchers 
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transferred  their  experience  throughout  the  military 
and  industry*  Increasing  interest  led  to  the  development 
of  theory  in  linear  and  dynamic  programming,  queuing, 
gaming,  network  analysis,  inventory,  scheduling,  and 
simulation.  These  theories  ha*fe  been  applied  to 
manufacturing,  transportation,  communications,  construction, 
health  care,  banking,  and  military  operations. 

There  is  often  confusion  concerning  the  names  of 
fields  which  are  variations  of  OR.  If  emphasis  is  given 
to  planning  and  design  of  new  industrial  or  military 
systems  to  increase  performance  of  existing  operations, 
then  the  term  "systems  research"  or  "systems  engineering" 
is  used.  Dealing  with  the  problems  of  efficient  management 
or  control  of  systems  is  called  "management  science." 
Attention  to  differences  in  costs  or  resource  requirements 
among  available  alternatives  is  referred  to  as  "cost 
effectiveness  analysis."  Finally,  the  term  "systems 
analysis"  is  applied  to  any  systematic  approach  to  the 
comparison  of  alternatives.  One  should  note  that  in  the 
military , establishment ,  systems  analysis  connotes  long- 
range  planning  and  is  associated  with  problems  where  one 
decides  what  ought  to  be  done,  not  just  how  to  do  it. 

The  total  analysis  is  more  complex,  usually  qualitative, 
and  seldom  suited  to  quantitative  optimization. 


In  any  event,  little  uirtinction  remains  between 
operations  research  and  its  variations.  The  differences 
afe  a  matter  of  emphasis  and  it's  not  worthwhile 
attempting  to  firmly  distinguish  among  them.  The  term 
operations  research  will  be  used  throughout  this  report 
and  is  considered  synonymous  with  systems  analysis. 

Basic  Themes  in  OR 

Three  pervasive  and  interrelated  themes  are  found  in 
OR  literature.  First,  there  is  emphasis  on  optimization. 
Typically  the  optimization  is  constrained  so  that  values 
of  the  decision  variables  which  maximize  the  objective 
function  are  restricted  so  as  to  satisfy  certain  technological 
restraints.  Secondly,  the  analyst  seeks  derivation  of 
analytic  properties  of  the  model.  These  could  include 
sensitivity  of  an  optimal  solution  to  model  parameters,  the 
structural  form  of  an  optimal  solution  such  as  an  (s,S) 
policy  in  inventory  theory,  or  operating  characteristics 
of  the  solution,  for  example,  the  probability  of  no  bombers 
reaching  their  targets.  Thirdly,  there  is  explicit 
recognition  of  system  interaction.  The  results  of  an  OR 
analysis  cannot  be  determined  or  applied  in  isolation 
from  the  surrounding  military  or  industrial  environment, 
but  are  c  part  of  and  are  affected  by  that  environment. 

Since  OR  cuts  across  many  fields  —  engineering  science, 
economics,  physical  science,  and  biological  science  —  it 

v- 
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must  be  regarded  as  a  systems  effort. 

Emphasis  on  making  decisions  or  taking  actions  is 
central  to  operations  research  applications.  We  will  now 
consider  the  decision  making  process  and  then  the  availability 
of  aids  which  the  decision  maker  can  employ. 

The  Decision  Making  Process 

Since  military  decisions  involve  the  security  of  the 
nation  and  its  scarcest  resources,  decision  makers  must 
make  the  best  possible  decision  every  time. 

The  decision  making  process  can  be  viewed  as  follows: 

(1)  A  review  of  the  goal  or  ends  to  be  served,  (2)  a  study  / 

of  proposed  alternatives,  (3)  the  ordering  or  ranking  of  / 

i 

the  alternatives  in  some  rational  arrangement,  and  (4)  j 

the  selection  of  one  or  some  combination  of  the  alternatives.; 

i 

l 

Thus,  decision  making  may  be  defined  as  selecting  a  course  i 
of  action  from  among  a  number  of  alternatives  according  j 

i 

to  some  criterion.  In  selecting  the  best  alternative,  it 
is  not  always  immediately  apparent  which  alternative  is 
most  desirable;  In  any  case  it  is  necessary  to  develop 
a  method  of  measuring  the  effectiveness  of  various  solutions. 
This  process  alone  will  be  an  aid  to  the  decision  maker 
in  clarifying  the  objective.  It  may  also  be  necessary  to 
work  toward  subgoals.  For  example,  the  main  objective  of 
a  logistics  system  is  to  maximize  the  effectiveness  of 

operational  units.  But  a  more  tractable  goal  may  be  to 
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minimize  the  supply  backorder  rate,  supply  fill  rate,  or 
aircraft  ready  rate. 

The  statistician  and  mathematician  sht  considerable 
light  on  decision  making  through  the  application  of 
statistically  derived  models  and  game  theory  to  military 
management  problems.  Economics,  with  its  well  developed 
analytical  tools,  addresses  the  problem  of  resource 
allocation.  Information  science  focuses  attention  on  the 
vital  ingredient  of  management  decision  making  -  the  flow 
and  organization  of  data.  All  of  these  fields  contribute 
to  the  various  techniques  available  to  the  operations 
research  analyst. 

Decision  Making  Aids 

Several  methods  are  available  to  aid  the  decision 
maker  in  determining  and  selecting  among  alternatives. 
Figure  1  shows  the  relationship  of  various  decision  making 
aids.  Based  on  previous  experience  alonq,  the  decision 
maker  may  use  his  intuition  to  decide  which  solution  to 
adopt.  A  small  business  man  such  as  the  corner  grocery 
store  (if  any  remain)  may  have  developed  a  rule  of  thumb 
or  may  simply  "know"  how  many  items  to  order  to  retain 
a  certain  level  of  inventory.  After  all,  it’s  always 
worked  in  the  past.  The  inefficiency  of  such  trial  and 
error  approaches  may  never  be  apparent  to  the  decision 
maker  especially  if  the  penalty,  lost  customers  say  for  an 
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out  of  stock  condition,  is  small.  It  may  be  the  closest 
store  so  the  customers  always  return. 

A  more  vigorous  approach  is  to  enumerate  all  possible 
alternatives  and  select  the  best  solution  according  to 
some  cricerion  or  objective  function.  Generally  enumeration 

is  too  time  consuming  or  physically  impractical  for  numerous 

r 

alternatives.  The  third  possibility  is  to  model  the 
problem  using  a  descriptive,  mathematical,  physical,  or 
simulation  model. 

Descriptive  models  are  expressed  verbally  in  one’s 
native  language  and  generally  are  used  in  the  humanities 
or  social  sciences.  For  example,  once  psychologists  con¬ 
ceived  of  human  beings  as  motivated  solely  by  their  need 
to  reduce  tension  or  discomfort.  Man's  behavior  was 
explained  in  terms  of  the  tension  reduction  model  of 
motivation.  With  the  model,  man  was  described  as  an 
organism  seeking  to  avoid  discomfort. 

Mathematical  models  use  concise  mathematical  symbols 
to  describe  the  status  of  variables  in  the  system  and  the 
way  in  which  they  interact  and  change.  Linear  programming 
theory  of  operations  research  is  a  prime  example  of 
mathematical  modeling.  Another  example  is  the  classical 
transportation  problem  which  can  be  generally  formulated 
as  follows:  A  product  is  available  in  known  quantities  at 
each  of  .a  origins.  Given  quantities  of  the  product  are 
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required  to  be  shipped  to  each  of  n  destinations.  The 
minimum  cost  of  shipping  a  unit  of  the  product  from  any 
origin  to  any  destination  is  known.  The  analyst  must 
determine  the  shipping  schedule  which  minimizes  the  total 


cost  of  shipment.  The  general  mathematical  model  is: 

m  n 

min  g  r  c . .  x-  • 

T«1  p-1 


13  13 


n 

such  that  for  i  ■  1,2,...,  m  (supply) 

3-1 


where 


m 

2*  x..  -  D,  for  j 
i-1  J 


1,2,. . . ,  n  (demand) 


jr  .  ■  quantity  of  product  to  be  shipped  from 
J  origin  i  to  destination  j 

c^j  ■  cost  of  shipping  one  unit  of  the  product 
J  from  source  i  to  destination  j 

S-  ■  amount  cf  the  product  available  from  source  i 

Dj  ■  amount  of  the  product  desired  at  destination  j 

x. .  are  non-negative  integers  where  all  Sj  and 
1  Dj  are  positive  integer  such  that 


Incidentally,  the  solution  is  provided  by  the  Transportation 
Simplex  Algorithm. 
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A  primary  objective,  in  developing  a  mathematical 
model  is  to  provide  a  realistic  representation  of  the 
behavior  of  the  real  system,  whether  certain  or  random. 

In  those  cases  where  the  values  of  all  parameters  and  data 
required  by  the  model  are  known  exactly,  the  model  is 
called  deterministic.  However,  decision  making  usually 
occurs  in  an  uncertain  environment,  i.e.,  one  does  not 
always  know  for  certain  whet  will  happen  as  a  result  of 
particular  actions.  By  associating  probabilities  with 
the  occurrence  of  a  particular  event,  one  can  use  the 
statistical  results  of  the  random  process  to  form  a 
probabilistic  model. 

There  are  basically  two  approaches  in  treating 
probabilistic  phenomenon.  One  approach  is  to  model  the 
behavior  in  terms  of  the  expected  value  of  different 
states.  The  expected  values  as  functions  of  time  are 
determined  from  experimentation  or  observation  and 
approximated  with  analytical  functions.  The  expected 
value  at  a  given  time  is  calculated  from  the  function  and 
used  as  input  to  determine  the  probability  of  occurrence 
of  the  event.  The  advantage  of  the  expected  value  approach 
is  that  only  a  single  execution  of  the  program  is  required 
to  determine  model  results.  For  example,  in  engineering 
reliability  studies,  an  equipment  component  often  is  found 
to  have  a  constant  mean  time  to  failure  throughout  the 
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equipment  lifetime.  If  E(t)  -  l/X  ;  then  it  can  be  shown 
that  f(t)  ■  Xe  ,  the  probability  of  failure  at  time  t. 

By  empirically  determining »  the  above  function  can  be 
used  to  determine  the  probability  of  the  event  "component 
fails"  at  time  t. 

The  seoond  approach  employs  stochastic  sampling  pro¬ 
cedures  and  is  often  mistakenly  called  Monte  Carlo.  It 
is  distinguished  from  the  expected  value  probabilistic 
model  by  the  use  of  statistical  sampling  and  random  numbers. 

A  typical  application  would  be  determining  the  single  scan 
probability  of  detection  of  a  target  by  a  radar.  The 
probability  of  detection  (Pd)  can  be  mathematically  described 
by  an  exponential  distribution  such  as 

-k  s 

Pd  -  Ae  n 

where  A  is  an  attention  factor  for  the  radar  operator 
K  is  a  factor  accounting  .jr  the  radar  type 
s^ 

n  is  the  signal  to  noise  ratio. 

To  determine  on  a  particular  scan  if  the  target  is  detected, 
a  uniform  ranOoT,i  number  (the  sample  value)  is  generated  and 
compared  to  Pd.  If  the  random  number  is  less  than  P<j,  the 
target  is  detected;  otherwise,  there  is  no  detection. 

Physical  models  are  scaled  replicas  of  the  real  system. 
Included  are  floor  plan  layouts,  one-eighth  scale  wind 
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tunnel  aircraft  models,  or  full  scale  mock-ups  of  newly 
designed  weapon  systems.  In  the  most  trivial  sense,  a 
street  map  can  be  considered  a  physical  model  of  a  city. 

While  on  the  other  extreme  the  petroleum  and  chemical 
industry  often  build  a  half  scale  fully  operating  version 
of  a  new  refinery  to  field  test  a  new  design  or  procedure. 

The  last  and,  for  our  purposes,  the  most  important 
model  is  simulation.  Simulation  is  the  representation 
of  certain  features  of  the  behavior  of  a  physical  system 
by  the  behavior  of  another  system.  In  many  cases, 
simulation  involves  mathematical/logical  models  of  real 
systems.  In  fact,  the  development  of  a  simulation  model 
usually  starts  with  a  mathematical  model  of  the  real 
world  system.  Simulation  has  received  increased  emphasis 
in  recent  years  since  it  allows  the  examination  of  the 
dynamic  interrelationships  of  variables  and  parameters 
especially  under  continuously  changing  situations. 

The  Structure  of  Simulation 

Figure  2  shows  the  structure  of  simulation  as  it 
relates  to  method,  approach,  and  objective.  The  simplest 
simulation  model  is  a  manual  simulation  in  which  model 
behavior,  bookkeeping,  and  any  other  actions  are  accomplished 
by  the  person  participating.  To  illustrate,  consider  a 
war  game.  A  war  game  is  a  simulation  of  a  military 
operation  involving  two  or  more  opposing  forces  using 
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rules,  data,  and  procedures  designed  to  depict  a  real* 
life  situation.  An  educational  war  game  provides  the 
participant  with  decision  making  experience  while  an 
analytical  war  game  provides  the  participant  or  an 
observer  with  decision  making  Information.  In  a  manual 
war  game  all  forces  and  equipment  are  represented 
artificially  and  participants  play  the  game  by  providing 
the  movement  of  the  simulated  forces  and  making  game 
decisions.  For  instance,  toy  soldiers,  artillery,  and 
aircraft  might  be  positioned  on  a  map  to  represent  a 
battle.  The  players  could  move  the  figures  and  equipment 
over  the  map  to  simulate  different  strategies  and  thereby 
gain  insights  into  their  relative  success  or  failure. 

The  most  frequently  used  method  for  simulation  is  an 
analog  or  digital  computer.  The  digital  computer  simula¬ 
tion  model  is  the  subject  of  this  report  and  will  be 
discussed  in  depth  throughout  the  remaining  sections. 

Again  consider  the  war  game  example.  Nowadays  war  games 
arc  often  executed  on  a  digital  computer  according  to 
some  fixed  logic  incorporated  into  the  program.  In  a  pure 
digital  model  there  are  no  real  players  and  all  physical 
world  characteristics  are  pre-programmed  or  are  input. 

The  computer  allows  the  war  game  to  be  rapidly  repeated 
with  different  inputs  and  generates  information  for  the 
comparison  of  various  alternative  actions. 
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Simulation  Yodel 
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MAN-IN-THE-LOOP 


Figure  3*  A  Simplified  View  of  Three  Types  of  Simulation 


If  the  simulation  is  complex  but  human  interaction  is 
still  desired,  the  model  could  be  manned  with  the  bookkeeping 
function  of  what,  where,  and  when  as  will  as  related  numeri- 
cal  calculations  relegated  to  the  computer.  These  manned 
simulations  are  also  called  "man- in- the- loop"  or  man- 
machine  games.  Depending  on  the  nature  of  the  simulation, 
the  human  interaction  may  be  continual  decision  making  or 
only  occasional  inputs.  A  well-known  manned  simulation  is 
the  Link  Trainer  which  simulates  an  aircraft's  responses 
to  the  actions  of  a  student  pilot. 

Finally,  the  simulation  may  be  hybrid,  in  that,  various 
combinations  of  the  other  methods  and  even  physical  models 
(replicas)  are  included.  For  example,  a  tactical  air 
training  exercise  might  involve  real  aircraft  against  real 
and  simulated  radars  which  control  digital  computer  simula¬ 
tions  of  surface-to-air  missiles  or  anti-aircraft  artillery. 

As  an  aside,  one  should  note  that  the  term  "computer 
simulation"  refers'to  simulations  which  use  either  analog 
or  digital  computers  entirely  or  in  part.  Thus,  the  term 
applies  to  all  types  of  simulation  except  manual. 

The  simulation  approach  may  be  to  examine  overall 
behavior,  such  as  the  large  scale  simulation  of  an  air 
defense  system,  or  subsystem  behavior  such  as  the  flight 
dynamics  of  a  surface-to-air  missile.  Directly  related 
to  the  approach  is  the  concept  of  open  and  closed  loop 
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systems.  Figure  4  illustrates  the  two  hasic  structures. 

The  difference  between  them  is  that  the  output  of  the 
closed  loop  system  feeds  back  into  the  input. 

In  the  open  loop  system,  variations  in  the  output  are 
caused  by  variations  in  the  input  or  system  parameters. 

A  parametric  study  could  determine  the  appropriate  cause 
and  effect  relationships.  On  the  other  hand,  closed  *oop 
systems  have  their  own  dynamics  and  tend  to  be  parameter 
insensitive.  In  a  closed  loop  system  the  input  depends 
on  the  output.  The  change  in  input  due  to  output  may 
oppose  the  original  output  (negative  feedback)  and  cause 
the  system  to  underrespond.  Oppositely,  the  output  may 
feed  back  in  such  a  manner  as  to  increase  the  output  (posi¬ 
tive  feedback). 

Large  scale  simulations  tend  to  be  closed  loop  systems 
and  simulations  of  subsystem  behavior,  open  loop.  When 
applied  to  the  management  of  large  cystems,  closed  loop 
systems  theory  is  called  systems  dynamics.  The  original 
application  to  industry  and  social  systems  was  accomplished 
by  Jay  W.  Forrester  at  Massachusetts  Institute  of  Technology 
[Reference  11  and  12]. 

Simulation  may  have  one  ot  two  objectives:  analysis  or 
synthesis.  In  many  cases  '’r.e  is  nble  to  describe  in  detail 
the  characteristics  of  a  system  and  acsircs  only  to  know 
how  that  system  reacts  to  a  given  input.  Hence,  in  analysis 
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OPEN  LOOP  SYSTEM 


Figure  4 .  Open  and  Closed  Loop  Systems 
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Figure  5,  The  Objective  of  Simulation 
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one  inputs  known  values  6 f  parameters  into  a  well-known 
system  to  determine  a  previously  unknown  output.  As  an 
example,  suppose  one  models  some  ECM  techniques.  The  unknown 
output  is  data  on  the  survivability  and  expected  cost  of  a 
mission  for  aircraft  carrying  the  given  ECM  which  penetrates 
the  threat  environment.  But,  in  synthesis,  one  knows 
the  input  and  output  but  desires  a  description  of  thp 

i 

model  which  produces  that  output.  For  example,  a  given 
amount  of  money  is  allotted  for  ECM  (known  input)  and  a 
given  acceptable  aircraft  survivability  is  stipulated 
(known  output),  for  which  one  must  develop  ECM  concepts  to 
provide  that  survivability  for  the  given  dollar  cost. 

Either  analysis  or  synthesis  can  be  equally  well  ap¬ 
plied  to  an  open  loop  simulation.  But  the  situation  is 
quite  different  for  a  closed  loop  system.  For  an  open 
loop  system,  the  relationship  between  input  and  output  is 
essentially  stable  so  that  synthesis  can  determine  a  system 
description  or  analysis  can  determine  cause  and  effect. 

In  a  closed  loop,  the  output  affects  the  input  and  one  can¬ 
not  be  sure  whether  the  output  is  determined  externally  by 
the  input  or  internally  from  output  feedback.  In  a  closed 
loop  system,  synthesis  usually  fails  so  that  the  analyst 
must  employ  other  appioaches  from  feedback  control  theory. 

Digital  Computers  and  Programming 

The  complexity  of  mathematical  models  used  in  OR,  the 
volume  of  data  to  be  manipulated,  and  the  magnitude  of 
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computations  necessitate  the  use  of  digital  computers.  The 
special  significance  of  the  computer  lies  in  its  ability 
to  remember,  process,  and  provide  information  at  high  speed. 
A  brief  discussion  follows  on  the  components  of  a  computer 
Ostein  and  the  requirements  for  algorithmic  representation 
instructions  to  the  computer. 

Components  of  a  Digital  Computer 

Any  computer  system  is  made  up  of  hardware  and  soft¬ 
ware.  The  hardware  is  the  physical  equipment  or  associated 
devices  and  peripheral  equipment  in  the  computer  system. 

The  software  is  all  programs  and  routines  which  control 
or  extend  the  capabilities  of  the  computer. 

Hardware 

The  hardware  components  may  be  grouped  into  three 
major  categories:  (1)  input  unit,  (2)  output  unit,  (3)  the 
central  processing  unit  (or  main  frame).  7 he  central 
processing  unit  consists  of  three  secti?.ti:  control, 
arithmetic  and  logic,  and  primary  storage  (memory). 

The  input  unit  transcribes  the  input  data  from  a  source 
medium  such  as  punched  cards  to  an  internal  medium  such 
as  magnetic  core  storage.  All  such  input  data  is  converted 
by  the  computer  to  a  binary  form  which  can  be  stored  in  the 
computer's  memory. 
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Figure  6 .  Computer  Components 


The  output  unit  performs  the  reverse  transcription 
process  for  information  which  has  been  derived  from  the  data 
processed  within  the  central  processing  unit.  It  transcribes 
the  results  of  computation  from  the  computer  storage  medium 
back  to  an  external  medium  usable  by  the  analyst. 

The  central  processing  unit  is  the  heart  of  the 
computer.  The  control  section  decodes  and  interpret.*; 
program  instructions  stored  in  memory  and  sends  commands 
to  the  other  computer  components  to  execute  those  instructions. 
It  also  performs  the  timing  and  sequencing  function  to 
provide  for  proper  manipulation  of  the  problem  data. 

The  arithmetic  and  logic  section  performs  the  basic 
operations:  addition,  subtraction,  multiplication, 
division,  and  certain  logical  tests  and  branching.  The 
execution  of  a  single  operation  is  measured  in  nanoseconds. 

The  primary  storage  or  memory  unit  stores  data  and 
program  instructions  for  instantaneous  access  by  other  com¬ 
puter  sections.  Memory  is  organized  in  a  hierarchy  ac¬ 
cording  to  iipeed  and  cost.  Magnetic  core  is  the  fastest 
type  of  memory  with  magnetic  drum,  magnetic  disc,  mag¬ 
netic  tape,  and  punched  cards  or  paper  tape  in  descending 
order  of  cost  and  speed.  A  large  system  would  include 
several  memories  of  each  kind. 

Software 

The  software  may  be  grouped  into  three  major  categories: 
(1)  executive  system,  (2)  command  language,  and  (3)  library. 
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Figure  7.  Language  Translation  in  Computer  Programming 
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The  executive  system  is  designed  to  organize  and  regulate 
the  flow  of  work  in  the  computer. 

The  software  most  directly  affecting  the  user  is  the 
computer  program.  A  program  is  a  set  of  instructions  written 
by  a  programmer  which  tells  the  computer  what  is  to  be  done 
at  each  step  in  processing.  The  initial  program  or  source 
is  expressed  in  a  computer  programming  language  such  as 
FORTRAN  or  ALGOL.  The  source  program  is  translated  or 
compiled  into  machine  language  for  execution  by  the  computer. 
In  some  cases,  the  analyst  may  write  all  or  part  of  the 
source  program  in  assembly  language.  The  computer  program 
in  machine  language  is  called  the  object  program  and  is 
read  into  the  memory  unit.  Machine  language  consists  of 
a  set  of  codes  built  from  the  binary  digits  used  by  the 
computer.  Each  procedure  such  as  add  or  subtract,  is 
coded  by  a  combination  of  0’s  and  l's. 

Algorithms 

A  program  is  defined  as  a  logical  sequence  of  operations 
to  be  performed  by  a  digital  computer  in  solving  a  problem 
or  in  processing  data.  Since  a  computer  program  usually 
represents  an  algorithm,  we  need  to  examine  the  nature  of 
algorithms . 

An  algorithm  is  a  set  of  rules  acting  effectively 
according  to  a  known  objective  on  some  input  to  produce 
an  output  in  a  finite  time  period.  An  algorithm  has  five 
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fundamental  properties:  finiteness,  definiteness,  input, 
output,  and  effectiveness.  An  algorithm  can  always  be 
performed  in  a  finite  number  of  steps  and  is  comprised  of 
steps  which  are  precisely  and  unambiguously  specified. 

Clearly  English  with  its  numerous  ambiguities  is  not 
suited  to  representing  algorithms.  The  algorithm  operate? 
on  a  set  of  inputs  to  produce  a  set  of  outputs.  In  being 
effective,  it  should  be  possible  to  perform  the  individual 
steps  of  the  algorithm  and  produce  exact  values. 

English  Language  Description  of  an  Algorithm 

Suppose  we  wish  to  determine  the  square  root  X  of  a 
number  A.  Most  digital  computers  use  the  Newton- Raphs on 
method  or  some  variation,  i.e.,  repeatedly  apply  the 
formula: 

Xi  -  1/2  (A/Xi.j  ♦  X..x  ) 

where  X0  is  initially  some  starting  value  from  which 
we  calculate  an  approximation  to  the  real 
square  root. 

X^  is  the  approximation  at  iteration  i  (repetition 
of  the  process) 

X4  .  is  the  previous  value  of  X  for  i  2 
1-1  X0  for  i  -  1 

The  first  step  is  to  set  Xi-l  to  a  starting /value  X0. 

The  starting  value  has  a  significant  influence  on  the  num¬ 
ber  of  iterations  required  to  converge  to  the  square  root. 

£or  convenience  let's  begin  with  XQ  ■  1.  We  must  calculate 
a  new  value  of  X  and  test  whether  the  absolute  value  of 
the  difference  between  X^  and  X^.j  is  less  than  some  pre- 
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determined  tolerance.  If  they  are  within  the  tolerance 


we  stop;  otherwise,  we  do  another  iteration  with 
assigned  as  the  value  of  X^.^.  The  flowchart  in  Figure 
8  concisely  describes  the  algorithm  in  a  clear  and 
unambiguous  fashion.  It  can  then  be  implemented  on  any 
computer  by  programming  in  a  language  such  as  FORTRAN  or 
ALGOL. 

A  FORTRAN  Program  of  the  Algorithm 

The  same  algorithm  as  a  FORTRAN  subroutine  would  be 
as  follows: 

SUBROUTINE  SQRT  (A,X) 

TOLERN  »  l.E-6 
XI  *  1 

3  X  -  . 5* (A/Xl  +  XI) 

IF  (ABS(X-Xl).LE. TOLERN)  GO  TO  7 
XX  -  X 
GO  TO  3 

7  PRINT  8,  A,  X 

8  FORMAT  (IX,  "THE  SQUARE  ROOT  OF,"  F10.5,"  IS",  F10.5) 
RETURN 

END 

The  subroutine  is  called  from  the  main  program,  calcu¬ 
lates  the  square  root,  writes  the  answer  on  the  printer,  and 
returns  to  the  main  program  with  the  calculated  value  of  X. 

Algorithms  are  necessary  since  we  must  communicate 
with  the  computer  in  an  unambiguous  fashion.  Therefore 
computer  programming  languages  allow  the  analyst  to  describe 
a  process  or  operation  in  algorithmic  form. 
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SECTION  II 


THE  ESSENTIALS  OF  COMPUTER  SIMULATION 

In  Section  I  we  sketched  an  outline  of  operations 
research  and  how  siaulation  fits  in  as  a  decision  aaking 
aid.  The  remainder  of  this  report  is  concerned  specific¬ 
ally  with  digital  computer  siaulation  models--* how  to 
build  them,  use  thea,  and  interpret  their  outputs.  In 

this  section  we  will  discuss  when  computer  simulations 

« 

are  used,  the  advantages  and  disadvantages  of  simulation 
models,  and  general  guidelines  for  a  simulation  study. 

Why  Simulate? 

On  his  first  exposure  to  computer  simulation,  an  in¬ 
dividual  may  validly  pose  the  question  -  "Why  should  I 
use  a  simulation  model?"  One  of  the  primary  reasons  for 
using  simulation  is  that  it  is  either  impossible  or  very 
costly  to  observe  detailed  behavior  of  the  real  world 
system.  For  instance,  without  siraulatirn  it  would  be 
virtually  impossible  to  observe  the  effects  of  electronic 
countermeasures  employed  onboard  aircraft  penetrating 
a  large  air  defense  system.  Another  reason  may  be  that 
the  real  world  system  is  so  complex  that  it  is  impossible 
to  describe  it  in  terms  of  a  set  of  mathematical  equations 
suitable  for  analytic  solution.  Simulation  makes  it 
possible  to  experiment  with  the  dynamics  of  the  system 
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and  study  the  complex  interaction  of  subsystems.  Through 
parameter  variations,  one  can  change  the  system's  environ* 
ment  and  observe  the  effects  on  the  model's  behavior.  A 
simulation  study  can  yield  valuable  insight  into  which 
variables  and  their  interrelationships  are  most  important. 
Dynamic  Systems  may  be  studied  in  real  time,  expanded 
time* or  compressed  time.  Lastly,  simulation  models  may 
be  used  as  teaching  devices,  for  experimentation  in  new 
situations  which  the  reel  world  system  has  not  yet  en- 
countered,  or  for  tests  of  new  policies  or  strategies  for 
operation  of  the  system.  In  a]'1  cases,  simulation  is  an 
effective  means  of  generating  data  which  would  be  other¬ 
wise  difficult  to  obtain. 


Advantages  and  Disadvantages  of  Simulation 

Computer  simulation  possesses  particular  characteristics 
which  provide  both  advantages  and  disadvantages  for  the  user 
[  8  ].  To  begin  on  an  optimistic  note,  let's  look  at  the 
advantages.  A  simulation  is  completely  repeatable  since 
the  analyst  has  complete  control  over  model  development, 
input  data,  and  execution  of  the  program.  The  model  is 
an  ideal  system  from  the  standpoint  of  collection  and 
processing  of  data.  Also,  a  simulation  model  is  free 
from  the  physical  limitations  of  the  system  being  studied. 
For  instance,  a  machine  repair  simulation  can  be  run 
continuously  for  several  hundred  thousand  seconds  of 
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simulated  time  without  worry  about  the  repairman  getting 
tired.  Finally,  the  analyst  is  able  to  build  realism 
into  the  model  only  constrained  by  computer  size  and  cost. 

Simulation  also  has  its  disadvantages.  A  simulation 
model  is  artificial  since  it  expresses  natural  phenomena 
in  purely  symbolic  terms.  It  can  be  inflexible  since 
slight  changes  in  the  purpose  of  the  model  could  cause 
drastic  changes  in  the  computer  program.  Lastly,  a  simu¬ 
lation  study  can  be  lengthy  and  the  models  quite  complex 
and  expensive. 

The  analyst  must  determine  if  simulation  is  the  ap¬ 
propriate  tool  for  his  particular  problem.  It  should  be 
not  only  applicable  but  the  lowest  cost  computational 
procedure  for  solution.  Furthermore,  the  type  of  model 
used  should  have  output  which  can  be  relatively  easily 
interpreted  by  those  who  will  use  the  results. 

General  Guidelines  for  a  Simulation  Study 

The  general  approach  used  in  a  siT  .ation  study  is 
basically  the  steps  of  the  scientific  method.  In  certain 
ways,  a  simulation  study  is  no  different  from  any  other 
system  study  using  some  other  technique.  A  flowchart  of 
the  steps  is  shown  in  Figure  9. 

The  approach  begins  with  formulation  of  the  problem. 

It  may  so  happen  that  the  initial  statement  of  the  problem 
differs  considerably  from  the  final  version  since  re- 


Figure  G.  Steps  in  a  Simulation  Study 
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finement  of  the  objectives  usually  occurs  throughout  the 
study.  Additionally  in  this  step,  the  analyst  must  de¬ 
cide  on  a  set  of  criteria  for  evaluating  the  degree  to 
which  the  objectives  are  fulfilled  by  simulation  experi¬ 
ments  conducted  in  the  study. 

The  second  step  is  to  collect  real  world  data.  Some 
preliminary  data  was  probably  already  collected  in  order 
to  formulate  t!-e  problem.  Additional  data  is  processed 
in  some  fashion  so  that  it's  in  a  suitable  form  to  aid 
in  formulating  a  mathematical  model. 

In  formulating  the  mathematical  model,  the  analyst 
must  consider  the  number  of  variables  to  include,  how 
complex  the  model  should  be,  the  computational  efficiency 
of  the  equations  or  mathematical  techniques  employed,  com¬ 
puter  programming  time  required  to  implement  the  model, 
the  degree  of  realism  required  or  permissible,  and  the 
compatibility  of  the  model  with  the  objectives  of  the 
study.  The  greatest  flexibility  of  mod'’  application  and 
ir^ight  into  the  operation  of  the  system  results  from 
models  which  are  modularly  constructed.  The  overall 
system  is  represented  by  a  set  of  submodels  of  the  in¬ 
dividual  components  of  the  system.  Through  such  an  ap¬ 
proach,  the  mind  is  able  to  grasp  the  basic  relationships 
between  subsystems,  and  the  programmer's  burden  is  some¬ 
what  lessened  since  he  can  independently  check  out  the 
submodels , 
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The  fourth  step  is  estimation  of  parameters.  The 
parameters  are  estimated  by  statistical  inference  from 
the  real  world  system  or  from  design  specifications  when 
the  system  does  not  yet  exist.  In  many  cases,  the  purpose 
of  the  study  is  parameter  sensitivity  analysis,  ice., 
determination  of  crucial  parameters  and  the  model's  reset 
tion  over  a  particular  range  of  values.  For  such  studies 
the  real  world  values  are  either  unknown  or  not  easily 
attainable  and  the  objective  is  to  determine  which  vari¬ 
ables  are  important  to  accurately  estimate. 

The  next  step  is  model  acceptance.  It  is  important 
to  examine  the  model  structure,  its  mnderlying  assumptions, 
and  the  estimates  of  parameters  for  accuracy  and  relevancy. 
The  analyst  must  determine  if  all  pertinent  variables  have 
been  included  and  that  none  currently  in  use  are  superfluous, 
if  relationships  between  variables  are  correctly  formula¬ 
ted,  and  the  level  of  detail  is  satisfactory  for  the 
study.  If  the  analyst  and  decision  maker  are  assured  that 
the  model  meets  their  requirements,  the  study  continues; 
otherwise,  they  must  return  to  step  1. 

When  the  mathematical  model  is  accepted,  it  is  pro¬ 
grammed  for  the  computer.  The  structure  of  the  program 
is  influenced  by  the  availability  of  programming 
languages:  either  general  purpose  languages  (GPL)  or 
simulation  programming  languages  (SPL).  The  choice  of 


33 


languages  is  affected  by  the  experience  of  the  progzammers, 
zompatibility  of  the  mathematical  model  and  the  language, 
and  the  type  and  form  of  output.  In  Section  VII,  we 
will  discuss  the  choice  between  GPL  and  SPL  in  more  depth.- 

Once  the  program  is  written,  it  must  be  tested  and  de¬ 
bugged.  Depending  on  its  complexity,  this  step  could  be 
the  most  frustrating  and  time  consuming.  At  this  time, 
the  analyst  must  also  define  the  initial  state  of  the  model 
and  determine  the  number  of  samples  needed  for  statistical 
analysis. 

The  eighth  step  is  program  verification.  Does  the 
program  represent  the  model?  Are  there  biases  which 
influence  the  results?  Is  the  output  suitable?  Does  the 
model  behave  the  way  the  analyst  intended? 

If  the  program  is  acceptable,  the  study  can  proceed 
with  experimentation,  analysis,  and  validation  of  model 
output.  Validation,  the  process  of  testing  the  agreement 
between  the  behavior  of  the  simulation  -«odel  and  the  real 
world  system,  should  be  a  continuing  effort  which  is  part 
of  all  simulation  studies.  Each  study  with  the  features 
unique  to  its  particular  objective  may  subject  the  model 
to  a  previously  unencountered  environment  and  test  the 
model's  validity  under  such  conditions.  The  results  of 
the  analysis  and  validation  can  feedback  vital  informa¬ 
tion  for  model  design. 
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The  ten  steps  above  outline  the  basic  approach  for  a 
simulation  study.  They  serve  as  a  general  guideline  for 
the  iterative  procedures  involved  in  building  a  computer 
simulation  model. 

The  Transformation  Effect 

Realism  and  fidelity  are  persistent  problems  for  the 
analyst.  The  model  must  not  only  behave  like  the  real 
world  system  but  do  it  well.  However  one  must  be  aware 
that  throughout  each  step  in  the  model  building  process , 
the  true  system  has  undergone  a  transformation. 

The  true  system  is  perceived  by  the  analyst  and  des¬ 
cribed  initially  as  a  natural  language  model.  The  mere 
description  of  the  true  system  is  a  transformation  from 


True  System  Perceived  System 

the  real  world.  Then  the  perception  must  be  expressed 
as  a  symbolic  or  mathematical  model. 

co — *<55) — kBSiS) 

True  System  Perceived  Expressed 

System  Perception  of 

True  System 

Then  finally  the  mathematical  model  is  written  as  a 
computer  program  or  code. 
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True  System  Perceived  Expressed  Coded  Ex- 

System  Perception  of  pression  of 

True  System  True  System 

The  true  system  has  undergone  three  transformations  from  its 
initial  state, to  the  final  representation- as  a  computer 
program.  Consequently,  a  simulation  model  is  necessarily 
an  approximate  representation  of  reality  -  an  abstraction 
of  the  real  world. 

Primary  Areas  of  Concern 

In  building  a  probabilistic  computer  simulation 
model,  the  analyst  has  four  primary  areas  of  concern: 

•  Representing  random  behavior 

•  System  representation,  including  event 

specification  and  the  time  flow  mechanism 

« 

•  Extract  and  communication  of  model  results,  and 
9  Statistical  analysis  of  model  output 

The  next  three  sections  of  this  report  discuss  each 
area  in  considerable  -epth. 


SECTION  III 


REPRESENTING  RANDOM  BEHAVIOR 

Most  Applications  of  siaulation  aodels  incorporate 
random  behavior.  For  example,  arrival  and  service  tiaes 
in  a  queuing  aodel,  target  detection  on  a  radarscope,  or 
vehicle  tracking  ermor  by  a  weapons  controller  require 
sampling  froa  probability  distributions. 

In  Section  I,  the  terms  Monte  Carlo  and  stochastic 
sampling  were  introduced  in  the  discussion  of  probabilistic 
aodels.  Although  soae  analysts  use  them  interchangeably 
there  is  a  clear  distinction.  Monte  Carlo  analysis  is 
used  to  solve  a  deterministic  analytic  problem  by  con¬ 
verting  it  to  a  .probabilistic  problem  having  a  similar 
mathematical  formulation.  Then  random  sampling  techniques 
are  used  to  estimate  the  solution  to  the  original  deter¬ 
ministic  problem.  On  the  other  hand,  stochastic  sampling 
techniques  in  a  computer  simulation  model  are  used  on 
dynamic  problems  which  have  no*  closed  form  mathematical 
representation.  The  simulation  model  is  required  to 
represent  the  observed  behavior  of  the  real  system  and 
stochastic  sampling  is  the  means  of  providing  certain 
aspects  of  that  behavior. 

The  confusion  between  Monte  Carlo  methods  and 
stochastic  sampling  in  simulation  arises  because  both 
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use  random  numbers.  To  illustrate  the  use  of  Monte  Car¬ 
lo  methods,  consider  the  area  beneath  the  irregular  curve 
below.  The  area  is  represented  by  the  integral  1  ffX). 

8t*: . /»*>  8 


Suppose  we  construct  the  rectangle  A  around  the  curve  and 

begin  to  generate  random  points  within  the  rectangl*. 

If  the  number  of  points  (P^)  within  the  region  bounded 

by  the  curve  is  compared  to  the  total  number  of  points 

generated  (P),  the  probability  of  a  point  landing  beneath 

the  curve  is  found  to  be  P  /  P.  Then, 

c 

f(x)  ■  Pc/P  x  Area  of  rectangle  A. 

We  have  taken  a  deterministic  problem  which  was  diffi¬ 
cult  to  evaluate  and  converted  it  to  a  probabilistic  analog. 
The  solution  was  then  found  using  sampling  techniques. 

In  simulation  models,  random  variables  are  used  to 
represent  the  behavior  of  those  factors  in  the  system  whose 
real  world  counterpart  fluctuates  in  an  unpredictable  but 
statistically  describable  way.  For  instance,  the  inter¬ 
arrival  time  (time  between  arrivals)  of  jobs  at  the 
base  computer  center  cculd  be  exponentially  distributed 
with  a  mean  which  can  be  estimated.  Random  variables 


38 


■ay  also  be  used  to  represent  a  deterministic  process 
when  too  detailed  a  model  would  be  required.  Adequate 
results  could  bo  obtained  with  a  simpler  model  by  ag¬ 
gregating  the  process  into  one  statistically  determined 
effect.  A  model  of  a  man  operating  a  radarscope  would 
be  too  complex. if  all  details  of  his  functions  and  work¬ 
load  were  incorporated.  So  in  large  scale  campaign  models 
an  exponential  distribution,  which  is  a  function  of  radar 
type,  signal-to-noise  ratio,  and  with  a  multiplier  for 
operator  attention,  is  used  for  determining  target  prob¬ 
ability  of  detection. 

Pseudorandom  Number  Generation 

Random  numbers  are  generated  either  by  a  device 
separate  from  the  computer  or  by  a  deterministic  calcu¬ 
lation  within  the  computer.  Only  the  latter  type  of 
generation  will  be  discussed  here  since  it  is  the  most 
frequently  used.  Pseudorandom  numbers  are  defined  as 
numbers  produced  by  a  non-random  process  which  demonstrates 
sufficiently  random  behavior  in  a  statistical  sense. 

For  simulation  purposes,  random  numbers  generated 
by  an  acceptable  method  should  be:  uniformly  distributed, 
stochastically  independent,  reprooucible ,  non-repeating 
for  any  desired  length,  capable  of  high  speed  generation, 
and  require  little  internal  storage.  Lehmer  las  been 
quoted  as  describing  such  numbers  as  "a  vague  notion 
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embodying  the  idea  of  a  sequence  in  which  each  term  is 
unpredictable  to  the  uninitiated  and  whose  digits  pass  a 
certain  number  of  tests,  traditional  with  statisticians 

l 

And  depwnds  somewhat  on  the  uses  to  which  the  sequence 
is  to  t>e  put"  {15, P.240  |. 

Mfcddle-of *the*Square*Method 

The  first  deterministic  procedure  proposed  to  generate 
pseudorandom  numbers  was  the  middle* of* the -square  method. 
Some  initial  value  Xq  is  specifiedaalong  with  a  number  of 
digits  n  in  a  base  such  ns  7  or  10.  The  algorithm 
follows : 

1.  X1  4-  xQ2 

2.  Xj  n(x|)  Where  n  selects  the  central  n  digits 


4.  If  more  values  are  needed,  go  to  1;  otherwise  STOP. 
The  middle-of-the-square  method  can  lead  to  very  short 
sequences  of  non-repeating  values.  Aq  -  30  base  10 

which  is  42  base  7.  The  square  (42) 2  is  2424  base  7 
which  is  42  base  7  [Reference  27], 

ITERATION 

step - 1 — - r 


1. 

x1*-  x2 

1  9 

2424 

2424 

2. 

X^nCxj 

)  *2 

42 

3. 

X 

o 

t 

X 

h-* 

42 

42 

The  sequence  has  length  1  thus  totally  unsuitable.  How¬ 
ever,  with  other  choices  of  arguments  very  long  sequences 
of  values  have  been  generated  [Reference  15]. 

Lineey  gongruential  Methods 

In  1947,  Lehaer  proposed  the  linear  congruential 
aethod  and  it  swiftly  replaced  the  niddle-of-the-square. 
The  linear  congruential  aethod  is  based  on  nodular  arith- 
aetic  froa  abstract  algebra  and  nuaber  tbeorv.  The  fol¬ 
lowing  discussion  presuaes  a  basic  understanding  of  aodu- 
lar  arithaetic.  Three  congruential  techniques  exist: 
aultiplicative ,  mixed,  and  additive. 

Multiplicative  Congruential  Method 

The  aultiplicative  congruential  generator  is  of 
the  form  X^+^  *  a  X  (nod  a)  n£0 

where  aq>  9,  a>  0. 

The  generator  has  a  full  period  when  the  starting  value 
Xq  (also  frequently  called  the  seed)  is  relatively  prine 
to  a  and  a  is  a  primitive  root  aodulo  m.  If  a  is  prine, 
the  period  will  be  of  length  a-1.  The  above  recursion 
means  to  multiply  the  last  random  nuaber  Xq  by  the  con¬ 
stant  a  and  take  the  result  aodulo  a. 

Any  number  sequence  generated  by  a  deterministic 
process  with  finite  input  will  eventually  repeat  when  the 
input  is  the  same  as  at  soae  previous  stage.  The  nuaber 
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of  entries  before  such  a  repeat  is  called  the  period.  • In 
selecting  a  pseudorandom  number  generator,  we  desire  a 
period  as  long  as  possible. 

A  natural  choice  for  the  modulus  is  the  computer 
word  length.  If  a  number  exceeds  the  word  length,  trun¬ 
cation  results  and  the  remaining  digits  are  the  residue 
modulo  m.  Table  1  lists  the  word  lengths  of  various 
computers  in  bits  where  the  sign  bit  is  neglected.  (Since 
overflow  will  probably  cause  the  sign  bit  to  change, 
some  type  of  bit  manipulation  such  as  the  FIELD  function 
in  UNIVAC’s  FORTRAN  V  should  be  used  to  zero  the  sign 
bit.  The  absolute  value  function  IABS  does  not  necessarily 
change  the  sign  bit.) 

Naylor  et  al  [Reference  24)  show  that  for  a  choice  of 

L 

modulus  of  the  form  2  ,  the  maximum  period  attainable  is 

2°  .  Jansson  [Reference  17)  demonstrates  that  the 

maximum  period  occurs  when  a=3  or  5  (mod  8)  with  Xq  odd. 

Finally,  Greenberger ' s  formula  [Reference  IS,  p . 238]  for  an 

approximation  to  serial  correlation  p  is 

1  6c 

p  =  a  -  am  (1  -  c/m)  +  K 

Values  of  a  near  4T  will  yield  small  values  of  p.  Suppose 

35 

the  computer  word  length  is  2  and  we  choose  as  the  initial 
value  X0=1907.  Since  the  square  root  of  2^  is  near  2^, 
a  is  chosen  as  2 ^  +  3  *  131075.  The  period  will  be  2^ 
which  should  be  sufficiently  large  for  most  uses. 
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TABLE  1.  COMPUTER  WORD  LENGTHS 


COMPUTER  MODEL  WORD  LENGTH  WITHOUT  SIGN 

(Bits) 


Burroughs  B550C 
CDC  1604,  3600,  3800 
CDC  6000  series 
PDP-6 

GE  200  Series 
GE  400  Series 
GE  600  Series 
Honeywell  800,  1800 
IBM  7040 
IBM  7090,  7094 
IBM  360,  370 
RCA  Spectra  70 
UNI VAC  1108 


239 

247 

2S9 

23S 

219 

223 

23S 

244 

2  35 
23S 

231 

231 
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C  *  *  *  MULTIPLICATIVE  CONGRUENT I AL  *  *  * 

t 

C  •  *  *  RANDOM  NUMBER  GENERATOR  *  *  * 

C 

C  ISTART  -  THE  INITIAL  STARTING  VALUE 
C  I A  -  THE  MULTIPLIER 

C  N  -  NUMBER  OF  RANDOM  NUMBERS  TO  BE  GENERATED 

SUBROUTINE  RMULT  (ISTART,  IA,  N) 

COMMON  RNUM(5000) 

I  -  ISTART 
DO  20  J  -  1,  N 
I  -  IA  *  I 

( 

FLD(0,1,I)  -  0 

20  RNUM  (J)  -  1/2.0  ••  35 

RETURN 
END 

Figure  FORTRAN  Subroutine  for  Multiplicative  Congru- 

ential  Generator. 
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Figure  10  shows  the  FORTRAN  IV  subroutine  for  e 
nultiplicative  congruentiel  generator  on  a  UN I VAC  1108 
with  word  length  2^. 

Mixed  Congruentiel  Generator 

The  nixed  congruentiel  .generator  is  of  the  forn 

xn*l  ■  a  *n  ♦  C  (nod  n). 

Jansson  [17]  shows  that  the  generator  for  nodulus  2^  has 
naxinun  period  when  a-1  (nod  4)  and  C  ■  1  (nod  2).  Hull 
and  Dobell  [16]  tested  over  a  thousand  different  nultipliers 
of  the  forn  2s  ♦  1  and  2s  ♦  3.  The  generators  were  ac~ 
ceptable  provided  that  the  nultiplier  did  not  satisfy 
either  a  ■  1  (nod  213)  or  a  <  30  and  a  4  2**  ♦  1.  The 
choice  of  C  had  snail  effect  but  nore  complicated  values 
of  C  tended  to  inprove  perfornance. 

Other  enpirical  evidence  has  indicated  that  a  equal 
to  27  ♦  1  has  been  a  satisfactory  choice  [17,24,30],  The 
initial  value  1907  night  be  used  with  a  equal  to  129  and 
C  equal  to  1.  The  nodulus  can  be  conveniently  chosen  to 
be  the  conputer  word  length.  Figure  11  shows  the 
FORTRAN  IV  subroutine* 


Additive  Congruentiel  Generator 

The  additive  generator  is  of  the  fern 


*n*l 


(nod  n) 
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C  *  *  *  SUBROUTINE  FOR  MIXED  CONGRUENT IAL 
C  *  *  *  RANDOM  NUMBER  GENERATOR 
C  I  START  *■  INITIAL  STARTING  VALUE  • 

C  IA  -  THE  MULTIPLIER 

C  IC  -  THE  CONSTANT  TERM 

C  N  -  THE  NUMBER  OF  RANDOM  NUMBERS  TO  BE  GENERATED 

SUBROUTINE  RMIX  (ISTART,  IA,  IC,  N) 

COMMON  RNUM  (5000) 

I  «  ISTART 
D020  .!  -  I,  N 
I  *  IA  J  I  ♦  IC 
FLD  (0,1, lj  -  0 
20  RNUM  (J)  -  1/2.0  **  35 
RETURN 
END 

Figure  11.  FORTRAN  Subroutine  for  Mixed  Congruential 
Generator 
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A  special  case  is  .he  generalized  Fibonacci  generator 

which  is  defined  as 

x  *  x  *  x  (nod  n) 

nU  n  ,  j>i. 

Soae  j  values  are  required  initially  before  the  generator 
can  be  used.  Jansson  U7J  has  examined  the  periods  of 
Fibonacci  generators  for  various  moduli  and  values  of  j . 
Depending  on  the  nature  of  the  polynomial,  i.e.,  primi¬ 
tive  ,  irreducible  or  reducible,  the  period  may  be  inde¬ 
pendent  or  dependent  on  initial  values.  Green,  Smith, 
and  Klem  [14]  also  examined  generators  for  j  values  be¬ 
tween  2  and  16  and  suggested  j  ■  16  as  the  smallest  value 
for  acceptable  pseudorandom  numbers.  They  described  the 

period  as  equal  to  k  2b  1,  where  k  is  a  constant  (a 

n 

table  of  kft  values  was  given).  For  j-16  and  b*35,  the 
period  is  255  x  234. 

Primitive  Polynomials  Modulo  Two 

R.  C.  Tausworthe  [32]  has  developed  a  method  for  generating 
random  binary  bit  patterns  by  the  addition  of  primitive  poly¬ 
nomials  modulo  2.  For  example,  using  polynomials  of 
degree  4,  the  recursion 

An  ♦  1  "  An-3  +  An-4  2) 

would  be  used.  Suppose  the  starting  values  are  1000,  0110, 

1101,  and  0111.  The  sequence  is 
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n  n 


C  *  *  *  ADDITIVE  RANDOM  NUMBER  GENERATOR  *** 

C 

C  N  •  NUMBER  OF  RANDOM  NUMBERS  TO  BE  GENERATED 

C  K  -  NUMBER  OF  INITIAL  VALUES  TO  BE  GENERATED  BY  THE 

MULTIPLICATIVE  GENERATOR 


SUBROUTINE  RADD  (N,K) 

COMMON  RNUM  (5000) 

DIMENSION  KORE  (79) 

11  -  0 

12  -  K-l 

C  MULTIPLICATIVE  GENERATOR  FOR  INITIAL  K  VALUES 

L  -  1907 
DO  10  J  -1,  K 
L  -  131075  *  L 
FLD  (0,1, L)  -  0 
10  KORE  (J)  -  L 

DO  20  J  -  1,  N 

C  INCREMENT  INDEXES  OF  VALUES  TO  BE  ADDED 

11  -  II  ♦  1 

12  -  12  ♦  1 

IF  (Il.GT.K)  II  -  1 
IF  (I2.GT.K)  12-1 
C  ADD  TWO  STORED  VALUES 

LAST  -  KORE  (II)  ♦  KORE  (12) 

FLD  (0,1, LAST)  -  0 
KORE  (II)  -  LAST 

C  TRANSFORM  TO  UNIFORM  (0,1) 

20  RNUM  (J)  -  LAST/2.0  **35 

RETURN 
END 

Figure  12  FORTRAN  Subroutine  for  Additive  Random  N umber 
Generator 
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1000  1 

1101  f  ^tart^n®  Values 

0111  J 

0110 

1011 

1010 

0001 

1101 

For  primitive  polynomials  of  degree  a  the  period  is  2*-l. 
Tests  for  Randomness 

If  the  numbers  generated  by  some  pseudorandom  tech¬ 
nique  are  random,  they  must  be  uniformly  distributed  and 
stochastically  independent.  The  following  tests  determine 
if  there  is  statistical  justification  for  rejection  of  a 
sequence  of  numbers  as  a  random  sample  from  a  uniform 
distribution.  A  test  may  be  based  on  either  the  property 
of  independence,  or  uniformity,  or  possibly  both.  Each 
test  is  briefly  explained  and  results  of  the  test  on 
each  of  the  three  types  of  congruential  generators  is 
discussed.  The  generators  are  used  with  the  initial 
values  suggested  above  for  a  bit  word.  The  additive 
generator  used  79  starting  values  which  were  initially 
produced  by  a  multiplicative  generator. 

All  of  the  following  tests  are  also  described  in 
Knuth  [19] and  many  are  available  in  statistics  packages 
associated  with  most  computer  library  files.  The  particular 
programs  used  here  are  described  in  Reference  25. 
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Chi  Square  Test 

The  Chi  Square  Test  is  one  of  the  earlier  Methods 
of  statistical  inference,  originally  proposed  in  1900 
by  Karl  Pearson. 

Consider  a  finite  number  X  of  mutually  disjoint 
sets  'categories  )A,  A2,  ...»  A^.  Let  P(A^)  ■  P^,  the 
probability  that  the  outcome  of  the  random  experiment  is 
an  element  of  A^.  The  random  experiment  is  to  be  repeated 
n  independent  times  and  e^  is  the  number  of  times  the  out* 

come  is  an  element  of  set  A^,  i.e.,  e^,  e2 . e^  ■ 

n-e^-e2. . . -e^.j  are  the  frequencies  with  which  the  out¬ 
come  is  an  element  of  ,  A^,  A^.  The  joint  pdf  is 

multinomial  with  parameters  n,  p^,  P^^. 

If  the  null  hypothesis  is  Ho:  P.  »  P,  for  all  i 

Ha:  all  alternatives 
7  n  2 

If  H  is  true,  the  random  variable  X*  ■  rv  (ei-np) 

0  L  an 

i-1 

has  approximately  a  Chi  Square  distribution  with  K-l 

2  2 

degrees  of  freedom *with  critical  region  X4 > X  d,  K-l. 
Pearson's  Chi  Square  criterion  is  equivalent  to  the 
likelihood  ratio  test  for  large  samples  since  -2  In  X 
(where  X  is  the  likelihood  ratio)  has  an  approximate 
Chi  Square  distribution  with  K-l  degrees  freedom.  The 
Chi  Square  test  is  a  test  of  uniformity. 
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All  generators  were  subjected  to  a  Chi  Square  test 

for  a  sequence  of  5000  nuabers  as  a  whole  to  test  global 

behavior,  then  for  groups  of  1000  and  500  to  test  local 

2  2 

behavior.  The  null  hypothesis  was  taken  as  Hfl:X  <  *JL  «, 
■-1  and  the  alternative  H^:  X2a-lJ  ^C2a,  a-i.  Table  2 
below  shows  the  Chi  Square  statistics  for  global  behavior 


of  all  generators 

and  for  which  the  nul.< 

i  hypothesis 

can- 

not  be  rejected  at 

the  .05  significance 

level. 

TABLE  2 

CHI  SQUARE  TEST  ON  5000  NUMBERS 

Number  of  Intervals 

n-19 

X2 

n«25 

X2 

MULTIPLICATIVE 

5.46  .5002 

21.83 

,5002 

MIXED 

10.03  .5012 

27.70 

.5012 

ADDITIVE 

14.61  .4987 

35.03 

.4987 

TABLE  Chi  Sq 

16.92 

36.42 

a* .  65 

Table  3  summarizes  the  results  for  the  local  behavior 
of  each  generator.  The  multiplicative  and  aixed  generators 
pass  all  local  tests  but  the  additive  generator,  for  a 
saaple  size  of  500,  does  not  pass  one  of  the  Chi  Square 
tests.  Such  behavior  should  be  noted  but  is  not  of  itself 
due  cause  for  rejection. 

The  means  for  the  saaple  size  of  500  were  subjected  to 
a  Sign  Test.  Consider  the  aeans  as  a  randoa  saaple  froa  a 
binomial  population  with  probability  .5  of  obtaining  a 
value  greater  than  the  aean  of  a  uniform  distribution  on 
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TABLE  3,  CHI  SQUARE  TESTS  OF  LOCAL  BEHAVIOR 
(9  degrees  of  freedon) 

MULTIPLICATIVE  GENERATOR 


SAMPLE  SIZE 

:  500 

1000 

TEST  NO. 

X2 

MEAN 

X2 

MEAN 

1 

4.32 

.5067 

8.36 

.5129 

2 

9.28 

.5191 

9.12 

.5088 

3 

12.68 

.5147 

4,46 

.  3  335 

4 

7.28 

.  5028 

7.84 

.4^95 

5 

14.20 

.5170 

9.28 

.4864 

6 

7.76 

.4899 

7 

6.64 

.5009 

8 

15.32 

.4781 

9 

8.88 

.4856 

10 

12.96 

.4871 

MIXED  GENERATOR 

1 

8.88 

.4994 

6.84 

.5046 

2 

6.52 

.5099 

12.44 

.4923 

3 

13.20 

.4971 

6.30 

.4936 

4 

8.20 

.4875 

11.40 

.5080 

5 

6.76 

.4912 

9.12 

.5074 

6 

3.12 

.  4959 

7 

13.80 

.4  869 

8 

10.16 

.5292 

9 

10.48 

.5093 

10 

4.48 

.  5054 

ADDITIVE  GENERATOR 

1 

4.76 

.4928 

11.64 

.4977 

2 

17.40 

.5026 

10.32 

.4997 

3 

7.08 

.4997 

11.02 

.4964 

4 

9.52 

.4997 

8.00 

.5097 

5 

13.52 

.4860 

8.20 

.4901 

6 

8.52 

.5069 

7 

5.16 

.5129 

8 

8.52 

.5065 

9 

6.16 

.4922 

10 

5.36 

.4880 

TABLED  CHI 

SQ  16.919 

■ .  05 
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.5  and 


(0,1).  The  null  hypothesis  was  taken  as  H^:  p  • 
the  alternative  H^:  p  /  .5.  The  critical  region  is  for 
values  more  than  8  or  less  than  2  above  the  mean.  The 
null  hypothesis  cannot  be  rejected  for  any  generator. 

The  means  could  be  froa  a  uni fora  population. 

KolBogorov-Smirnov  Test 

The  second  test  of  uniformity  is  the  Kolmogorov* 

Smirnov  Goodness  of  Fit  Test  (KS  test).  The  test  is 
based  on  the  difference  between  the  theoretical  cumula- 
tive  distribution  function  and  the  empirical  distribution 
function  which  is  derived  from  the  saapled  values.  The 
KS  test  may  be  used  when  the  theoretical  cumulative  dis¬ 
tribution  function  is  continuous.  For  a  detailed  dis¬ 
cussion  the  reader  is  referred  to  Knuth  119, p. 41]  or  a 
statistics  teat. 

The  test  was  first  applied  to  5000  numbers  in  100 
blocks  of  50  each  and  then  again  applied  to  the  100  suc¬ 
ceeding  results.  The  null  hypothesis  is  that  the  numbers 
are  uniforaly  distributed.  The  critical  region  is  K  >  Kn ,<^/a 
and  K<  K^,  where  KR,  n  ■  100  are  tabled  values 

[19,p.44].  Global  behavior  is  described  by  the  block 
statistics  given  in  Table  4.  The  null  hypothesis  can¬ 
not  be  rejected  for  any  generator.  For  local  behavior 
in  samples  of  size  50,  all  generators  failed  soae  tests. 
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The  multiplicative  failed  thirteen  out  of  100,  and  the 
mixed  and  additive  nine  each. 


TABLE  4  KS  GOODNESS  OF  FIT  TEST 
on  5000  numbers  in  100  blocks  of  50  each 


GENERATOR: 

MULTIPLICATIVE 

MIXED 

ADDITIVE 


TABLED  K  VALUE 


1.0651  .5020 

.7899  .3513 

.6054  .1591 


97.51  2.5% 

.0961  1.3879 


.7287  .1948 

.7683  .1999 

1.0791  .8621 


Runs  Above  and  Below  the  Mean 

A  run  is  a  succession  of  identical  symbols  which  is 
followed  and  preceded  by  different  symbols.  For  a  sequence 
of  random  numbers  r2,  . ..,  r  ,  there  corresponds  a  se¬ 
quence  sj,  . . . sn  where  S£  is  the  letter  "a"  if  r ^ 

is  greater  than  .5  and  "b"  if  less  than  .5.  The  resulting 
sequence  might  look  like  a  a  b  a  b  b  a... etc. 

Assuming  independence  between  the  random  numbers,  the 
probability  and  expected  numbers  of  runs  of  various  lengths 
can  be  calculated  and  a  Chi  Square  statistic  formed.  The 
null  hypothesis  is 
critical  region  is 


that  the  numbers  are  random  and  the 
X* 


TABLE  5  RUNS  ABOVE  AND  BELOW  MEAN 

GENERATOR  CHI  SQUARE  STATISTIC  (8  deg.  freedom) 

MULTIPLICATIVE  8.423 

MIXED  6.900 

ADDITIVE  11.208 

TABLED  CHI  SQ  15.507  *-.05 

The  null  hypothesis  cannot  be  rejected  for  any  generator. 

Correlation  Coefficient 

If  a  set  of  numbers  is  randoa,  a  given  nuaber  should 
not  depend  on  its  predecessors.  The  correlation  for  any 
given  lag  L  should  be  zero.  Since  we  do  not  have  truly 
randoa  nuabers,  the  aagnitude  of  any  correlation  should 
be  snail  and  have  equally  distributed  positive  and  nega- 
tive  signs. 

A  correlation  test  for  lags  froa  1  to  15  was  nade  on 

a  set  of  4985  numbers  froa  each  generator.  The  results 

are  displayed  in  Table  6.  Of  special  interest  is  the 

correlation  of  lag  1  (serial  correlation).  Anderson  [1] 

has  shown  that  the  serial  correlation  for  large  samples 

(N  >  75)  is  approximately  normally  distributed  with 

2 

mean  -1/(N-1)  and  variance  (N-2)/(N-l)  .  The  single  tail 
significance  points  are  -1  i  1.645  V  N-2  i.e.,  for 

- n - * 

N«5000j -0.0235  <  Pi  <  +  0.0231.  No  generator  is  reject¬ 

ed  on  the  basis  of  the  serial  correlation. 
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TABLE  6 

CORRELATION  TEST  OF  4985  NUMBERS 


LAG 

MULTIPLICATIVE 

MIXED 

ADDITIVE 

RANGE 

1 

-.0101 

.0218 

.0024 

-.0176 

2 

.0053 

.0043 

-.0008 

-.0113 

3 

.0028 

.0012 

-.0021 

.0063 

4 

-.0027 

.0076 

.0130 

.0155 

5 

-.0*98 

.0143 

.0044 

.0016 

6 

-.0127 

-.00004 

-.0226 

.0065 

7 

.0089 

-.0030 

.0053 

-.0079 

8 

-.0003 

-.0177 

-.0105 

-.0162 

9 

.UUO/ 

.0280 

-.0040 

.0042 

10 

.0018 

.0191 

-.0023 

.0007 

11 

-.0044 

.0091 

-.0073 

-.0083 

12 

.0033 

-.0054 

.0023 

-.0004 

13 

10149 

.0116 

-.0055 

-.0030 

14 

.0001 

-.0061 

.0130 

-.0161 

15 

.0093 

-.0011 

-.0114 

.0295 

A  sign  test  can  be  applied  to  the  other  correlations  of 
lag  L.  Consider  the  signs  of  the  correlations  as  a  ran¬ 
dom  sample  from  a  binomial  pvpulation  with  probability 
.5  of  obtaining  a  positive  correlation  The  null  hy¬ 
pothesis  is  Hq:  p  ■  .5  and  the  alternative  H^.  p  t  .5. 

For  sample  size  N-15  and  0<  ■  .05,  the  null  hypothesis 
is  rejected  if  the  number  of  successes  (positive  signs) 

1?  or  less  than  3.  The  null  hypothesis  cannot 
be  rejected  for  any  generator. 
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Runs  Up 

A  staple  of  5000  numbers  from  etch  generttor  wts 
subjected  to  Knuth's  Runs  Up  test  |T9  P.60].  The  staple 
is  examined  for  lengths  of  sequences  which  tre  incretsing 
and  t  statistic  coaputed  which  has  a  Chi  Square  distri¬ 
bution.  The  critical  region  is  X*  >  'jt  .05,  Y  . 

TABLE  7.  RUNS  uP 

GENERATOR  CHI  SQUARE  STATISTIC  (6  deg.  freedoa) 

MULTIPLICATIVE  6.678 

MIXED  2.998 

ADDITIVE  4.210 

TABLED  CHI  SQ  12.59  •  .05 

The  null  hypothesis  cannot  be  rejected  fme  any  generator. 
Poker  Test 

The  Poker  Test  is  a  frequency  test  for  conbinations 

of  distinct  values  in  a  set  of  five.  The  categories  are: 

5  different  *  all  different 

4  different  =  one  pair 

3  different  *  two  pairs  or  three  of  a  kind 

2  different  *  full  house  or  four  of  a  kind 

1  different  *  five  of  a  kind 

Counts  in  each  category  are  compared  to  the  expected  numbers 

and  a  Chi  Square  test  is  made. 

The  Poker  Test  was  applied  to  1000  groups  of  five 

numbers  to  test  independence.  Table  8  below  shows  the 

Chi  Square  statistics.  The  null  hypothhsis  is  that  the 

numbers  are  independently  distributed  against  the  alternative 

that  they  are  not.  The  critical  region  is  *X2  £  .OS^Y  . 

57 


TABLE  8  ,  POKER  TEST 


GENERATOR 


CHI  SQUARE  STATISTIC  (3  deg.  freedom) 


MULTIPLICATIVE 

MIXED 

ADDITIVE 


3.439 
.937 
.  432 


TABLED  CHI  SQ 


.815 


*  -.05 


The  null  hypothesis  cannot  be  rejected  for  any  generator. 
Coupon  Collector  Test 

The  Coupon  Collector  Test  was  applied  to  5000  num¬ 
bers  for  a  set  of  integers  0  to  4.  The  lengths  of  sequences 
to  complete  the  set  of  integers  is  observed  and  a  Chi 
Square  statistic  calculated.  The  null  hypothesis  is  that 
the  numbers  are  independently  distributed.  The  critical 
region  is  X2  >  'jCVes,  11. 

IABLE  9  COUPON  COLLECTOR  TEST 
GENERATOR  CHI  SQUARE  STATISTIC  (11  deg.  freedom) 


MULTIPLICATIVE 

MIXED 

ADDITIVE 


10.123 

20.509 

9.656 


TABLED  CHI  SQ  19.675  o(  »  .05 

The  mixed  generator  failed  the  test,  but  the  null  hypothesis 
was  not  rejected  for  the  other  generators. 


Permutation  Test 

If  the  numbers  are  truly  random,  there  should  be  an 
equal  probability  of  any  permutation  of  digits.  The  re- 
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suits  of  the  test  are  given  in  Table  10.  The  null 
hypothesis  is  that  the  numbers  are  independently  distri¬ 
buted.  The  critical  region  is  X2  >  ^*.05,  Y  . 

TABLE  10.  PERMUTATION  TEST 

GENERATOR  CHI  SQUARE  STATISTIC  (119  deg.  freedom) 


MULTIPLICATIVE  122.48 
MIXED  120.08 
ADDITIVE  103.28 


TABLED  CHI  SQ  157.7 


The  null  hypothesis  is  not  rejected  for  any  generator. 
Spectral  Test 

The  spectral  test  fl9Qwas  applied  to  the  values  pro¬ 
duced  by  the  nultiplicative  and  nixed  congruential  gen¬ 
erators.  It  is  a  test  of  independence  among  n-tuples  of 
numbers  produced  by  linear  congruential  generator*  and 
successful  performance  depends  on  the  choice  of  the  multi 
plier  and  modulus.  The  results  for  values  of  n  -  2,3,4 
are  given  in  Table  11. 

TABLE  11  SPECTRAL  TEST 

GENERATOR  C(2)  C(3)  C(4) 

MULTIPLICATIVE  1.57  c  .125  x  10"5  .193  x  10’ J 

MIXED  .152  x  10’5  .262  x  10"3  .398  x  10’1 

PASSING  .1  .1  .1 

UPPER  BOUND  3.63  5.90  9.86 
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Knuth  [  I9>tates  that  a  generator  passes  the  spectral 
test  if  C(n)  *  9.1  for  n  ■  2,  3  and  4,  and  C(n)  >  1  is 
considered  very  good.  Based  on  the  criteria  neither 
generator  passes  the  spectral  test.  The  nultiplicative 
generator  has  a  good  value  of  C(2)  indicating  independ¬ 
ence  of  consecutive  pairs,  but  the  values  of  C(3)  and 
Ci.4)  are  too  low  as  are  all  the  values  for  the  nixed 
generator.  The  test  results  suggest  that  to  improve  in¬ 
dependence  a  larger  multiplier  should  be  used.  Knuth  [19] 
has  found  that  a  multiplicative  congruential  generator  with 
a  multiplier  of  515  has  very  good  results  on  the  spectral 
test. 

Serial  Test 

Pairs  of  successive  random  numbers  should  be  distri¬ 
buted  uniformly  and  independently.  In  the  serial  test, 
the  number  of  times  the  pair 

(r2i»  r2i  +  l)  *  (a»  b)  for  0  j  <  n 
J  **  1  04  v  b4d 

occurs  is  counted.  A  Chi  Square  test  is  applied  to  these 
2  2 

d  categories  with  a  1/d  probability  for  each  category. 

2 

The  expected  number  in  each  category  is  n/d  .  The  serial 

test  was  applied  to  a  sequence  of  5000  numbers  for  pairs 

whose  digits  are  between  0  and  9,  and  a  Chi  Square  statistic  is 

2  2 

computed.  The  null  hypothesis  is  that  X*<  .05  ^  V  . 

The  results  are  displayed  in  Table  12. 
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TABLE  12 
SERIAL  TEST 


GENERATOR  CHI  SQUARE  STATISTIC  (99  deg.  freedom) 


MULTIPLICATIVE  91.46 
MIXED  125.86 
ADDITIVE  83.30 
RANGE  120.98 


TABLED  CHI  SQ  134.51 

The  null  hypothesis  cannot  be  rejected  tor  any  generator. 


Gap  Test 

The  Gap  Test  incorporates  both  independence  and 
uniformity  in  examining  the  lengths  of  subsequences  (gaps) 
separating  values  in  a  specific  range.  Let  <X,  be 
such  that  0  £  c<  4  ^  £  1.  Consider  lengths  of 

subsequences  separating  specific  values  Xj.^,  Xj  ,  xj  +  i* 
....  Xj+r,  so  that  Xj.j,  xj+r  belong  to 

j  °(  4  Xj_lt  Xj  +  r  |  Xj+k  <  <*}  ^|xj+k  >/$  j  1 

where  k»0,  1. . . ,  r-1 

The  r+2  numbers  give  a  gap  of  length  r.  Suppose  the  in¬ 
terval  is  .4  to  .5.  Then  the  subsequence 
X-i  X0  Xj  x2  x3  x4  x5 

.41  .32  .22  .78  .61  .10  .48 

not  in  the  interval 
represents  a  gap  of  length  5. 

The  distribution  of  gap  lengths  is  geometric 

h(j)  -  fp(l-p)^*1  0  £  n 

)  0  otherwise 
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where  p  ■(&-©(  ,  the  probability  that 

<*  £  %«.  i  £ 

A  Chi  Square  test  is  applied  to  the  t+1  values  of  counts 
of  length  r,  Of  r  it,  The  gap  test  has  aspects  of  uniformity 
and  independence:  i .e.  /’uniformity"  (or  equally  likely) 
because  the  same  probability  p  is  assigned  to  the  occur¬ 
rence  of  in  (®<j£  ),  and  "independence"  because 
independent  geometric  trials  are  assumed.  For  each 
generator  from  among  5000  numbers,  500  gaps  were  requested 
for  the  intervals  . 4  to  .5  and  .5  to  .6.  Table  13  presents 
the  results.  The  null  hypothesis  is  that  X^  <T:  flS,Y 
The  mixed  generator  failed  both  gap  tests  while  the  null 
hypothesis  cannot  be  rejected  for  the  other  generators. 

TABLE  13 
GAP  TESTS 

500  gaps  on  the  intervals 
. 4  to  . 5  . 5  to  . 6 


GENERATOR 

DEG 

FREEDOM 

X2 

DEG  FREEDOM 

MULTIPLICATIVE 

20 

25.505 

20 

29.506 

MIXLD 

19 

35.969 

*  ~ 

33.321 

ADDITIVE 

20 

21.036 

20 

28.438 

TABLED  CHI 
TABLED  CHI 

SQ  (19 
SQ  (20 

Deg  f) 
Deg  f) 

30.144 

31.410 

<X  :  .05 

Sensitivity  to  Initial  Values 

The  multiplicative  and  mixed  congruential  generators 
were  testeu  with  three  different  starting  values  with 
respect  to  the  Chi  Square,  Kolomogorov-Smirnov,  and  Runs 
Above  and  Below  the  Mean  tests.  The  Chi  Square  was  a 


62 


test  of  global  behavior  on  all  5000  embers.  The  nuiiber 
corresponding  to  KS (local)  is  a  count  of  the  number  of 
failures  of  the  KS  test  in  100  repetitions  of  sample  size 
50.  Table  14  summarizes  the  results. 


TABLE  14  INITIAL  VALUE  TEST 


MULTIPLICATIVE 


starting  value 

1  203  1307 


MEAN 

CHI  SQUARE 
KS  (global) 
KS  (local) 
RUNS 


4976 

.5016 

.5002 

N 

N 

N 

N 

R 

N 

12 

11 

9 

N 

N 

N 

**  NOTE:  N  •  not  reject  null  hypothesis;  R  ■  reject 


The  initial  value  does  have  some  effect  on  the  properties 
of  the  sequence  of  numbers  generated.  A  very  small  prime, 
say  one,  produces  detrimental  effects,  such  as  failure 
to  pass  the  global  KS  test  for  a  uniform  distribution. 
Larger  initial  values  moved  the  mean  closer  to  .50. 

The  additive  generator  was  tested  with  initial  values 
generated  by  a  multiplicative  with  its  initial  value  equal 
to  one.  The  additive  generator  failed  the  KS  test  when 
applied  to  the  block  statistics.  The  additive  type  of 
all  generators  should  be  most  sensitive  to  starting  values. 


Generation  Time  and  Storage 

All  pseudorandom  number  generator  subprograms  were 
timed  on  a  UNIVAC  1108  for  the  generation  of  5000  numbers. 
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The  mean  times  for  five  runs  each  are  displayed  in  Table  15. 
The  fastest  are  the  multiplicative  and  mixed  while  the 

additive  is  substantially  slower. 

» 

TABLE  15  .  MEAN  GENERATION  TIME 


GENERATOR  SECONDS 

MULTIPLICATIVE  .07 
MIXED  .07 
ADDITIVE  (79)  .23 
ADDITIVE  (16)  .22 


There  is  no  significant  difference  between  storage  require¬ 
ments  of  the  multiplicative  and  mixed  generators.  Both 
depend  only  on  the  last  number  generated  and  have  no  require 
ment  for  additional  values;  also  the  programs  themselves 
are  the  same  lenyth. 

However,  the  additive  generator  requires  at  least  16 
initial  values  to  start,  which  also  must  be  updated.  For 
example,  the  additive  congruential  generator  used  in 
this  report  had  79  initial  values.  Such  requirements  could 
be  a  hindrance  if  storage  is  a  severe  restriction  as  in  a 
mini- computer. 

Conclusion  Based  on  the  Statistical  Results 

The  multiplicative  congruential  generator  produces 
the  most  favorable  statistical  results,  has  the  fastest 
generation  time,  and  has  the  least  storage  requirements 
of  all  generators  considered.  The  multiplicative  did 
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have  a  higher  number  of  failures  of  the  Kolmogorov- 
Smirnov  test  for  local  uniform  distribution  behavior 
than  the  other  generators.  The  difference  did  not  ap¬ 
pear  significant.  Unsatisfactory  local  behavior  is  pos¬ 
sible  even  in  a  sequence  with  good  global  behavior  [  15]. 
The  generator  did  fail  the  spectral  test  overall,  but 
independence  for  consecutive  pairs  was  good. 

The  mixed  congruential  generator  has  a  theoretical 
advantage  of  longer  period  and  an  easier  basic  theorem 
to  satisfy.  However,  its  statistical  behavior  is  not  as 
good.  In  particular,  the  mixed  generator  did  not  pass 
the  spectral,  coupon,  and  gap  tests.  Time  and  storage 
requirements  are  equivalent  to  those  of  the  multiplica¬ 
tive.  The  mixed  generator  as  formulated  for  this  report 
is  not  recommended  for  general  use.  The  failure  of  the 
spectral  test  indicates  a  conspicuous  lack  of  independence 
for  consecutive  pairs,  triples,  and  quadruples. 

The  additive  type  is  markedly  slower  but  has  favor¬ 
able  statistical  behavior.  It  appears  that  these  genera¬ 
tors  are  sensitive  to  initial  starting  values.  Storage 
requirements  are  the  major  disadvantage  for  additive 
generators  with  acceptable  random  numbers. 

A  sumary  of  results  for  selected  tests  on  all 
generators  is  given  in  Table  16,  xabie  17  also  pre¬ 
sents  the  probabilities  P(  X*  £  x*)  for  a  better 
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comparison  within  a  given  set  of  test  results.  The 
serial  correlation  is  the  actual  computed  value. 

Among  the  three  congruential  generators,  the  multiplica* 
tive  has  the  advantage  in  uniformly  distributed  numbers. 
For  independence  the  results  are  close  between  the  addi¬ 
tive  and  multiplicative.  Serial  correlation  is  lowest 
for  the  additive,  possibly  giving  it  the  edge  in  con¬ 
siderations  of  independence. 


TABLE  16  SUMMARY  @F  STATISTICAL  TESTS 


MULTIPLICATIVE  MIXED  ADDITIVE 


Uniformly  Distrib. 


CHI  SQUARE 
KS 


N  N  N 

N  N  N 


Independence 


POKER  N 
COUPON  N 
PERMUTATION  N 
RUNS  ABOVE  §  BELOW  MEAN  N 
SERIAL  CORRELATION  N 
RUNS  UP  N 
SPECTRAL  R 


N 

R 

N 

N 

N 

N 

R 


N 

N 

N 

N 

N 

N 


SERIAL 

GAP 


N  N  N 

N  R  N 


TESTS 

REJECTED  1  30 


Note:  N  ■  Null  hypothesis  not  rejected;  R»  rejected. 
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TABLE  17 


COMPARISONS  OF  TESTS 


MULTIPLICATIVE 

MIXED 

ADDITIVE 

RANGE 

Uniformly  Distrib 

\  CHI  SQUARE 

20 

70 

90* 

86 

l  POKER 

67* 

20 

8 

89 

|  COUPON 

SO* 

96* 

43 

33 

Independence 

RUNS  ABOVE  6  BE¬ 

LOW  MEAN 

60 

45 

80* 

65 

RUNS  UP 

63* 

20 

37 

49 

SERIAL  COR¬ 

RELATION 

-.0101 

.0218 

.0024* 

-.0176 

GAP  (.4  to  .5) 

80* 

98X 

60 

50 

Note:  *  The  highest  non-failing  probability  for  that  test 
(disregarding  RANGE) 
x  Failed  at  .05  significance  level 


A  Recommended  Pseudorandom  Number  Generator 

The  above  results  indicate  that  a  multiplicative 
con^ruential  generator  produces  the  most  favorable  results. 
For  computers  having  a  word  length  of  2^,  the  generator 
with  initial  value  1907,  and  multiplier  515  demonstrates 
very  good  statistical  behavior.  The  FORTRAN  IV  subroutine 
is  Figure  13. 

Random  Variates 

The  pseudorandom  number  generators  produced  numbers 
between  zero  and  m-1  in  modulo  m.  A  keen  observer  would 
have  noticed  that  in  the  FORTRAN  subroutines  the  random 
number  was  divided  by  the  word  length  2^  to  produce  a 
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C  *  *  *  PSEUDORANDOM  NUMBER  GENERATOR  *  *  * 

C 

C  FUNCTION  RANG (DUMMY ) 

I  -  IA  *  I 

FLD  (  0,1,1)  «  0  .  _ 

RANG  -  1/2.0  **  35 
RETURN 
C 

C  INITIALIZATION  OF  MULTIPLIER  AND  INITIAL  VALUE 

C 

ENTRY  RAN  (IN1TAL) 

I  -  INITAL 
IA  »  j  **  15 
RETURN 
END 


Figure  13.  FORTRAN  Subroutine  for  Pseudorandom  Number 
Generator 
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Mati'Miir  'I*."  - 


uniform  random  variable  on  (0,1).  To  generate  other 
random  variates  we  must  consider  some  transformation 
techniques  from  the  uniform  to  other  distributions. 

Transformations  are  broadly  categorized  as  table 
look  up  routines  or  mathematical  techniques.  Many  talle 
look  ups  are  based  on  an  inverse  transformation.  Suppose 
F (X)  is  the  cumulative  distribution  function  correspond¬ 
ing  to  a  desired  probability  density  function  f(X). 

From  probability  theory  we  know  that  the  probability  dis¬ 
tribution  for  the  cumulative  distribution  function  F(X) 
is  uniform  on  (0,1).  Therefore  we  can  use  the  uniformly 
distributed  random  numbers  to  generate  the  distribution 
of  F  (X)  regardless  of  the  functional  form  of  distri¬ 
bution.  It  follows  then  that  for  r  a  random  variable  on 

(0,1)  r  -  F (X) 

and  . 

X  -  F'-^r) 

F  (r)  is  the  inverse  of  the  original  cumulative  distri¬ 
bution.  In  the  computer,  an  array  would  represent  the 
functional  values  of  F(X)  and  a  subroutine  would  inter¬ 
polate  to  determine  the  value  of  X  corresponding  to  the 
F(X)  produced  by  the  pseudorandom  number  generator. 

Mathematical  Techniques 

The  inverse  transformation  can  also  be  applied 
mathematically  to  simple  distributions  such  as  the  ex¬ 
ponential. 
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F  (Z)  -  1  -  e-2/0 

--  "  In  (l-F(Z) ] 

e 

Z  -  -91n(l-f)  let  r  -  F(z) 

or  using  symmetry, 

Z  ■  -0  In  r,  where  r  is  the  uniform  random  number. 

C  *  *  *  FORTRAN  FUNCTION  FOR  EXPONENTIAL  *  *  * 

C  *  *  *  RANDOM  VARIATES 
C 

FUNCTION  EXPRV  (XMEAN) 

EXPRV  «  -  XMEAN  *  ALOG  (RANG (RANUM) ) 

RETURN 
END 

Figure  14.  A  FORTRAN  Function  for  Exponential  Random 
Variates . 

Geometric 

The  probability  density  function  for  the  geometric 
distribution  is  f(x)  *  pqx’*,  X  *  1,  2,... 
where  p  ■  1-q 

Then  FCX)  .£  pq1"1 

y-l  70 


» 


1-F  (X)  -  Z*  ,  .y-1  "  £  „  .  v.! 

y-X+l^^  y»X+l  1 


i-f(x: 

q 


-  £ 

y«X+l  (qy  ^-qX)  ■  q* 


-  qX_1 


X-l  -  log^r 

X-l  -  In  r 
In  q 

X  -  In  r  ♦  1 
In  q 

where  r  is  the  uniform  random  number  and  q  is  the  prob¬ 
ability  of  failure  on  an  individual  trial. 


C  *  *  *  FORTRAN  FUNCTION  FOR  *  *  * 

C  *  *  *  GEOMETRIC  RANDOM  VARIATES  *  *  * 


FUNCTION  GEOM  (Q) 

GEOM  -  ALOG  (RANG (RNUM) ) /ALOG (Q) ♦ 1 


RETURN 


Figure  15  A  FORTRAN  Function  for  Geometric  Random  Variates 


Negative  Binomial 


A  negative  binomial  distribution  can  be  described 
as  a  series  of  Bernoulli  trials  which  are  repeated  until 
K  successes  have  occurred.  The  pdf  is 
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X-0,  1,  2... 


(where  probability  of  success  i;  p)  which  is  the  sun  of 

K  geometric  variates, 
k 

X  ■  In  r. 

i-1 


X  - 


*n  lit. 

In  q 


ri  is  random  number  i»l,...,  k 
q  is  probability  of  failure 


C  *  *  *  FORTRAN  FUNCTION  FOR  NEGATIVE  *  *  * 

C  *  *  *  BINOMIAL  RANDOM  VARIATES  *  *  * 

C 

FUNCTION  NEGBN  (K,  Q) 

PROD  -  1.0 
DO  10  I  -  1,  K 
10  PROD  -  PROD  *  RANG  (RANUM) 

NEGBN  -  ALOG  (PROD)/ALOG  (Q) 

RETURN 

END 

Figure  16  .  A  FORTRAN  function  for  Negative  Binomial 
random  variates. 


Uniform  on  (a,  b) 

The  pseudorandom  number  generator  produces  uniform 
random  variates  on  (0,  1).  But  frequently  uniform  ran¬ 
dom  variates  are  required  on  other  intervals,  say  (a,  b)\ 
The  uniform  density  function  is 


£(x)  - 


a  <  x  <  b 
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Otherwise 


The  cumulative  distribution  function  is 


F  (x) 


r 


1 


dt 


x-a 


a 


Let 

so 


F(x)  -  r 

s-a  »  r 
F=* 


and  x  ■  a  ♦  (b  -  a)  r 

where  r  is  a  uniform  random  variate  m  (0,  1) 


C  *  *  *  BORTRAN  FUNCTION  FOR  UNIFORM  *  *  * 

C  *  *  *  RANDOM  VARIATES  ON  (A,  B)  *  *  * 

C 

FUNCTION  UNIFM  (A,  B) 

UNIFM  ■  A  *  (B  »  A)  *  RANG  (RANbM) 

RETURN 

END 

Figure  IT  A  FORTRAN  Function  for  uniform  Random  variates 
on  (a,  b) 


Rejection  Method 

The  rejection  method  can  be  applied  to  any  probability 
distribution  for  which  upper  and  lower  bounds  can  be 
placed  on  the  range  of  values.  For  example,  consider  the 
binomial  distribution. 


Binomial 


The  binomial  can  be  described  as  X  successes  in  n 
trials  with  probability  P  of  success: 
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f(x)  ■  £)  pX  qn’X  X  -  0,  1,  2,  ...  n 
The  following  algorithm  effectively  reproduces  the  n 
bernoulli  trials  and  rejects,  i.e.,  does  not  count,  the 
failures. 

1.  I  »  0 
KOUNT  »  0 

2.  I  -  I  +  1 

3.  IF  I  >  n  ,  STOP 

4.  GENERATE  UNIFORM  RANDOM  NUMBER  RI 

5.  IF  RI  ^  P,  KOUNT  »  KOUNT  +  1 

6.  GO  TO  2 


The  value  of  KOUNT  is  the  number  of  successes  in  n  trials. 


Poisson 


Generating  Poisson  random  variates  depends  on  the 
known  relationship  between  the  exponential  and  Poisson 
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distribution.  The  Pdf  for  the  Poisson  is 


f(x) 


'V 


x  -  0,  1,  2. . . 

o 


If  the  tine  interval  between  events  has  an  exponential 
distribution  with  nean  1/7*  ,  then  the  nunber  of  events 

x  occurring  during  a  unit  tine  interval  has  a  Poisson 
distribution  with  nean  X  .  Thus ,  the  nethod  is  to 
generate  exponential  randon  variates  t^,  t^ ,  ...  where 

t.  *  -  In  r.  and  accumulates  then  until^  .  *1* 

1  1  *  li 4  * 4  £ 

i-1  i.i 

the  sum  exceeds  ^  .  Then  X  the  number  of  exponential 

variates  is  the  desired  Poisson  random  variate. 


Gamma 

The  Gamma  distribution  is  also  related  to  the  ex¬ 
ponential.  The  pdf  for  the  gamma  is 
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C  *  *  *  FORTRAN  FUNCTION  FOR  GAMMA  *  *  * 

C  *  *  *  RANDOM  VARIATliS  *  *  * 

C 

FUNCTION  GAMMA  (IALPHA,  BETA) 

PROD*  0. 

DO  10  I  =1,  IALPHA 
10  PROD  =  PROD  *  RANG  (RANUM) 

GAMMA  *  -  BUT A  *  ALOG  (PROD) 

RLTURN 

LND 

Figure  2(  A  FORTRAN  Function  for  Gamma  Random  Variates. 


Normal 

Since  the  normal  distribution  is  so  frequently  used, 
several  techniques  have  been  developed  to  generate  normal 
random  variates.  Muller  [22]  discusses  six  approaches  in 
current  use.  Three  methods  will  be  discussed  below. 


The  most  common  approach  is  table  look  up.  In  GPSS 
language,  for  example,  values  of  the  cumulative  distribution 


C  *  *  *  FORTRAN  FUNCTION  FOR  NORMAL  *  *  * 

C  *  *  *  RANDOM  VARIATES  *  *  * 

C 

FUNCTION  XNORMAL  (MEAN,  STD) 

REAL  MEAN 
SUM  -  0. 

DO  10  I  -  1,  12  » 

10  SUM  -  SUM  ♦  RANG  (RANUM) 

XNORMAL  -  STD  *  (SUM  -  6.)  ♦  MEAN 

RETURN 

END 

Figure  21  A  FORTRAN  Function  for  Normal  Random  Variates 
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function  for  the  particular  normal  distribution  are 
stipulated  as  a  table  which  is  then  searched  to  determine 
the  random  variate. 

Box  and  Muller  f  4  ]  have  developed  a  direct  approach 
to  generate  a  pair  of  random  variates  from  the  S4\me  nor¬ 
mal  distribution.  Let  r^  and  r2  be  uniformly  distributed 
independent  random  variables  on  (0,  1).  Then, 

l/2 

Xj  ■  (-2  In  r^)  cos  2  r 

v  1/2  2 
X2  ■  (-2  In  rj)  sin  2  r 2 

Where  (X^,  X2)  will  be  a  pair  of  independent  random  variables 
whose  joint  pdf  is  a  standard  normal  distribution. 

The  last  technique  is  based  on  the  Central  Limit 
Theorem  which  states  that  the  mean  of  identical  uniform 
independent  random  variables  is  asymptotically  normally 
distributed.  Consider  a  set  r^  of  K  uniformly  distributed 
independent  variables  on  (0,  1).  Then 


x  .  p==== 

y  k/12 

is  normally  distributed  with  mean  0  and  variance  1.  Thus, 
to  generate  a  normal  variate  with  mean  and  standard 

deviation  <T 

*  *!-!>♦/* 

i-l 

The  values  of  X  are  reliable  within  three  standard  devia¬ 
tions.  The  value  of  K  should  be  no  less  than  10 » and  12  is 
frequently  used. 
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SECTION  IV 


SYSTEM  REPRESENTATION 

A  model  is  a  representation  of  a  real  world  system 
which  the  analyst  uses  for  prediction  or  control.  The 
model  should  reasonably  approximate  the  real  world  system 
and  incorporate  the  important  features  of  the  real  system 
but  not  become  so  complex  that  it  is  impossible  to 
understand  cr  manipulate. 

In  this  section  we  will  look  in  more  detail  at  model 
building  and  the  components  of  the  model  which  represent 
the  real  system.  It  is  important  to  realize  that  model 
building  is  as  much  an  art  as  a  science  and  that  success 
depends  on  the  analyst's  experience  as  well  as  technique 
and  considerable  luck. 

Model  Structure 

A  basic  problem  is  how  to  describe  the  system  using 
a  simulation  model.  For  our  purposes,  we  will  define 
four  elements  which  comprise  .the  model  structure: 
entities,  attributes,  functional  relationships,  and  a 
time  flow  mechanism. 

Entities  are  the  objects  by  which  a  system  can  be 
defined,  i.e.,  the  components  of  „he  system  represented 
in  the  model.  They  are  further  divided  into  permanent 
and  temporary  entities.  Permanent  entities  are  entities 
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whose  presence  in  the  model  are  unaffected  by  time. 
Temporary  entities  are  entities  whose  presence  in  the 
model  is  affected  by  time. 

Using  set  notation,  the  entities  are  described  as 


|  Pi  *?2 *  *  *  *  *  ^n*  tj,*  *  2 '  *  *  *  *m  j * 


The  permanent  entities  form  a  subset  P  •  jpj,  P2,...,Pn| 
PS  E  and  similarly,  the  temporary  entities  form  a  subset 
T  ■  jtl»  t2»***»  tm|*  T£  E.  As  a  simple  example  consider 
a  queuing  system  consisting  of  a  single  server  and  single 


CUSTOMERS 


SERVED 


CUSTOMERS 


waiting  line.  The  permanent  entity  is  the  single  server- - 
he  is  always  there  whether  customers  are  there  or  not. 

The  temporary  entities  are  the  customers  which  arrive, 
wait,  are  served,  and  then  leave  the  system. 

The  next  element  of  the  model  is  the  set  of  attributes. 
Attributes  are  parameters  or  variables  associated  with 
the  entities.  Attributes  may  be  further  classified  as 
exogenous,  endogenous,  or  status  variables.  Exogenous 
variables  are  the  independent  or  input  variables  which 
are  assumed  to  be  predetermined  and  given  independent  of  the 
system  modeled.  They  have  values  which  affect,  but  are 
unaffected  by  the  system.  Endogenous  variables  are  the 
dependent  variables  of  the  system.  Their  values  are 
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determined  by  the  other  variables  in  the  system. 


The  attribute  set  is 
A 


( 


1 


v(ej)  :  e  j  6  E 
where  3j  is  some  permanent  or  temporary  entity  in 
the  model. 

Functional  relationships  are  the  expressed 
relationships  among  permanent  and  temporary  entities. 
They  describe  the  interaction  of  the  variables  and 
attributes  of  the  model. 

The  functional  relationship  set  is 
E  “  I  fi(Pj  t  t^)  :  jj,  6  P,tkeT 


! 


for  some  i. 

Returning  to  the  queuing  example ,  the  exogenous  variables 
are  the  interarrival  time  of  the  customers  and  the 
interservice  time  of  the  server.  The  status  variables 
include  the  amount  of  time  a  customer  spends  waiting  and 
the  amount  of  time  the  server  is  idl*.  An  endogenous 
variable  is  the  total  time  a  customer  is  in  the  system. 
Parameters  include  the  expected  interarrival  time  and 
expected  interservice  time.  The  probability  density 
functions  for  the  interarrival  time  and  interservice 
times  express  the  functional  relationship  between  the 
entities . 

The  first  three  sets  --  entities,  attributes, 
and  functional  relationships  --  comprise  the  static 
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description  of  the  system.  The  static  model  (M)  is 
then  tho  union  of  the  three  sets 


M  -  E  IT  A  V  F 


The  state  of  the  model  at  time  T  (s  )  is  the  value 
assigned  at  that  time  to  the  variables  associated  with 
all  entities  in  the  model  i.e., 

ST  •  sr(ej)  :  for  all  j,  ej  E  at  timo  t| 

The  state  of  entity  e^  may  also  be  written  sj  when 
the  time  t  is  obvious  from  the  context  or  usage.  The 
state  space  is  the  set  of  all  states  of  the  model 


S  -  {V  :  T  €  3  ] 

where  $  is  the  time  horizon  of  the  simulation. 

The  principle  type  of  simulation  discussed  in  this 
report  is  discrete  event  simulation  ,  i.e.  »  simulation 
of  a  system  which  has  discrete  units  such  as  customers 
or  machine  failures  flowing  through  a  sequence  of  stations 
or  locations.  The  basic  simulation  structure  is  illustrated 
in  F  igure  Z2. 

A  primary  task  in  simulation  is  the  method  of 
representing  processes,  which  are  activities  that  proceed 
over  time.  The  initiation,  alteration,  or  conclusion  of 
an  activity  is  called  an  event .  Since  events  are 
associated  with  the  system's  entities,  the  state  of  the 
system  changes  if  and  only  if  an  event  occurs.  The 
significance  of  all  this  is  that  processes  are  not  modeled 
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1 


SIMULATION  CONTROL  PROGRAM 


DETERMINATION  OF  NEXT 
EXECUTABLE  SUBPROGRAM 


SUBPROGRAM 
MODULE  2 


J 


SUBPROGRAM 
MODUUS  n 

1 


Figure  22.  Basic  Structure  in  Discrete 
Event  Simulations  (After  Reference  18) 


explicitly  but  rather  are  represented  by  the  modeling  of 
the  events  which  affect  the  status  of  the  process. 

In  building  a  simulation  model  the  analyst  must 
construct  an  event  list  which  contains  descriptions  of 
events  computed  to  occur  at  some  future  time  in  the 
simulation.  The  list  is  updated  as  other  events  are 
added  and  as  the  system  state  changes.  The  event  time 
prediction  routine  or  time  flow  mechanism  is  the  heart 
of  the  simulation  model  since  it  determines  the  system's 
dynamics . 

Time  Flow  Mechanisms 

Time  flaw  mechanisms,  timing  routines,  event 
scheduling  procedures,  or  simulation  executive  routines 
are  widely  used  terms  which  will  be  treated  as  synonymous. 
Regardless  of  the  name,  its  function  is  to  advance 
simulation  time  and  select  a  subprogram  for  execution 
that  performs  a  specified  activity.  Historically,  time 
flow  mechanisms  have  been  categorized  as  belonging  to 
one  of  two  general  methods:  fixed  time  increment  or 
next  event  increment.  We  will  now  compare  the  logic 
involved  in  each  formulation. 

Fixed  Time  Incrementing 

Consider  a  sequence  of  events  Ej_,  E2,  ...,  En 
as  they  actually  occurred  with  real  world  times 

rl*  r2*  •••»  rn* 
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(a) 


ACTUAL 

EVERTS 


(b) 

FIXED 

TIME 

INCREMENT 


(c) 


NEXT 

EVENT 

INCREMENT 


*1  *2  *3  ** 


tlae 


tlM 


tlae 


*0  *1  *2  *3 


Figure  23.  A  Cjm^ipon  of  Real  World  Events  and 
Fixed  Time  increment,  and  Next  Event  Increment 
Simulated  Time 


They  are  depicted  as  (a)  in  Figure  23.  In  fixed  time 

4, 

incrementing  the  analyst  defines  &n  incremental  time 
value  Lt  and  the  fiodel  progresses  through  simulated  time 
by  At  increments: 

for  t0  ■  o,  tj*At,  t2*2At,  tj«3At,  t^AAt,...,  tj^-KAt 
Line  (b)  of  Figure  23  illustrates  the  fixed  time 
increments.  Time  advances  uniformly  from  t0  to  tj, 
tj  to  t2«and  continues  advancing  until  t^ .which  is  the 
end  of  the  simulation.  Since  time  is  simulated  only  at 
the  time  points  t, ,  t^,  t£»...»  t^,  the  events  Ej»  E2,..., 
En  do  not  occur  at  their  "real  world"  times,  fivmnt  is 
shifted  or  translated  in  time  to  occur  at  t^,  events 
E2  and  E3  are  shifted  to  occur  tt  tv  and  the  shifting 
continues  until  event  En  occurs  at  t^.  Note  that 
events  E2  and  E3 .which  really  occurred  at  different  times 
appear  to  occur  simultaneously  at  the  end  of  the  time 
interval.  Thus,  all  events  Having  end  times  r^,  r^j, 
r.^2»‘*‘>  r;+n  9UC^  that 

tk-l<rt*  r  i+2'  *  ’  ’ '  rt+n^tk 

are  treated  as  if  they  occurred  at  tj,.  The  analyst  must 
be  careful  in  choosing  At  so  that  the  interrelationship  of 
event  occurrences  are  not  radically  altered  solely 
because  of  the  time  increment  size.  In  practice,  this 
is  not  a  difficult  problem  since  At  can  be  chosen  small 
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enough  such  that  at  most  only  one  event  is  likely  to 
occur.  Some  experienced  researchers  use  the  rule  of 
thumb  that  a  reasonable  time  interval  is  one  tenth  of 
the  shortest  expected  inter-event  time  [7,P.  162]. 

Next  Event  Incrementing 

In  next  event  incrementing ,  time  points  are  defined 
only  at  the  occurrence  of  an  event  or  events.  Line  (c)  of 
Figure  23  illustrates  the  ne*t  event  increments.  Simulated 
time  started  at  time  t0  and  advanced  immediately  to  t\ 

(the  time  of  event  E|)  then  to  t2  (the  time  of  event  E2) 
and  continued  to  tn  (the  time  of  event  fin).  Thus,  rj-t,, 
r2-t2»  •••  rn-tn*  T**e  si*ula*i°n  *°del  updates  the 
status  of  those  variables  associated  with  event  E^,  processes 
Ej ,  determines  the  time  and  subsequent  events  which 
occur  as  a  result  of  Ej ,  and  then  advances  simulated  time  to 
the  time  of  the  next  event  Ej+i.  Two  events  are  not 
processed  as  simultaneous  events  unless  they  bear 
identical  occurrence  times . 

A  Comparison  of  Fixed  Tine  and  Next  Event 

Figures  24  and  25  depict  the  general  structure  of 
next  event  and  fixed  time  increment  time  flow  mechanisms. 

In  next  event  incrementing  the  simulation  master  clock 
is  advanced  by  the  amount  necessary  to  t uuse  the  next 
most  imminent  event  to  take  place.  When  a  particular 
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TIME  AT 
T  -  T  +  AT 


Figure  25.  Flowchart  For  Next  Event 
Tine  Flow  Mechanism 
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event  has  been  executed,  clock  time  is  advanced  to  the 
time  at  which  the  next  significant  event  is  to  occur, 
whether  seconds  or  hours  away.  The  intervening  periods, 
when  no  changes  occur  in  the  system,  are  skipped  over. 

Thus,  the  passage  of  time  is  dependent  on  event 
occurrences.  In  fixed  time  incrementing,  the  master 
clock  is  updated  in  uniform  discrete  intervals  of 
time.  The  system  is  scanned  every  unit  of  clock  time 
to  determine  whether  any  events  are  due  to  occur  at  that 
particular  clock  time.  Consequently,  the  occurrence 
of  events  is  dependent  on  the  passage  of  time. 

The  casual  observer  might  foolishly  conclude  that 
a  next  event  formulation  is  always  preferred  over  the 
fixed  time  increment  method.  After  all  with  next  event  method 
one  does  not  need  to  choose  an  arbitrary  and  artificial 
time  increment  and  events  are  not  processed  simultaneously 
unless  they  have  identical  occurrence  times.  The  above 
argument  is  overly  simplistic  and  a  hast/  generalization. 

Let's  look  at  the  tradeoffs  entailed  in  the  use  of  each 
method.  With  fixed  time  increments,  some  information  is 
always  lost  since  translation  to  the  end  of  the  time 
interval  always  results  upon  the  occurrence  of  an  event. 

False  simultaneity  can  be  induced  due  to  translation  of 
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events  within  the  sane  interval.  Furthermore,  accuracy 
can  be  lost  if  at  is  too  large  and  execution  can  be 
slowed  if  at  is  too  small.  Lastly,  the  model  must  always 
search  for  the  next  event  even  during  intervals  when  no 
event  occurs.  With  the  next  event  method,  execution  time 
could  be  longer  than  fixed  time  incrementing  for 
interdependent  events  due  to  increased  computations.  In 
addition,  next  event  is  cumbersome  when  one  or  more 
events  change  in  a  continuous  fasnion.  Conway  et  al  l6J 
has  concluded  that  the  desirability  of  the  fixed  time 
increment  method  increases  with  the  increase  in  number 
of  entities  and  the  desirability  of  the  next  event 
method  increases  with  the  mean  length  of  event  times. 
Naylor  [  24]  suggests  using  fixed  tiiM  incrementing  when 
events  occur  in  a  regular  manner  and  next  event  when  they 
occur  unevenly  in  time.  For  some  models,  one  method  will 
be  move  t  ffective  than  the  other, and  one  must  choose 
depending  on  which  is  best  for  the  given  situation. 

Further  Classification  of  Time  Flow-Mechanisms 

Earlie*  in  this  section  we  stated  that  historically 
time  flow  mechanisms  were  categorized  as  fixed  time  or 
next  event  incrementing.  Kiviat  [18]  extended  the 
classification  to  three  methods:  event  scheduling, 
activities  scan,  and  process  interaction,  and  Nance  [23] 
proposed  a  concept  of  a  continuum  of  algorithms  for  time 
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flow  in  digital  simulation.  Since  typically  one  is  not 
able  to  conclusively  separate  tine  flow  mechanisms  into 
just  two  methods,  we  will  examine  the  concept  of  a 
continuum  in  more  depth. 

The  continuum  has  fixed  time  incrementing  at  one 
pole  and  next  event  at  the  other.  A  particular  algorithm 
would  lie  somewhere  along  the  continuum  depending  on  the 
degree  to  which  it  possesses  characteristics  of  each 
pole.  Considering  the  innumerable  discrete  event 
simulations  possible,  there  are  no  doubt  an  in',...  te 
number  of  time  flow  mechanisms.  For  any  specific 
simulation  application,  the  most  efficient  algorithm 
will  be  somewhere  along  the  continuum. 

Nance  [23]  examined  several  time  flow  mechanisms 
for  the  patrolling  repairman  problem  to  determine  their 
relative  effectiveness.  The  patrolling  repairman  problem 
is  a  classical  machine  interference  problem  from  queuing 
theory.  A  single  repairman  is  assigned  to  service  a 
group  of  N  semiautomatic  machines  which  fail  intermittently. 
We  shall  assume  that  the  failure  rates  of  the  machines 
are  identical,  the  went  of  failure  for  any  machine  is 
independent  of  the  state  of  any  other  machine,  and  the 
time  between  failures  of  a  single  machine  and  the  repair 
time  are  distributed  negative  exponentially  with  respective 
means  X  and  yu  .  The  machine  layout  is  a  rectangular 
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pattern  with  two  rows  of  the  same  number  of  machines. 

The  repairman  requires  time  T  to  walk  between  any  two 
adjacent  machines.  The  repairman  unidirectionally  patrols 
the  perimeter  of  the  machines  repairing  any  failed  machine. 
An  event  is  defined  as  the  failure  of  a  machine. 


STEr 


OrtMEHT 


1.  J  •  J+l 

2.  IF  (.T>N) ?  J*l 

3.  Clock  -  Clock  ♦  T 


Advance  t  e  next  facility  in  path 

Cheer,  for  correct  facility 
number  >c>  /;  is  not  exceeded. 

If  true,  set  i - 1 ;  otherwise,  go 
to  .3 

Increment  ,  ime  ny  Constant 
vain.: 


4.  IF  (CLOCK Jffj) 
next  event  at  j 


5  Go  to  J . 


Test  *  .auCiii.»e  nas  failed. 

Let  ii-  failure  time  of  machine 
j.  if  true  Uie  next  imminent 
event  (cents  at  j;  otherwise, 

go  to  J,  . 


The  algorithm  records  which  macniiv  i  cuing  inspected  (j ) , 
increments  CLOCK  by  Tf  and  compares  1  «itn  tne  failure 
time  of  machine  j  (fi).  If  the  mas  lev  clock  time  is 
greater  than  or  equal  to  f  j  ,  then  the  j  * machine  has 
failed  and  the  next  imminent  event  occurs  at  machine 
j.  J.f  CLOCK  is  less  than  f  j  ,  then  the  machine  is 
operational  and  j  is  incremented  1/  one  On  the  continuum 
the  time  flow  algorithm  can  be  depicted  follows: 


•— - - - K . . . — X . . a 

FIXED  TIME  <-:XT  EVENT 

INCREMENTING  OA  MCNTINC 

Execution  time  as  a  functi  on  ,  s .  ume tars  /\  and  M  , 

the  number  of  machines,  ana  uu..m.-r  :>r  repairs  was 
studied  for  each  algorithm.  As  lie;  »  ot  aiachirc*- 

increased  the  differences  among  tin?  a  1  K,<»  i  c ; .  in  s  became 
pronounced.  The  conclusions  were  i  ai  .  i  as  mox' 
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efficient  than  MICA  as ^  increased,  MICA  was  store  sensitive 
and  less  efficient  than  Cl  as  intermachine  transit  time 
(T)  increased,  and  regardless  of  the  number  of  repairs 
generated  (within  the  test  values  of  N,  T,  and  t  ci 
was  more  efficient  than  MICA.  The  study  indieated  that 
MICA  suffered  from  a  requirement  for  too  much  computation 
time  per  test  and  an  increasing  number  of  machines  tested 
due  to  an  increase  in  machines  or  intermachine  transient 
time. 

The  most  important  result  from  this  study  is  that 
the  procedure  by  which  the  passage  of  time  is  structured 
most  advantageously  depends  on  the  system  being  modeled. 

It  may  display  in  various  degrees  characteristics  of  the 
two  poles,  i.e.,  fixed  time  or  next  event. 
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SECTION  V 


OUTHJT  GENERATION  AND  STATISTICAL  ANALYSIS 

Frequently  the  areas  of  output  generation  and 
statistical  analysis  of  the  model  results  are  given  little 
emphasis  by  the  model  builder.  But  the  analyst  should 
keep  in  mind  that  decisions  about  structural 
characteristics  of  the  model  must  be  made  in  parallel  with 
decisions  about  how  the  model  will  be  used.  This  Implies 
that  decisions  about  hov  detailed  the  simulation  should 
be  are  affected  by  what  output  is  desired.  The  output,  in 
turn,  should  'satisfy  the  requirement?  of  the  decision 
mailer .  who,  therefore,  should  be  involved  in  output  design. 

Output  Requirements 

The  analyst  must  address  some  fundamental  questions 
before  devising  the  output  routines.  Consistent  with 
the  objectives  of  the  study,  he  must  decide: 

•  What  information  is  required, 

•  What  statistics  provide  this  information, 

•  What  analysis  techniques  provide  sufficient 
statistics  on  model  behavior,  and 

•  What  format  best  commun'  _ates  the  information 
derived  from  model  results  and  analysis. 

The  analyst  must  also  consider  whether  to  combine  the 

data  analysis  routines  in  the  simulation  model  program 

or  provide  information  in  some  form  to  be  used  by  a 
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separate  analysis  program.  In  addition,  the  output  itself  could 
be  punch  cards,  tape  or  disc  files,  or  printouts. 

Lastly,  diagnostic  information  about  the  model  should 
be  incorporated  *.  .th  output  information  derived  from 
the  model.  In  this  way  data,  say  from  the  random  number 
generator,  can  be  compared  to  theoretical  values  of  the 
distribution  assumed. 

If  the  simulation  is  deterministic,  there  are  no 
problems  of  statistical  inference.  The  output  of  the 
simulation  at  the  end  of  the  computer  run  is  the  required 
measure  of  system  performance.  Thus,  we  will  henceforth 
treat  probabilistic  models. 

Basic  Terminology 

Before  beginning  the  discussion  of  statistical 
analysis  of  the  model  results,  it  would  be  wise  to  define 
some  basic  terminology.  A  simulation  run  is  an 
uninterrupted  exercise  of  the  model  for  a  specified 
combination  of  controllable  variable-*  cx  parameters.  A 
replication  (for  stochastic  models)  of  a  run  is  an 
exercise  for  the  same  combination  but  with  different 
random  variations.  An  observation  from  the  model  is  a 
segment  of  a  run  sufficient  for  estimating  the  value  of 
some  measure  or  atatistic. 
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Preliminaries  to  Analysis 

Once  the  general  outline  of  the  output  requirement 
has  been  established,  certain  statistical  considerations 
aust  be  recognized.  A  simulation  model  is  concerned 
with  producing  a  stochastic  process,  i.e.,  a  family  of 
random  variables  |x(t),  tioj  indexed  by  a  time  parameter 
t.  The  analyst  is  interested  in  the  behavior  of  some 
variable  X  as  a  function  of  time.  This  time  dependent 
behavior  produces  results  which  are  not  independent  but, 
in  fact,  interdependent  to  some  degree.  For  example,  in 
a  queuing  model  for  any  single  simulation  run,  the  waiting 
times  of  successive  customers  will  be  autocorrelated 
since  there  is  a  greater  likelihood  that  the  (j+l)st 
customer  will  be  delayed  if  the  jth  customer  waits,  than 
if  the  jth  cas tome r  were  served  immediately.  To  determine 
a  statistically  valid  measure  of  system  performance,  the 
analyst  must  replicate  experiments  with  varying  sequences 
of  pseudorandom  numbers.  Since  the  behavior  of  the  time 
dependent  process  may  be  irregular,  e'.:i„nded  run  lengths 
may  be  required. 

After  the  model  is  constructed,  the  analyst  is  faced 
with  designing  the  simulation  experiment  --  the  set  of 
runs  for  the  model.  He  must  determine  the 

•  initial  state  definition,  i.e.,  statting 
conditions  for  the  model 

•  criteria  for  recognition  of  steady  state. 
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•  parameter  settings  to  expose  different  system 
responses 

•  length  of  each  run,  and 

•  sampling  procedures  including  sample  size  and 
aethod  of  selection. 

Initial  State  Definition 

Unlike  the  rail  world  systea,  the  siaulation  aodel 
is  not  in  continual  use.  The  analyst  is  faced  with  the 
problea  of  how  to  start  the  aodel  and  obtain  aeasureaents 
which  are  not  biased  by  the  aethod  of  starting  or  stepping. 

Simulation  is  frequently  used  to  study  the  perforaance 
of  a  systea  which  operates  under  steady-state  conditions. 

The  analyst  attempts  to  determine  fche  limiting  distribution 
of  the  state  of  the  systea  or  aore  precisely  empirical 
estimates  of  the  distribution's  moments  such  at  mean  and 
variance. 

A  system  whose  behavior  does  not  satisfy  steady-state 
conditions  is  considered  to  be  in  a  transient  stars. 

When  starting  conditions  are  not  near  steady-state,  a 
transient  period  exists  until  the  state  of  the  model 
approaches  the  steady-state  condition.  Implied  in  the 
above  is  the  assumption  that  a  steady-«tate  exists.  For 
some  real  world  systems,  a  transient  phenomenon  may  actually 
occur  prior  to  the  steady  state  and  the  transient  period 
may  be  pertinent  to  the  study.  For  other  systems,  no 
steady-state  exists  and  it  is  the  transient  phase  itself 
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which  is  to  be  studied.  In  a  strategic  air  defense  model, 
it  is  the  transient  state,  where  aircraft  initially  attempt 
to  penetrate  the  air  defense  system  and  destroy  their 
assigned  targets,  that  is  of  prime  importance,  hither 
no  steady-state  exists  because  of  the  system's  dynamic 
and  reactive  capabilities,  or  it  is  not  reached  because  the 
air  battle  does  not  last  sufficient  time  or  because  there  is 
not  an  infinite  supply  of  aircraft  or  other  resources. 

Let  us  consider  some  strategies  for  setting  initial 
conditions.  An  appealing  method  is,  to  use  queuing 
terminology,  empty  and  idle,  i.e.,  the  system  has  no 
activity.  The  simulation  is  then  run  until  the  transient 
effects  are  apparently  insignificant.  However,  Conway 
[ 5 J  suggested  that  empty  and  idle  is  a  poor  state  to 
consider  sirice  almost  anything  else  is  better.  In  fact, 
unless  empty  and  idle  is  a  typical  state  of  the  system  we 
snould  choose  another  condition.  Ideally,  we  would  like 
to  select  starting  conditions  corresponding  to  the  steady- 
state  condition.  But  an  analyst  who  has  already  determined 
what  the  equilibrium  state  of  the  system  is  either  does 
not  need  simulation  or  may  bias  the  results  to  his 
preconceived  conclusions  if  data  is  collected  before  the 
transient  effects  are  removed.  Thus,  we  can  only  conclude 
that  the  analyst  probably  knows  something  about  the 
equilibrium  behavior  of  the  system  and  should  use  this 
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information  to  select  a  starting  value.  Realistic  initial 
conditions  will  reduce  the  cost  to  reach  «  *.’y  state 

and  will  affect  the  precision  of  statistical  estimates. 
Simulation  models  are  frequently  used  to  compare 
twc  alternative  system  configurations  cr  disciplines.  To 
avoid  biasing  the  results  and  retain  a  basis  for  comparison, 
the  initial  state  definition  for  both  configurations  should 
be  the  same,  such  as  a  compromise  or  average  of  the  expected 
steady-state  levels  for  each. 

Two  methods  are  commonly  used  to  remove  transient 
effects  in  simulations.  The  first  is  to  use  long  simulation 
runs  so  the  data  from  the  transient  period  is  insignificant 
compared  to  the  amount  of  steady-state  data.  Although 
this  method  is  simple  to  employ,  it  can  be  costly  in  terms 
of  total  running  time  and  computer  charges. 

The  second  method  is  to  run  the  simulation  until 
steady-state  conditions  are  achieved  and  then  truncate 
atl  sampled  data  up  to  that  point.  The  conditions  at  the 
end  of  the  transient  state  become,  in  effect,  the  initial 
state  for  the  data  collection  phase  of  the  simulation. 
Theoretically,  one  could  start  the  model  anywhere  and 
merely  run  the  simulation  for  a  sufficient  period  of  time 
until  transient  effetts  are  removed  and  truncate  data. 
However,  Fishman  {  10]  has  shown  that  if  mean  square 
error  (MS£)  is  used  as  a  criterion,  where 
MSE  •  variance  ♦  (bias)2 
10, 


then  truncation  decreases  the  sample  size  needed  to 
reduce  bias  but  can  inflate  variance  so  that  mean  square 
error  it  not  reduced.  We  cannot  say  categorically 
that  the  truncation  strategy  reduces  mean  square  error. 

Later  in  this  section  we  will  discuss  variance  reduction 
techniques  which,  when  used  together  with  truncation, 
provide  an  acceptable  treatment  of  transient  effects. 

Recogniti on  of  Steady-State 

Whether  we  are  interested  in  transient  behavior 
or  steady-state  conditions,  we  require  a  criterion  for 
recognition  o i  steady-state.  Equilibrium  is  a  limiting 
condition  which  * ay  be  approached  but  never  actually 
attained.  There  is  no  single  point  in  the  simulation 
beyond  which  the  system  is  at  steady-state.  We  will 
assume  that  the  difference  between  the  temporal  and  limiting 
distribution  decreases  with  time  and  that  bey-ond  a 
given  point  one  is  willing  to  neglect  the  error  made  by 
considering  the  system  to  be  at  equilibrium.  Thus,  we 
assume  that  Lj.  (t),  k*l,2,...,a  is  the  limiting 
distribution  for  an  n  state  process  and  Sk  (t)  is  the 
simulated  temporal  distribution.  Then 
|Lk(t)-Sk(t',j«Ek(t) 

where 

lira  Ek(t)  ■  0  foi  k  *  1,2,. . . ,  n 
t-»oo 

Although  there  are  no  fixed  rules  for  determining 
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when  transient  behavior  ends,  several  methods  are  in 
current  use.  A  simple  method  is  to  examine  a  sequence 
of  observation*  and  assume  steady-state  begins  at  the 
first  point  which  is  neither  the  maximum  nai*  the  minimum 
of  the  preceding  sequence.  A  second  method  is  to 
select  the  print  which  is  the  first  value  to  repeat  itself 
after  the  start.  The  third,  and  probably  thr,  ircst 
frequently  used  method,  is  to  select  the  time  urh«r;  three 
consecutive  points  differ  hy  no  more  than  E,  say  E>».01r 
Each  technique  will  differ  in  the  point  selected  for 
s?*ady-stste  and  all  may  be  good  values.  The  analyst 
must  choose  fhich  method  is  most  appropriate  for  his 
simulation. 

Sampling  Procedures 

It  is  not  difficult  to  .estimate  the  means  of 
attributes  associated  with  permanent  and  temporary  entities. 
But  is  is  much  more  difficult  to  state  a  level  of 
confidence  that  the  mean  is  the  true  mean.  This  requires 
the  variance  of  the  data.  In  the  following  discussion 
we  will  first  examine  the  measurement  of  simulation 
results  for  permanent  and  then  for  temporary  entities. 

Permanent  Entities 

Permanent  entities  exist  throughout  the  simulation 
and  their  attributes  take  on  different  values  at  each 
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As  we  previously  mentioned  problems  arise  with  the  variance 
estimate  due  to  autocorrelation  introduced  in  calculating 
values  associated  with  permanent  entities. 


Independent  Replicates 

Ono  obvious  method  of  eliminating  problems  of 

autocorrelation  is  to  replicate  each  run  several  times. 

The  mean  from  each  run  is  treated  as  one  observation  in 

the  total  sample  used  to  calculate  the  appropriate 

statistic.  Of  course,  on  each  r'*n  the  transient  data 

prior  to  steady-state  is  lost  leading  to  higher 

computer  run  costs.  For  sample  size  n,  the  estimate  of 

n 

the  mean  is  X  *  I*  xi(/n 


and  the  estimate  of  the  variance  of  the  sample  mean  is 

Var(X)-l  £  (X.-X)2 

n  j.-i  1 

or  equivalently, 

Var(X)-  i  (LxJ  -  n  X2). 

The  reader  should  note  that  the  latter  formulation  ef  the 
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sample  variance  estimate  is  easier  to  implement  in  a 

2 

computer  program  since  only  the  sums  of  and  need 
be  accumulated  and  no  individual  X^  values  stored. 


Stratified  Sampling 

To  avoid  the  loss  of  data  over  the  interval  (Tr,Tt), 
consider  replications  by  dividing  a  single  simulation  run 
into  n  parts  over  the  interval  (Tr,Tt).  This  prevents  the 
waste  of  (n-1)  periods  of  (0,Tr)  but  succeeding  sample 
values  are  not  independent. 

Let's  define  a  measurement  process  in  terms  of 
measurement  periods  where  an  estimate  E[X(t)]  is  determined 
and  nonmeasurement  periods  in  which  no  statistical  data  is 
initially  gathered  [5],  Suppose  we  have  somehow 
determined  an  interval  of  length  so  the  measurement  periods 
are  independent.  Then  there  are  n/2  measurement  periods 
for  which  attribute  values  X2,  X^,  ...,  Xn  can  be 
determined.  We  will  assume  that  these  values  are  n/2 
independent  identically  distributed  random  variables  with 
mean  and  variance  0‘x.  We  have  then  for  the  measurement 
periods 


n/2  X2k 

k-1 

Var  (X)  -  iLX-  2  £* 
n/2  n 

and  for  the  nonmeasurement  periods 
n/2 

Y  "  £  x  tj, 

i-i  2k+i  -  r 
n/2  10S 


Var  (y)  -  2g^ 


n 

We  can  pool  the  means  so  that 

2  “  "  %  *i 

TP  l*z  Tr¬ 
an  d 


E[  1  ]  »  E[  X+Y]  -  E[X1  ♦  E  [Y] 
2  2 


•  !*-  /* 


The  lack  of  independence  prevents  a  directly  additive 
form  for  the  pooled  variance. 

Var  (X+Y)  -  1/4  [Var  X  ♦  Var t  ♦  Cov  (X,Y)] 

T~ 

™  1/4  [Var  X  ♦  Y]  ♦  l/4  Cov  (X,Y) 

-  1/#  [Var  X  *  Var  ?L  ♦  i  XYj  •?) 

■  1/4  [Var  5T  ♦  Var  YJ  ♦  aCn-ljc 

n*  ' 

where  C  is  the  covariance  term  (X-X)(Y-Y) 

between  adjacent  periods.  The  covariance  between 

nonadjacent  periods  is  assumed  to  be  zero.  In  fact  the 

length  of  the  period  is  selected  so  l i  this  is  true 

The  importance  of  this  technique  lies  in  the  fact 

that  the  pooled  variance  is  smaller  than  the  variance 

from  either  se.  of  periods.  Since 

Var  fl  )-l/4(Ver  X  ♦  Var  Y)  ♦  2(n-l)c 

n'2 

-1/4  ♦  2T2)  ♦  2  (n^l ) c 

n  n  n* 

"£?.  +  2  (n-l)c 
n 
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We  have  Var  ( £ )  ■  n(T^+2 (n~l)c 

i>2 


It  can  be  shown  thnt  C40^ 

“7“ 

so  Var  (I  )  ■  n<r2+2fn»l)C  ^q2  1*  (n- 1) 

IT2-  n~  n 

<  2C2  *  Var  (X) 

n 

This  variance  reduction  technique  is  called  stratified 
sampling.  It  is  also  sometimes  known  as  blocking  since  the 
data  in  essentially  divided  into  old  and  even  blocks  from 
which  sample  means  are  calculated. 

Estimating  Autocorrelation 

In  the  above  version  of  stratified  sampling  the 
length  of  each  measurement  period  wis  selected  so  that 
consecutive  measurement  periods  were  independent,  i.e., 
the  correlation  was  zero.  In  some  cases  nonadjacent 
periods  may  have  a  positive  covariance  that  is 
significant.  Then  the  autocorrelation  function  can  be 
estimated  and  included  in  estimates  of  the  variance. 

The  autocorrelation  function  between  an  observation 
at  simulated  time  t  and  an  observation  at  t+s  is 
R<5)“ 

t»i 

for  s-0,1,2, . . . ,  n-1 
where  Xt"  observation  at  time  t 
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X  *  sample  mean 
q  •  total  time  of  simulation 
s  -  lag 

The  variance  is  then 

i  f  n- 1 

Var(Z)  «  nfc  Ig2^  T\  (n-s)  R(s) 

L  s-i 

Antithetic  Variates 

Another  variance  reducing  method  is  Antithetic 

variates.  The  aim  is  to  introduce  negative  correlation 

between  two  separate  replications  of  the  simulation,  so 

the  variance  of  the  combined  averages  is  less  than  if 

the  replications  were  independent.  The  variance  of  the 

mean  of  two  replications  is 

var  (k)  -  <T4(1*P) 

i  ^ 

where  a  ■  population  variance 

p  ■  correlation  between  pairs  of  observations. 

If  the  observations  are  independent,  P-0.  But  if  there  is 
negative  correlation  (P<0) ,  then  the  variance  of  the 
sum  of  observations  will  be  reduced. 

The  procedure  most  commonly  followed  to  generate 
negatively  correlated  variates  is  to  use  uniform  pseudorandom 
numbers  [Vj]  for  probabilistic  events  on  one  run  and  to 
use  [1-Vj]  for  the  equivalent  event  on  the  second  run. 

It  is  important  to  maintain  event  equivalence  between  runs 
since  otherwise  the  processes  will  be  out  of  phase  and 
proper  negative  correlation  will  not  result. 
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Sample  Size  Determination 

The  determination  of  an  appropriate  sa.^iple  size  for 
simulation  is  no  different  from  sample  size  determination 
in  ordinary  statistical  problems.  The  added  complication 
is  that  two  sizes  must  be  determined:  the  run  length 
and  number  of  replications. 

If  stratified  sampling  is  used,  the  length  of  each 
block  or  measurement  period  can  be  determined  by 
examining  the  autocorrelation  function.  The  number  of 
,  replications  or  number  of  blocks  is  dependent  on  what 
the  analyst  believes  is  a  valid  sample  size,  Since 
computer  time  is  not  free  and  execution  times  for  complex 
simulations  can  be  great,  sample  sizes  tend  to  be  modest. 
The  clear  preference  is  for  large  sample  sizes  since 
large  sample  statistics  (n£30)  allow  the  Central  Limit 
Theorem  to  be  invoked  and  normal  approximations  made  for 
hypothesis  testing. 

Fishman's  work  with  the  spectral  density  function 
has  led  to  procedures  for  determining  the  sample  size  for 
equivalent  independent  observations  [9].  The  crucial 
factor  in  deciding  to  use  sophisticated  techniques  is 
whether  the  benefits  are  sufficient  to  warrant  the  extra 
computations  required. 

Temporary  hntities 

All  of  the  above  procedures  were  discussed  in  the 
context  of  measuring  attributes  associated  with  permanent 
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entities.  The  problems  relating  to  temporary  entities 
are  somewhat  different. 

Temporary  entities  are  active  for  only  a  portion 
of  the  total  simulation  and  usually  exist  in  such  numerous 
quantities  as  to  prevent  maintaining  records  on  each 
individual  entity.  Typically,  only  the  final  value  of 
an  attribute  is  measured  for  a  temporary  entity  before 
it  becomes  inactive.  Thus  each  entity  contributes  a 
single  value  and  the  sample  size  is  the  number  of  entities. 
Temporary  entities  existing  at  the  same  time  are  subject 
to  the  same  system  conditions  so  that  the  attributes  tend 
to  be  correlated.  Another  difficulty  involves  which 
entities  to  include  since  temporary  entities  often  are 
not  cheated  and  destroyed  in  the  same  order.  For  example, 
in  a  simulation  of  customers  shopping  in  a  commissary  the 
people  usually  do  not  leave  in  the  same  order  in  which 
they  entered. 

Conway  {5  J  suggested  three  strategies  for  collecting 
a  sample  of  temporary  entities.  The  first  was  to 
designate  an  interval  of  time  and  include  in  the  sample 
the  final  value  of  the  relevant  attribute  for  those 
entities  whose  active  history  terminated  during  the 
interval.  In  comparing  say  two  alternatives,  the  sample 
size  N  may  vary  but  the  time  interval  is  fixed.  The 
secodd  strategy  was  to  specify  a  beginning  point  and  a 
sample  size  and  include  the  required  number  for  consecutive 
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terainations  after  this  point.  Here  the  staple  size  is 
fixed  but  'Jie  duration  of  the  run  aay  vary.  Furtheraore, 
in  both  of  the  above  strategies  entities  which  were 
created  during  the  transient  period  aay  be  included  in 
the  staple.  The  third  stra >egy  is  to  specify  a 
beginning  point  and  staple  s$ze  and  include  all  attribute 
values  for  teaporary  entities  creates  consecutively  after 
this  point.  Although  no  transient  effects  are  included, 
a  problea  of  teraination  arises  since  the  run  must 
continue  until  all  teaporary  entities  in  the  saaple  are 
no  longer  active. 


Confidence  Intervals 


In  the  statistical  development  thus  far,  we  have 
been  concerned  with  reducing  the  variance  of  the  sample 
mean  since  it  expresses  in  some  fashion  the  reliability 
or  variability  of  our  estimate  of  the  mean.  To  give  better 
definition  to  our  statistical  estimates,  we  need  to 
examine  the  idea  of  confidence  intervals. 

Confidence  intervals  are  ‘based  on  statistical  properties 
of  the  sample  statistic  and  express  the  probability  1-oc 
that  the  random  variable  will  take  on  values  within  the 
constructed  interval.  The  preassigned  probability  l-o( 
is  called  the  decree  of  confidence. 

Below  we  will  aonsider  confidence  intervals  for 
some  standard  distributions.  Let  X1#xx,  ...,  Xn  be  a 
random  sample  with  the  sample  mean  denoted  by  X. 

Mean  of  a  Normal  Distribution  with  Known  Variance 

Suppose  we  are  sampling  from  a  noruai  population 

2 

with  mean  and  known  variance  CT .  The  random  variable 

X  ■/» 

r*  <r4nr 

is  distributed  normally  with  mean  0  and  variance  1.  Let 
^<*/2  be  such  *****  *he  integral  of  the  standard  normal 
d^.aSity  from  ^^to  infinity  equals  */z..  Then  the  random 
variable  R  will  take  on  a  value  between -Z^and  2# ^ 

with  probability  l-c< .  We  then  have  that 
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— 


X-/< 

*747T 


L  ^  <&/% 


cr 


X-  JL  <  /*  <■  *  +  2*4  X. 

>fn“  'f  n 


where  X  is  the  settle  mean,  For  given  c(,  the  value 
can  be  obtained  from  a  standard  normal  cumulative 
probability  table.  The  last  expression  above  is  the 
confidence  interval  for  the  mean  when  the  variance  is 
known; 

Mean  of  a  Normal  Distribution  with  Unknown  Variance 

In  some  instances  we  may  know  that  we  are  sampling 
from  a  normal  population  but  do  not  explicitly  know  the 
variance.  It  can  be  shown  [13,  p.222]  that  for  random 
samples  of  size  n  from  normal  populations  the  random 
variable 


has  a  t  distribution  with  n-1  degrees  of  freedom.  The 

random  variable  will  take  on  a  value  between  -t^  ,  n-1 

and  t  till  ,  n-1  with  probability  l-°<  .  Let  the  (unbiased) 

sample  variance  be  defined  as 
1  ^ 

s2  "  n71  Z  (V*)2 

\*  I  — 

Then  for  a  sample  with  mean  X  and  standard  deviation  S. 

z«lz,  n-1  4  X 2  j±  4  t«(/i  ,n-l 
or  S/AJTT 

The  last  expression  is  the  l-<  confidence  interval 
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for  J+n  when  fl-' 


is  unknown  for  normal  populations. 


Large  Sample  Confidence  Intervals  For  the  Maan 

Often  we  are  able  to  deal  with  samples  sufficiently 
large  (n>30)  to  justify  the  use  of  the  Central  Limit 
Theorem,  If  the  variance  of  the  population  is  known ,  then 
the  1*  ot  confidence  interval  for  f*  is 

X-  Z.C/2..JT  <  U  <  X  +  & 

%Tn 

If  the  variance  is  unknown ,  then  an  approximate  large 
sample  confidence  interval  for  jl\  is 


X’  A  <  U,  <  J  ♦  1 

TO  ' 


TO 


where  S  is  the  sample  standard  deviation. 


Variance  of  a  Normal  Distribution 

Suppose  we  have  a  random  sample  of  size  n  from 
a  normal  population  and  we  desire  a  1-tf  confidence  interval 
for  (T2.  It  can  be  shown  [  13,  p . 214 ]  that  the  random 
variable 

V'  -  \2 

L  (xt  -  x ) 
ui  <r*- 


has  a  Chi  Square  distribution  with  n-1  degrees  of  freedom. 

'V1 

The  random  variable  will  take  on  a  value  between  ri-ttfzin-i 
t 


and 


^  d!% ,  n* 


v/ith  probability  1-  .  Then  we  have 


^  cr1  L 


*  0-*  <  2  <*-*>■ 

■V1-  «Y  ** 

f  i-d/a.n-i  A-  oi/2>n>i 

The  last  expression  is  the  l-o(  confidence  interval  for 
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the  variance  of  a  normal  population. 


>othesis  Testini 


A  simulation  study  is  often  employed  to  compare  two 
or  more  strategies  or  alternatives.  One  means  of  comparison 
is  to  use  statistical  tests  of  hypotheses  where  somm 
assumption  is  made  about  the  probability  distribution  of 
the  population.  It  is  the  nature  of  statistical  testing 
that  rejection  of  a  hypothesis  can  be  more  positively 
stated  than  acceptance.  The  hypothesis  formulated  for 
testing  is  called  the  null  hypothesis  and  is  denoted  H0. 

The  hypothesis  which  contradicts  the  null  hypothesis  is 
called  the  alternative  hypothesis  and  is  denoted  H],  or  Ha. 

Two  types  of  errors  may  be  committed  while  testing 
an  hypothesis.  If  the  null  hypothesis  is  trui  and  the 
analyst  rejects  it,  he  is  6aid  to  have  committed  a  Type  I 
Error.  If  the  null  hypothesis  is  not  rejected  when  in 
fact  it  is  false,  he  is  said  to  have  committed  a  Type  II 
Error. 

In  order  to  test  an  hypothesis  the  analyst  must 
construct  regions  of  acceptance  and  rejection  for  sample 
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values.  The  set  of  values  for  which  the  null  hypothesis 
is  rejected  is  called  the  critical  region.  The  size  of 
critical  region  is  the  probability  of  committing  a  Type  1 
Error  and  is  usually  denoted  as  0(  .  The  size  of  the 
Type  II  Error  is  the  probability  that  a  sample  value  will 
be  outside  the  critical  region  (therefore  in  the  region 
of  acceptance)  when  the  null  hypothesis  is  false.  As 
long  as  the  sample  size  is  fixed,  an  inverse  relationship 
exists  between  Type  I  and  Type  II  errors:  if  the 
probability  of  one  error  is  reduced,  the  probability  of  the 
other  is  increased.  The  only  way  to  reduce  the  probabilities 
of  both  errors  is  to  increase  the  sample  size. 

Tests  are  generally  classified  as  one-sided  or  two- 
sided,  In  a  one-sided  test,  one  tail  of  the  distribution 
of  the  test  statistic  is  the  critical  region.  In  a  two- 
sided  test,  both  tails  are  used.  Selection  of  one  or 
two-sided  tests  depends  on  the  nature  of  the  alternative 
hypothesis.  Usually  two-sided  alternatives  lead  to  Owe- 
tailed  tests  and  one-sided  alternatives  lead  to  one-tail 
tests.  For  instance,  for  the  null  hypothesis  Ho;  s  /*>o 
and  the  alternative  :  jA4J*0t  we  would  have  a  two-sided 
test.  If  the  alternative  were  H^:  the  test  would 

have  been  one-sided. 

Often  in  hypothesis  testing,  the  alternative  to  rejecting 
the  null  hypothesis  is  veserving  judg  ment.  Since  the  null 
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hypothesis  is  not  accepted,  we  cannot  commit  a  Type  II 
error.  Such  tests  are  called  tests  of  significance 
and  the  probability  of  a  Type  I  error  (0< )  is  called  the 
level  of  significance.  The  most  frequently  used  values 
of  <X.  are  .  05  and  .  01. 


Sampling  from  Normal  Populations 

Let  X},  X2,  ...»  Xn  be  a  sample  of  size  n  from  a 
normal  population.  Let  X  be  the  sample  mean  and  o(  be  the 
size  of  the  critical  region. 

The  Mean  of  Normal  Population  with  Known  Variance 
Ho  JA.  — 

Hi  :  >  /*o  ;or  ^  fio  f  or  jJL  #  fJc 

The  test  statistic  is 

l  -  x 
<t/J7T 

and  the  critical  regions  are 


i  * 

** 

for 

^  >  ju, 

for 

yu  4  }xt 

1*1 1 

for 

/*  1  po 

The  Mean  of  a  Normal  Population  with  Unknown  Variance 
Ho  -•  o  /t„ 

H1  «  ^>^0;  or  ^UyH.  1  or  yU*yU0 

The  test  statistic  is 
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t.  X 

S/JT 

2 

where  S‘  is  the  sample  variance. 


:ritical  regions  are 

t  >  t 

”  n-1  for 

yu  >  /*. 

1  £  -*«  ,n-l  for 

/*  <  /» 

Mi.  tor 

/»  *  A 

where  t  .  .is  such  that  the  integral  of  the  t  distribution 
ct  ,n-l 

of  n-1  degrees  of  freedom  from  negative  infinity  to 
equals  l-o(  . 

Difference  Between  Means  with  Known  Variances 

Let  ,  X2 ,  . . . ,  Xn  and  Y ^ ,  Y2  , . . . »  Ya  be  random 
sap’iles  of  size  n  and  m  from  two  normal  populations  with 
variances  0T2,  and  CT^  respectively. 

Ho:  /^i  —  ^-*2.  -  ci 

Hl:  -  jix  >  <f:  or  pr  n%‘&  ,or^px4£ 


The  test  statistic  is 

2  ■  *  -  *  -  £ 
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The  critical  regions 

are 

for 

2s-2< 

for 

|i|  2. 

for 

Difference  Between  Means  for 


The  hypotheses  are  the 
The  test  statistic  is 


/*.  -^1  ><f 

/*\  ~  <f 

/A, 

Unknown  Variances 

axe  as  in  the  previous  case. 


where  S 


t-  X  -  Y  -  <f _ 

]  [  ( i/r\H  C«/m)  J  [  (n-l ) 5,z  4  (*-» )  S*  ]  /  (n+vn-Z ) 

2  2 

,  and  S*  axe  the  sample  variances. 

2 


The  critisal  regions  are 


for 

^  n-rvl 

for 

|  >  jkw/i ,  n-w>i 
Variance  for  Unknown  mean 

for 

/<!  _/<l  ^<f 

Ho:  cr1^  (T1 

Hi:  cr1  ■>  g,1  » °r 

The  test  statistic  is 

or1*  cr*  » or 

<r V  (r.1 

X2-  S  (Xi-X)2 


o;2 

The  critical  regions  are 


X2 

> 

t' 

for 

X2 

4 

for 

X2 

< 

(  H</z)  ,n-l 
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O'2  >  < 

g1  4  cr* 


or 


for 


X2 


«(/l,  n-l 


(T%  Q% 

I 


Where  2  is  from  a  Chi  Square  distribution  with  N-l  degrees 
of  freedom. 


Large  Sample  Tests 

When  sample  sizes  are  large  enough  (n^30)  for  the 
Central  Limit  Theorem  to  be  applied,  the  following  tests 
may  be  used. 


Mean  When  Variance  Is  Known 


H, 


h  =  A 


Hi:  jA  >  ^Iq  \ot  c  jko  ,ot 
The  test  statistic  is 

i  - 


<r/JTT 

the  critical  regions  are 


W  i  po 


2* 

for 

fA  >  ^Ac 

i  i  -2* 

for 

ft  ^  ft « 

1^1  ?  Zm/i 

for 

h 

where  ^4/2,  is  from  a  normal 

table. 

Mean  When  Variance  Is 

Unknown 

X 

0 

*> 

11 

0 

Hl:  H  *  /*o  ;or 

1 A  <^c;or 

h  *  /*° 

121 

The  test  statistic  is 

i  --  *-/*. 

3  /fr\ 


where  is  the  sample  variance. 


The 

critical 

regions  are 

1 

z 

for 

h 

>  /*• 

2 

£ 

for 

h 

CyM  0 

1*1 

> 

Icllz 

for 

h 

^  /^0 

Differences 

in 

Mean  When 

Variances 

are  Unknown 

Ho: 

/«■- 

h 

1  -  <T 

Hi: 

A*> 

h 

a  >  <T  ; 

or  K'fr 

^<£or 

{ 

<f 

The  test  statistic  is 

-  y,  -  y,  -S' 


2  ' 


y  2 

where  S“,  and  S2  are  the  sample  variances.  The  critical 
regions  are 


i  $r 

M  - 
ru 


ii 

r\% 


2  >  2* 

for 

/*■  -  2  2  <f 

2  i  -4* 

for 

yU,  -  c  <T 

IZl  2  2«/t 

for 

/>!  "  pi 
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Nonparaaetric  Methods 

When  a  normal  assumption  cannot  be  aet  or  when  large 
sample  testing  is  not  feasible ,  the  analyse  may  adopt 
nonparametric  methods.  One  such  nonparaaetric  test,  the 
sign  test,  was  introduced  in  the  discussion  of  pseudorandom 
numbers. 

Sign  Test 

Suppose  we  have  a  symmetrical  distribution  so  that  the 
probability  that  a  value  is  greater  than  the  mean  equals 
the  probability  that  it’s  less  than  the  mean.  We  will 
take  as  our  hypothesis 
H0:  /A  =  ^4, 

V  >  yu  , ;  or  f4</v° r  h  *  f*' 

Then,  we  will  test  the  null  hypothesis  ,  by  replacing 

each  sample  value  exceeding  ^ ,  by  a  plus  sisn(*)  and 
each  value  less  than yU,  ,  by  a  minus  (-)  and  testing  a 
new  null  hypothesis  that  the  plus  and  minus  signs  represent 
a  random  sample  from  a  binomial  population 

K0:  p  -  ,5 

Hj:  p  >  .5  (or  p  < ,5  or  p  /  .5) 

The  test  statistic  is  S,  the  number  of  plus  signs. 

The  critical  regions  are 

S  £  kj  for  p  >  ,  5 

s  -  ^d  for  p  <  .5 


123 


for  p  i  .5 


s*  kd/2 
or  s  £  K 


d/2 


where  kdis  the  smallest  integer  such  that 

b(y;n,.5)  ^  d 

y*k 

and 

b(y;n,.5)  is  the  binomial  probability  of  y  successes 
in  n  trials. 

Similarly,  k^is  the  lirgest  integer  for  which 
k*  d 

y*  b  (y  ;n, .  5)  £  d 

y-0 

the  values  k^and  k,  can  be  quickly  determined  from  a 
table  of  binomial  probabilities  for  the  given  sample  siae 
n  by  respectively  moving  down  the  table  of  Y  values  from 
y«0  until  the  cumulative  probability  is  d  *nd  by  moving 
up  from  y«n  until  the  cumulative  probability  is  d  . 

One  should  note  that  if  symmetry  is  not  a  valid  ?.s sump  cion, 
the  population  median  can  be  used  in  lieu  of  the  mem. 

Also,  if  n>100  the  normal  approximation  to  the  binomial 
distribution  can  be  used  and  the  test  statistic  becomer 

r . 

aI  np  (i-b) 

which  for  p*.5,  1  ' 
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I  -  2Vn 

JT 

The  critical  regions  are 


2  2 

He* 

for 

p  >.5 

l  i 

for 

P  *  .5 

121  2 

for 

p  4 .5 

Tables  of  binomial  probabilities  from  n*2  to  n*49  can  be- 
found  in  reference  34  and  from  n-SO  to  n-100  in  reference 
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Wilcoxon  Test 

The  Wilcoxon  test  also  known  as  Mann- Whitney  or  the 
U  test)  is  used  to  test  the  null  hypothesis  that  two  samples 
come  from  identical  populations.  The  two  sets  of  sample 
values  of  sizes  n^  and  ri£  are  taken  as  if  they  were  one 
sample  of  size  n^*^  and  the  vilues  are  arranged  in  order 
of  increasing  magnitude.  The  values  are  assigned  ranks 
1,2,3,...,  n1*n2*  We  observe  which  ranks  are  occupied  by 
each  sample.  If  there  are  ties,  each  is  assigned  the  rank 
which  is  the  mean  of  the  ranks  they  occupy,  i.e.,  if  ranks 
sixteen  and  seventeen  are  the  same,  each  is  assigned  the 
value  16,5, 
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The  rank  inforaation  used  to  determine  the  value  of  U, 

N,N,.  N1  <N1*U  -R, 

U  -  1  - 2 -  1 

where  nj,  n2  are  the  sample  sizes,  Rj  is  the  sum  of  the 
ranks  assigned  to  one  of  the  samples,  say  the  first. 

For  sample  sizes  nj  and  n£  of  less  than  eight,  special 
tables  of  the  U  distribution  are  required  as  in  reference 
26  •  For  nj,  and  n£  values  both  greater  than  eight,  a  normal 
approximation  to  U  can  be  used.  Then, 

H„:  /*.  */»J- 

Hi:  yu,  >  {*  i  ;  or  f*,  *  h,* 

The  test  statistic  is 

2  r  V  -  E(V) 

''I  V«r  (V) 

where  E  (U)  « (n^  n2  )/2 

2 

Var(U)  ■  n1n2(ni+n2+l) 

— n - 

The  critical  regions  are 


-  2* 

for 

i 

> 

2«( 

for 

c  M  i 

1*1 

> 

for 
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Verification  anc3,  Validation 

In  the  discussion  of  model  building  in  Section  II, 
we  concluded  that  a  simulation  aodel  is  an  approxiaate 
representation  of  reality  and,  in  fact,  an  abstraction  of 
the  real  world.  The  analyst  still  aust  deteraine  how 
accurately  the  aodel  portrays  reality.  He  hears  the 
recurring  questions  "Is  the  aodel  valid?"  and  "Is  your 
aodel  telling  ae  the  truth?".  Now  let's  exaaine  the  veri¬ 
fication  and  validation  process  which  has  been  called  "the 
most  elusive  of  all  the  unresolved  problems  associated 
with  computer  siaulation"  [24,  p.310]. 

First,  we  will  define  some  basic  terminology. 
Verification  is  ensuring  that  a  simulation  model  behaves 
as  the  analyst  intends.  Validation  is  testing  the  agree¬ 
ment  between  the  behavior  of  the  simulation  model  and  the 
real  world  system  [2], 

Verification  is  not  a  predetermined  number  of  steps 
but  a  continuing  effort  throughout  the  model  building 
process  to  ensure  that  the  model  behaves  as  intended.  In 
Figure  9  verification  of  the  mathematical  model  and  the 
computer  program  are  discretely  represented  in  the  flow¬ 
chart.  For  stochastic  models,  one  suggested  method  of 
verification  is  to  replace  all  probabilistic  effects  of 
the  model  by  a  deterministic  sequence  (constants)  so 
that  hand  computations  of  expected  results  can  be 
compared  to  aodel  output  [2],  Also,  using  the 
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statistical  tests  discussed  in  this  section,  the 
stochastic  model  output  can  be  compared  to  theoretical 
or  historical  data  to  determine  the  existence  of  significant 
difference.*1  Throughout  the  verification  process  the 
analyst  is  checking,  rechecking,  and  debugging  the  simu¬ 
lation  model  until  he  is  assured  that  it  performs  the  way 
he  intends  it. 

Once  the  internal  consistency  of  the  model  is  established, 
the  analyst  can  begin  the  complex  and  everlasting  process 
of  validation.  Face  validity  is  the  surface  impression  of 
a  simulation's  realism  and  is  obtained  by  asking  people 
who  know  the  real  system  [7  ].  In  actuality  it  is  an 
examination  of  the  credibility  or  reasonableness  of  the 
model.  The  "expert"  on  the  real  system  could  compare  real 
world  and  simulation  model  output  to  see  if  he  can 
differentiate  between  them.  However,  the  analyst  must  be 
cautious  since  Turing  has  shown  that  one  can  never  prove 
two  finite  state  systems  f'are  identical  simply  by  comparing 
a  finite  sample  of  input-output  transformations."  [  l  ] 

Some  one  knowledge ible  about  the  real  world  system  could 
also  examine  the  model  in  detail  and  review  its  structure 
and  parameters. 

Parameter  Validity  is  the  agreement  between  the  values  of 
variables  and  parameters  in  the  model  and  their  real 
world  counterparts.  Sensitivity  analysis,  where  variations 
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are  made  in  the  values  of  parameters  to  determine  the 
effect  on  output,  is  one  test  of  parameter  validity. 

Factors  which  are  insensitive  to  variation  do  not  need 
to  be  closely  estimated.  As  part  of  verification,  the 
parameters  should  have  been  accurately  estimated  from 
existing  data.  Various  statistical  tests,  such  as  Chi 
Square,  Kolmogorov- Smirnov,  and  analysis  of  variance,  can 
be  used  to  evaluate  how  well  the  parameters  fit  the  real 
world  process. 

Validity  of  functional  relationships  can  be  partly 
examined  when  face  validity  is  checked.  But  more  detailed 
investigation  is  required  to  determine  if  the  subsystem 
models  are  valid.  Simple  empirical  tests  on  means  and 
variances  or  more  complex  statistical  tests  may  be 
required.  If  a  real  world  system  is  available,  some 
special  data  collection  may  be  arranged  to  acquire  the 
necessary  data  base  for  the  submodels. 

As  we  have  emphasized  on  several  occasions,  the 
decision  maker  should  be  actively  involved  throughout 
the  model  building  process.  In  the  validation  process, 
he  is  able  to  explore  the  interactions  in  the  simulation 
model  and  determine  if  there  is  sufficient  agreement 
between  the  model  output  and  actual  data. 

The  following  quotation  indicates  the  elusiveness  of 
validation  alluded  to  in  the  opening  remarks.  "It  is  never 
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possible  to  completely  validate  a  decision-aiding  model 
since  there  is  never  real  data  about  alternatives  not 
implemented;  this  proolem  is  common  to  any  decision-aiding 
procedure,  not  only  to  simulation.  If  the  decision  makers 
believe  the  model  is  useful  and  use  it,  the  analyst  has 
done  his  job" [7,  p.  206].  Well,  he  has  almost  done  his 
job.  One  significant  point  not  mentioned  in  the  above 
quotation  is  the  analyst's  inherent  obligation  to  advise 
the  decision  maker  of  the  assumptions,  limitations,  or 
simplifications  in  the  moiel.  These  should  be  clearly 
documented  for  the  decision  maker  as  well  as  future  users 
of  the  model. 
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SECTION  Tl 

A  SIMULATION  MODEL  FOR  A  CLASSICAL  PROBLEM 

In  this  section  we  will  develop  a  simulation  model 
to  compare  two  different  strategies  for  the  repair  of  a 
system  of  machines.  In  doing  so,  we  will  utilize  the 
concepts  of  computer  simulation  advanced  in  this i report 
and  apply  some  techniques  of  statistical  analysis. 

Consider  a  single  system  where  a  repairman  is 
assigned  to  service  a  group  of  N  semi-automatic  machines 
which  fail  intermittently.  We  shall  assume  that  the 
failure  rates  of  the  machines  are  identical,  the  event 
of  failure  for  any  machine  is  independent  of  the  state 
of  any  otner  Line,  and  the  time  between  failures  of  a 
singld  machine  and  the  rapair  time  are  distributed  negative 
exponentially  with  respective  means  ^  and  ft  .  The  machine 
layout  is  a  rectangular  pattern  with  two  rows  6f  the  same 
number  of  machines.  The  repairman  requires  time  T  to 
walk  between  any  twm  adjacent  machines. 

Two  strategies  for  tending  the  machines  will  be 
considered: 

(1)  Uni directionally  patrol  the  perimeter  of  the  machines 
repairing  any  failed  machine,  or 

(2)  Repair  the  machines  in  the  order  of  their  failure 
with  the  repairman  located  at  a  central  point  when  no 


machines  are  failed.  In  the  latter  strategy,  we  will 
assume  that  if  a  machine  fails  while  the  repairman  is 
enroute  to  the  idle  location,  the  repairman  continues  to 
the  central  location  before  starting  to  the  machine.  In 
the  patrolling  strategy,  if  a  machine  he  has  just  passed 
fails  before  he  reached  the  next  machine,  the  repairman 
cannot  turn  around  and  return  to  the  failed  machine. 

Development  of  the  Model 

The  general  structure  of  the  simulation  is  shown  in 
the  block  diagram  of  Figure  i'8  .  The  primary  difference 
between  the  models  for  the  two  strategies  is  the  time 
flow  mechanism.  For  the  patrolling  vepairman,  the 
Constant  Increment  Algorithm  from  reference  [23]  was  used 
as  a  basis  for  the  time  flow  mechanism.  The  Cl  algorithm 
was  shown  to  be  very  efficient  for  this  particular  problem. 
For  the  centrally  located  repairman,  a  next  fevent  method 
was  devised. 

The  following  variables  are  defined  for  use  in  the 
simulation: 

TDT  -  total  time  down  for  all  machines 

EFF  -  total  system  efficiency 

DWAIT(J)  -  total  time  machine  J  is  down  awaiting  repair 

DMT(J)  -  individual  time  down  for  machine  J 
Total  system  efficiency  is  further  defined  as  the  total 
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operating  time  for  all  machines  divided  by  total  possible 
operating  time.  The  efficiency  is  then  expressed  as  a 
percentage. 

The  simulation  program  consisted  of  the  following 
subprograms ; 

(1)  MAIN  program  read  the  input  parameters  lambda,  mu, 
number  of  machines,  model  to  be  used,  initial  value  for  the 
random  number  generator,  number  of  repairs  to  be  used  in 
each  block  of  data  (stratified  sampling  is  used),  number 

of  blocks,  indicators  for  determining  if  steady  state 
is  to  be  found  ,and  the  epsilon  to  be  used  in  comparing 
consecutive  blocks.  MAIN  also  calls  the  initialization 
subroutine  and  selects  the  time  flow  mechanism  to  be  used. 

(2)  PATROJ -  is  the  time  flow  mechanism  for  the  patrolling 
repairman  which  determines  the  next  machine  to  be  repaired, 
calle  the  event  processing  subroutine,  and  checks  if  the 
necessary  number  of  replications  have  been  completed. 

(3)  FIFO-  is  the  time  flow  machanisra  for  the  centrally 
located  repairman,  and  its  functions  are  similar  to  PATROL. 

If  ties  for  failure  times  exist,  they  are  broken  by 
choosing  the  closest  machine,  or,  if  they  remain  tied, 

by  lowest  machine  number. 

(4)  MATH-  is  the  event  processing  subroutine  which 
determines  a  service  time  and  new  failure  time,  calculates 
wait  time  and  down  time,  and  updates  the  system  variables. 
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If  steady  state  has  rot  yet  been  reached,  SSTEST  is  called 
to  test  the  values  of  efficiency.  If  the  number  of 
repairs  in  the  block  equals  the  number  required,  the  STAT 
subroutine  is  called  to  store  the  values  of  the  variables. 

(5)  SET-  is  the  initialization  subroutine  which  sets 

the  necessary  variables  to  zero,  initially  assigns  failure  times 
to  all  machines,  and  calls  the  TMATRX  subroutine  for  the 
centrally  located  repairman  for  the* distances  between 
machines.  At  the  end  of  each  block  SET  is  called  again 
and  entered  at  the  multiple  entry  point  ZERO  to  initialize 
the  appropriate  variables, 

(6) =  SSTEST-  is  the  subroutine  to  determine  when  steady 
state  is  reached.  The  last  three  efficiencies  are  stored  and  th 
their  differences  compared  according  to  the  steady  state 
criterion  discussed  later.  The  steady  state  clokk  time, 

last  three  efficiencies ,  and  number  of  repairs  are  printed 
out.  The  appropriate  variables  are  initialized  and  data 
collection  for  simulation  can  begin. 

(7)  TMATRX  -  is  the  subroutine  used  for  the 
centrally  located  repairman  to  calculate  the  distances 
between  machines  and  between  machines  and  the  central 
point. 

(8)  RANG-  is  the  random  number  generator  function. 

The  multipliaative  ccngruential  generator  with  multiplier 
515  and  was  also  tested  with  an  initial  value  of  1907  using 
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the  Pseudorandom  Number  Generator  Test  Package  described 
by  Overstreet,  et  al  I  25].  It  passed  all  tests 
satisfactorily  including  very  good  results  on  the  spectral 
test. 

(S)  F.XTIME  -  is  the  function  to  determine  negative 
exponentially  distributed  random  variates  with  a  given 
mean  failv^e  time, 

(10)  ACOVAR*  is  the  subroutine  calculate  covariance 
term  C  between  consecutive  blocks  of  data. 

(11)  STAT-  is  the  statistics  subroutine  which 
accumulates  values  of  the  variables  at  the  end  of  each 
block,  and  calls  ACOVAR  if  the  necessary  number  of  blocks 
have  been  reached.  At  end  of  the  required  number  of 
replications,  STAT  is  entered  at  the  multiple  entry 
point  FINISH  to  compute  the  final  statistics  using 
stratified  sampling. 

Modular  development  of  the  simulation  program  was 
achieved  without  difficulty.  Virtually  all  subroutines 
are  oommon  to  the  simulation  of  each  strategy.  The 
major  difference  is  in  the  time  flow  mechanisms  which 
were  expressly  chosen  and  tailored  for  the  particular 
strategy.  The  other  difference  is  in  the  initialization 
subroutine  where  the  distances  between  various  locations 
must  be  generated  for  the  centrally  located  repairman. 
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SIMULATION  ANALYSIS  IN  FORTRAN  IV 
MACHINE  INTERFERENCE  PROBLEM 


INPUT  INSTRUCTIONS:  THE  VALUES  OF  THE  FOLLOWING  PARAMETERS 
SHOULD  BE  TYPED  IN  ORDER  ON  A  DATA  CARD  ACCORDING  TO  FORMAT  1. 


C 

C 

C 

C 

C 

C 

C 

C 


C 

C 


C 


LAMB I A  -  (REAL)  MEAN  FAILURE  TIME 
MU  -  (REAL)  MEAN  SERVICE  TIME 
N  -  (INTEGER)  NUMBER  OF  MACHINES 
T  -  (REAL)  TIME  REQUIRED  TO  WALK  BETWEEN  MACHINES 
MODEL  -  (INTEGER)  1  -  PATROLLING  REPAIRMAN 

2  -  CENTRALLY  LOCATED  REPAIRMEN 
INITIAL-  (INTEGER)  INITIAL  VALUE  TO  BE  USED  IN  RANDOM  NUMBER 
GENERATOR 

ITEST  -  (INTEGER)  0  -  IF  NO  TEST  FOR  STEADY  STATE  REQUIRED 
1  -  IF  A  TEST  FOR  STEADY  STATE  DESIRED 
EPSLON  -  (REAL)  VALUE  TO  BE  USED  IN  DETERMINING  STEADY  STATE 
REACHED  WHEN  3  CONSECUTIVE  VALUES  DIFFER  BY  NO  MORE 
THAN  EPSLON 

NBLOCK  -  (INTEGER)  NUMBER  OF  OBSERVATIONS  PER  BLOCK 
IAUTO  -  (INTEGER)  0  -  IF  NO  CALCULATIONS  OF  COVARIANCE  REQUIRED 
N  -  SOME  NONZERO  VALUE  OF  THE  NUMBER  OF  BLOCKS 
TO  BE  USED  IN  CALCULATION  OF  COVARIANCE 
BETWEEN  CONSECUTIVE  BLOCKS.  FOR  EXAMPLE,  5'0. 
NREPET  -  (INTEGER)  NUMBER  OF  BLOCKS  (REPLICATIONS)  OF  SIZE  NBLOCK 
TO  BE  USED  IN  OVERALL  SIMULATION.  ANALYSIS 
-  (REAL)  VALUE  OF  COVARIANCE  BETWEEN  CONSECUTIVE  BLOCKS 


C  DEFINITION  OF  VARIABLES: 

C  CLOCK  -  MASTER  CLOCK  TIME  FROM  BEGINNING %r  SIMULATION  RUN 
C  TSTART-  TIME  WHEN  DATA  COLLECTION  BEGAN  IN  CURRENT  BLOCK  OF  EVENTS 


C  FOR  EACH  BLOCK: 


C  IBLOK  -  BLOCK  COUNTEh 
C  TDT  TOTAL  TIME  DOWN  FOR  ALL  MACHINES 
C  EFF  -  TOTAL  SYSTEM  EFFICIENCY 

C  DWAIT(J)  -  TOTAL  TIME  MACHINE  J  IS  DOWN  AWAITING  REPAIR 
C  BBT(J)-  INDIVIDUAL  TIME  DOWN  FOR  MACHINE  J 
C  DIST(I,J)  -  TIME  REQUIRED  TO  WALK  FROM  MACHINE  I  TO  MACHINE  J 
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c  *****  MAIN  program  ***** 

c 

c 

COMMON/C 1/N,FAIL(18) ,T , CLOCK , TSTART , NBLOCK , NREPET , IBI  OK 

1  /C2/TDT  ,EFF ,DWAIT( 18) ,DMT( 18) ,EF(3 ) , KOUNT , IAUTO, C 

2  /C3/DIST(19,19) 

4  /C4/  LAMBDA ,MU , ITEST , MODEL , EPSLON 

REAL  LAMBDA, MU 

C  tiEAD  INPUT  PARAMETERS 

100  READ(5,1jEND=99)  LAMBDA  ,MU ,  N  ,T  ,MODEL ,  INITAL ,  ITEST , EPSLON , NBLOCK , 
1  IAUTO, NREPET  C 

1  FQHMAT(2F5.^,I5,F5.1,3I5,F5.3,3I5,F5.3,3I5,F5.3) 

WRITE  (6,2)  LAMBDA, MU, N,T, MODEL, INITAL, ITEST, EPSLON, NBLOCK, 

1  IAUTO, NREPET, C 

2  FORMAT  ClNPUT:,,F5.1,5X,F5,l,5X,I2,5X,F4,0,5X,I8,5X, 

1  I2,5X,F3,2,5X,I2,5X,I2,5X,I3,5X,F5.3) 

C  INITIALIZE  RANDOM  NUMBER  GENERATOR  AND  VARIABLES 

CALL  RAN( INITIAL) 

CALL  SET 

C  BRANCH  TO  APPROPRIATE  TIME  FLOW  MECHANISM 
IF(MODEL  -  1)  10,10,20 
10  CALL  PATROL 
GO  TO  100 
20  CALL  FIFO 
GO  TO  100 
99  STOP 
END 


C  *****  TIME  FLOW  MECHANISM  FOR  ***** 

C  *****  PATROLLING  REPAIRMAN  STRATEGY  ***** 

C  SUBROUTINE  PATROL 

C0MM0N/Cl/N,FAIL(18) , T , CLOCK , TSTART , NBLOCK , NREPET , IBLOK 
1  /C2/TDT , EFF , DWA IT ( 18 ) , DMT  ( 18 ) , EF( 3 ) , KOUNT , IAUTO , C 

J  =  0 

TEST  FOR  END  OF  REPLICATIONS 

1  IF ( BLOK, GT, NREPET)  GO  TO  99 

2  J  =  J  +  1 
IF(J.GT.N)  J  =  1 
CLOCK  =  CLOCK  +  T 

IF( CLOCK. LT.FAIL(J) )  GO  TO  2 
C  K  =  J 

CALL  MATH  SUBROUTINE  TO  UPDATE  VALUES  OF  VARIABLES 
CALL  MATH(K) 

GO  TO  1 

C  WHEN  ALL  REPLICATIONS  COMPLETED  CALL  SUBROUTINE  FOR  CALCULATIONS 

OF  FINAL  STATISTICS 
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99  CALL  FINISH 
RETURN 
END 


C  *****  TIME  FLOW  MECHANISM  FOR  ***** 

C  *****  CENTRALLY  LOCATED  REPAIRMAN  ***** 

C 

C 

SUBROUTINE  FIFO 

COMMON  /Cl/  N , FAIL( 18 ) , T , CLCOK , TSTART , NB LOCK , NREPET , IBLOK 

1  /C2/  TDJ.‘,EFF,EWAIT(l8),DMr(l8),EF(3),  KOUN^IAUTO^ 

2  /C3/  DIST(19,19) 

N1  *  N  +  1 

K  =  N1 

C  TEST  FOR  END  OF  REPLICATIONS 

1  IF(IBL0K,GT, NREPET)  GO  TO  99 

DETERMINE  MINIMUM  FAILURE  TIME,  IF  A  TIE,  USE  CLOSEST  ONE, 

2  EMIN  =  2. t*30 
DO  10  J  *  1,N 

IF( FAIL( J )  -EMIN)  5-3.10 

3  IF(DIST(K,JMIN).LE.DIST(K,J))  GO  TO  10 
5  FMIN  =  FAIL(J) 

JMIN  =  J 
10  CONTINUE 
C 

C  IF  FAILED  ALREADY,  GO  TO  THAT  MACHINE:  OTHERWISE  GO  TO  CENTRAL 

C  LOCATION. 

IF( CLOCK. GE. FMIN)  GOTO  20 
CLOCK  =  CLOCK  +  DIST  (K,Nl) 

K  =  N1 

IF( CLOCK. GE. FMIN)  GO  TO  20 
CLOCK  =  FMIN 
C 

C  UPDATE  CLOCK 

20  CLOCK  =  CLOCK  +  DIST (K, JMIN) 

K  =  JMIN 
C 

C  UPDATE  VARIABLES 
CALL  MATH  (K) 

GO  TO  1 

99  CALL  FINISH 
RETURN 
END 
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*****  EVENT  PROCESSING  ROUTINE  ***** 


SUBROUTINE  MATH(K) 

C0MK0N/C1/N ,  FAIL(  18),T,  CLOCK,  TSTART ,  NBLOGK ,  RTuSPET ,  IBLOK 
1  /C2/TDT ,EFF ,DWAIT( 18) ,EF(3), KOUNT , I AUTO , C 

4  / C4/LAMBDA ,MU , ITEST , MODEL , EPSLON 

REAL, LAMBDA,  MU 
DETERMINE  SERVICE  TIME 
ST  =  EXTIME(MU) 

DETERMINE  TIME  MACHINE  WAITED  FOR  REPAIRMAN 
WAIT  =  CLOCK  -  FAIL(K) 

DETERMINE  DOWNTIME  AND  UPDATE  TOTAL  DOWNTIME  FOR  ALL  MACHINES 
DOWN  «  WAIT  +  ST 
TDT  =  TDT  +  DOWN 

UPDATE  CLOCK  AND  VARIABLES 
CLOCK  =  CLOCK  +  ST 

EFF  =  (N*( CLOCK  -  TSTART)  -  TDT)  /  (N*( CLOCK  -  TSTART) 

1  )*100. 

DWAIT(K)  =  DWAIT(K)  +  WAIT 
DMT(K)  =■  DMT(K)  +  DOWN 

CALCULATE  NEW  FAILURE  TIME 
FT  =  EXTIME(LPMBDA) 

FAIL(K)  =  CLOCK  +  FT 

INCREMENT  COUNTER  FOR  NUMBER  OF  REPAIRS 

KOUNT  =  KOUNT  +  1 

IF  REQUIRED,  CALL  STEADY  STATE  TEST  SUBROUTINE 
IF  (ITEST.  LE-  0)  GO  TO  99 

IF  STEADY  STATE  HAS  BEEN  REACHED,  CALL  STATISTICS  SUBROUTINE 
CALL  SSTEST 

99  IF(  ITEST- LE  0.  AND.  KOUNT.  GE.  NBLOCK)  CALL  STAT 
RETURN 
END 


******  VARIABLE  INTIALIZATION  ****^s* 
******  SUBROUTINE  ****** 


SUBROUTINE  SET 

COMMON/ Cl/N , FAIL( 18 ) ,T , CLOCK , TSTART , NBLOCK , NREPET , IBLOK 
1  /C2/TDT , EFF , DWAIT( 18 ) , DMT ( 18 ) , EF ( 3 ) , KOUNT , I AUTO , C 

4  /C4/  LAMBDA, MU, ITEST, MODEL, EPSLON 

REAL  LAMBDA,  MU 
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C  INITIAL  ENTRY  POINT  TO  INITIALIZE  VARIABLES 

CLOCK  =  0 
IBLOK  =  1 
DO  10  K  =  1  N 

10  FAIL(K)  =  EXTIME  (LAMBDA) 

C  IF  MODEL  2  USED  CALL  SUBROUTINE  TO  CALCULATE  VARIOUS  DISTANCES 
C  BETWEEN  MACHINES 

IF(MODEL.  LE.l)  GO  TO  11 
CALL  TMATRX(N.T) 

C 

C 

C  SUBSEQUENT  ENTRY  POINT  TO  ZERO  VARIABLES  IN  EACH  BLOCK 
ENTRY  ZERO 
C 

11  TSTART  =  CLOCK 
EF(1)  =  199. 

EF(2)  =  299. 

EF( 3 )  =  399. 

KOUNT  =  0 
IDT  =  0. 

EFF  =  0 . 

DO  15  K=  1,N 

dwait(k)=o. 

15  DMT(K)  =  0. 

RETURN 

END 

C  *****  STEADY  STATE  TEST  ***** 

C 

C 

SUBROUTINE  SSTEST 

COMMON/ Cl/N, FAIL( 18 ) ,T , CLOCK , TSTART , NBLOCK , NREPET ,  IBLOK 
1  / C2/TDT , EFF , BWAIT ( 18  ) , DMT ( 18 ) ,  EF(3),  KOUNT, IAUTO,C 

4  /C4/LAMBDA,MU  ITEST  MODEL  EPSLON 
DATA  1/0/ 

C  STORE  LAST  THREE  EFFICIENCIES 

1  =  1  +  1 
IF(I.GT.  3)1=1 
EF(I)  =  EFF 
C 

C  IF  ALL  OF  LAST  THREE  EFFICIENCIES  DIFFER  BY  EPSLON  OR  LESS, 

C  STEADY  STATE  REACHED. 

IF  (ABS(EF(l)-EF(2)).LE.  EPSLON.  AND.  ABS(EF(2)  -  EF(3)) 

1  .  LE. EPSLON.  AND. ABS(EF(1)  -  EF(3)).  LE.  EPSLON) 

2  GO  TO  90 
GO  TO  99 

90  ITEST  =  -1 
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C  PRINT  INFORMATION  ON  STEADY  STATE 

WRITE(6,9l)-CLOCK,  EF(l),  EF(2),  EF(3),  KOUNT 
91  F0RMAI( 1H0 ' STEADY  STATE  AT  ' , F15 .3 , ' CLOCK  TIME ' , 5X, ' EFFICIENCIES 
1  :',3(F8,3,3X),' AND' ,110, 'REPAIRS*) 

CALL  ZERO 
99  RETURN 
EHD 


SUBROUTINE  TMATRX(N,T) 

*****  SUBROUTINE  FOR  CALCULATION  OF  ***** 
*****  DISTANCES  BETWEEN  MACHINES  ***** 


C0MM0N/C3/DIST(19,19) 

CENTRAL  LOCATION  IS  POSITION  N+l.  WHERE  N  IS  NUMBER  OF  MACHINES 
N1  =  N  +  1 
DO  10  I  =  1,  N 
DO  10  J  =  I,  N 
K  =  N/2 

TEST  FOR  MACHINES  IN  SAME  ROW:  OTHERWISE  USE  DISTANCE  FORMULA 
IF(I.LE.K.AND.J.GT.K)  GO  TO  5 
C  =  J  -  I 
GO  TO  6 

5  C  =  SQRT((N  =  1  -  J  -  I)**2  +  l) 

MULTIPLY  DISTANCE  BY  TIME  BETWEEN  MACHINES 
6  DI3T(I,J)  =  C  *  T 
10  DIST(J,I)  =DIST(I,J) 

CALCULATE  DISTANCE  TO  CENTRAL  LOCATION 
DO  15  J  =  1,K 

DIST(N1,J)  =  T*SQRT((J-N/4-l)**2+,25) 

M  =  J  +  K 

15  DIST(N1,M)  =  DIST(N1,J) 

D)  20  J  =  1,N 

20  DIST(J,N1)  =  DIST(N1,J) 

RETURN 

END 


*****  SUBROUTINE  FOR  COVARIANCE  ***** 
*****  BETWEEN  CONSECUTIVE  BLOCKS  ***** 


SUBROUTINE  ACOVAR 

COMMON/ C2/TDT ,  EFF,  DWAIT(l8),  DMT(lB),  EF(3),  KOUNT,  IAUTO,C 
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1  /C5/BTDT(100),  BEFF(lOO),  BDWAIT( 100,18),  BDMT(lOO,18) 

DIMENSION  X(  100 ) , C0VAR( 100 ) 

XBAR  =  0. 

C  COPY  VALUES  OF  VARIABLE  IN  TO  X  ARRAY. 

DO  5  I=1,IAUT0 
10  XBAR  =  SBAR  +  X(l) 

AVG  =  0. 

SBAR  =  SBAR/IAUTO 
N  =  IAUTO  -  1 
C 

C  CALCULATE  COVARIANCE  TERM  C  BETWEEN  CONSECUTIVE  BLOCKS 
DO  15  I  =  1,N 

COVAR(I)  =  (X(I)  -  XBAR)*(X( I*-l)-XBAR) 

15  AVG  =  AVG  +  COVAR(I) 

AVG  =  AVG/N 

C  PRINT  CALCULATIONS 

WRITE(6,20)  (COVAR(I),  1=1, N),  AVG 
20  FORMAT(  ( 1H0, 10(  IX.,  Fll .  8 ) )  ) 

RETURN 

END 


*****  STATISTICAL  CALCULATIONS  SUBROUTINE  ***** 


SUBROUTINE  ST AT 

COMMON  /Cl/N , FAIL( 18 ) , T , CLOCK , TSTART , NBLOCK , NREPET , IBLOK 

1  /C2/TDT,EFF, DWAIT( 18) , DMT( 18) , EF(3 ) , KOUNT , IAUTO ,  C 

2  /C3/DIST(19,19) 

4  / C4/LAMBDA ,MU , ITEST , MODEL, EPSLON 

5  /C5/BTDT( 100 ) , BEFF( 100 ) , BDWAIT( 100 , 18 ) , BDMT ( 100 , 18 ) 
REAL  LAMBDA  MU 

DIMENS ION  AWAITL( 18 ) , ADMTl( 18 ) , AWA IT2 ( 18 ) , ADMT2 ( 18 ) , AWAIT( 18 ) , 

1  VWAIT( 18 ) , SWAIT( 18 ) , VWAIT1( 18) , VWAIT2( 18 ) ,  VDMT(l8), 

2  VDMT  (18),  VDMT2(18),  SDMT(lS),  ADMT(l8) 

C 

C  ASSIGN  VALUES  TO  -VARIABLES  FROM  EACH  BLOCK 
BTDT( IBLOK)  =  TDT 
DEFF(  IBLOK)  =  EFF 
DO  10  1=1, N 

BDWAIT( IBLOK  i)  =  DWAIT(l) 

10  BDMT  (IBLOK,  I)  =  DMT(l) 


143 


c 

C  CALL  COVARIANCE  SUBROUTINE  IF  APPROPRIATE 

IF  ( IAUTO.GT.O. AND. IBLOK. GE.TAUTO)  CALL  ACOVAR 
C  INCREMENT  BLOCK  COUNTER 

IBLOK  =  IBLOK  +  1 
C 

C  ZERO  VARIABLES  FOR  NEXT  BLOCK 

CALL  ZERO 
RETURN 
C 

C  ENTRY  POINT  FOR  FINAL  CALCULATIONS 

ENTRY  FINISH 
ATDT1  =  0. 

ATDT2  =  0. 

AEFF1  =  0. 

AEFF2  -  0. 

VTDTi  =  0. 

VTDT2  =  0. 

VEFF1  =  0. 

VEFF2  =  0. 

IBLOK  =  IBLOK  -  1 
IF(IAUTO.GT.O)  RETURN 
IHALF  =  IBLOK/2 
C 

C  CALCULATE  MEAN  TOTAL  TIME  DOWN  AND  EFFICIENCY  FOR  EVEN  AND  ODD 
C  BLOCKS  AND  OVERALL  MEAN 
DO  20  1=1,  IBLOK, 2 
ATDT1  =  ATDT1  +  BTDT(l) 

ATDT2  =  ATDT2  +  BTDT  (1+1) 

AEFF1  =  AEFF1  +  BEFF  (I  +  l) 

20  AEFF2  =  AEFF2  +  BEFF  (i+l) 

ATDT1  =  ATDTl/ IHALF 
ATDT2  =  ATDT2/IHALF 
AEFF1  =  AEFFl/lHALF 
AEFF2  =  AEFF2/IHALF 
ATDT  =  (ATDTl  +  ATDT2)/2, 

AEFF  =  (AEFF1  =  AEFF2)/2. 

C 

C  COMPUTE  VARIANCES  FOR  EVEN  AND  ODD  BLOCKS  AND  OVERALL  VARIANCES 
DO  30  1=1, IBLOK, 2 

VTDTI  =  VTDTI  +  (BTDT(l)  -  ATDTl)**2 
VTDT2  =  VTDT2  +  (BTDT(+l)  -  ATDT2)**2 
VEFF1  =  VEFF1  +  (BEFF(l)  -  AEFFl)**2 
30  VEFF2  =  VEFF2  +  (BEFF(l+l)  -  AEFFl)^2 
\TDT2  =  VTDT2/IHALF 
VTDTI  =  VTDTi/ IHALF 
VEFF1  =  VEFF1/ IHALF 
VEFF2  =  VEFF2/IHALF 

VEFF  =  .25-;KVEFFl  +  VEFF2)  +C^2.  IB  LOK-1 ) /lBL0K^2 
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VTDT  =.  25*(VTDT.+VTDT2)  +0*2.  *(  BLOK-1  )/lBL0K**2 
C  DETERMINE  STANDARD  DEVIATIONS 
SEFF  =  SQRT(VEFF) 

STDT  =  SQRT(VTDT) 

CALCULATE  MEAN  WAITING  TIME  AND  TOTAL  TIME  DOWN  FOR  EACH  MACHINE 
FOR  EVEN  AND  ODD  BLOCKS  AND  OVERALL 
DO  40  J  =  1,N 
AWAITI(J)  =  0. 

AWAIT2(J)  =  0. 

ADMTI(J)  =  0. 

ADMT2( J)  -  0. 

DO  40  I  =  1,  BLOK, 2 
AWAITl(J)  =  AWAITl(J)  +  BDWAIT(I,J) 

AWAIT2(J)  =  AWAIT2(J)  +  BDWAIT(I+1, j) 

ADMTI(J)  =  ADMTl(J)  +  BDMT(I, '  J) 

40  ADMT2(J)  =  AEMT2( J)  +  BDMT(I+1,J) 

DO  50  J=1,N 

AWAITl(J)  -  AWAITl( J)/IHALF 
AWAIT2(J)  =  AWAIT2(J)/IHALF 
ADMTI(J)  =  ADMTl(J)/lHALF 
ADMT2(J)  =  ADMT2(J)/IHALF 
AWAIT(J)  =  (AWAITI(J)AHALF 
50  ADMT(J)  =  (ADMTl(J)  +  ADMT2(j))/2 
C 

C  COMPUTE  VARIANCES  AND  STANDARD  DEVIATIONS 
DO  60  <1=1,  N 
DO  60  1=1, BLOK, 2 

VWAITI(J)  =  VWAITI(J)  +(BDWAIT(I,J)  -  AWAITl( J))**2 
VWAIT2(J)  =  VWAIT2(J)  +(BDWAIT(I+1, j)  -  AWAIT2( J))**2 
VDMTl(J)  =  VDMTl(J)  +  ( B  DMT  (l,j)-ADMTl(<J))  *-*2 
60  VDMT2(J)  =  VDMT2(J)  +  (BDMT(l+l, j)-ADMT2( J) )**2 
DO  70  J=1N 

VWAITI(J)  =  VWAITl( J )/lHALF 
VWAIT2(J)  =  VWAIT2(J)/IHALF 
VDMTI(J)  =  VDMTl( J )/ IH ALF 
VDMT2(  j)  =  VDMT2(j')/lHALF 

VWAIT(J)  =  .25*(VWAITl(J)  +  VWAIT2(  J))+C*2,*(BL0K-l)/BL0K**2 
SWAIT(J)  =  SQRT(VWAIT( J) ) 

VDMT(J)  =  .25*(VDMT1(J)  +  VEMT2(  J))+O>2(lBL0K-l)/lBL0K^2 
70  SDMT(J)  =  SQRT'( VDMT( j) ) 

C 

C  PRINT  OUTPUT  DATA  AND  STATISTICS 

WRITE  (6,80)  MODEL, LAMBDA  MU,  ATDT, STDT, AEFF, SEFF, 

1  CLOCK,  IBLOKjNBLOCK 

80  FORMAT  (1H)/1H0,38X, 'OUTPUT  VALUES  FOR  MODEL' ,12, 'WITH  LAMBDA', 

1  '=' ,FY,1, 'AND  MU  =  ',F5.l///lHO,45X,'MEAN',15X, 

2  ' STANDARD  DEVIATION ' /lHO, ' TOTAL  TIME  DOWN ' , 

3  'FOR  ALL  MACHINES ',8X,F15, 3, 13X,F11.3/1H0, 'TOTAL', 
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4  'SYSTEM  EFFICIENCY', 17X,F15.3,6X,F15.3,L3X/lHO,/lHO, 

5  'FINAL  CLOCK  TIME',  13X,F25.3/  1H0,  NUMBER  OF  REPLICATIONS', 

6  18X, 110/  1H0, 'NUMBER  OF  REPAIRS  PER  REPLICATION' ,7X,  ILO) 
PRINT  90 

90  F0RMAT(1H0, 'MACHINE',  15X, 'TIME  DOWN' ,23X, 'TIME  DOWN  AWAITING' 
1  'REPAIR'/l8X,  "MEAN', 14X, 'STD  DEV' ,15X, 'MEAN' ,L3X, 'STD  DEV'/ 
DO  91  J  =  1,N 

91  WRITE(6,92)  J,  AEMT(J),SDMT(J),AWAIT(J),SWAIT(J) 

92  FORMAT ( 3X , 13 , 2 ( 10X , F10 . 3 , 9X , F10 , 3 ) ) 

RETURN 

END 

*****  RANDOM  NUMBER  GENERATOR  ***** 


FUNCTION  RAND(DUMMY) 

I  =  IA  *  I 
FLD(0,1  I)  =  0 
RANG  - 1  /2. 0**35 
RETURN 

INITIALIZATION  OF  MULTIPLIER  AND  INITIAL  VAI3JE 
ENTRY  RAN(  INITAL) 

I  =  INITAL 
IA  =  5  **  15 
RETURN 
END 


*****  EXPONENTIAL  INTEROCCURRENCE  ***** 
*****  TIME  FUNCTION  ***** 


FUNCTION  EXTIME(MEAN) 

REAL  MEW 

EXTIME  +  -  MEAN*  ALOG(RANG(RANUM)) 
10  RETURN 
END 
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Convergence  to  Steady  State 

The  simulation  for  both  strategies  began  with  all 
machines  working.  The  repairman  began  at  machine  one  in 
the  patrolling  case  and  at  the  idle  position  for  the 
centrally  located  strategy. 

The  total  system  efficiency  was  chosen  as  a  variable 
which  would  be  indicative  of  the  model's  convergence  to 
steady  state.  The  steady  state  criterion  used  was;  Steady 
state  was  reached  when  three  consecutive  values  differed 
by  no  more  than  .01.  Table  18  summarizes  the  steady 
state  times  and  number  of  repairs  for  each  set  of 
parameters. 


TABLE  18  CONVERGENCE  TO  STEADY  STATE 

-30 


N 

X 

T 

PATROLLING 

CENTRALLY  LOCATED 

CLOCK 

REP 

EFF 

CLOCK 

REP 

EFF 

6 

600 

10 

17822 

177 

87.22 

6579 

62 

91.58 

300 

10 

17332 

285 

75.05 

11138 

195 

77.79 

600 

100 

50738 

306 

55.96 

32921 

194 

53.32 

300 

100 

56784 

411 

36.47 

64932 

572 

28.72 

18 

600 

10 

20019 

404 

63.38 

32068 

518 

54.49 

300 

10 

9785 

240 

59.16 

36497 

572 

27.34 

600 

100 

36665 

280 

26.15 

262554 

698 

10.00 

300 

100 

31145 

240 

14.80 

244965 

397 

6.58 
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Statistical  Measurements 

Statistical  Measurements  were  obtained  on  total  system 
efficiency,  total  time  down  for  all  machines,  individual 
time  down  for  each  machine,  and  total  time  each  machine 
was  down  awaiting  repair.  Total  system  efficiency  was 
considered  the  variable  of  primary  interest  and  the 
estimates  for  each  strategy  are  given  in  Table  22.  Sample 
model  printouts  are  shown  in  Figures  29  and  30. 

The  first  concern  in  data  collection  was  the  sample 
size,  i.e.,  the  number  of  repairs  contained  in  each 
replication.  Several  techniques  are  available  for 
determining  length  of  the  sample  needed  for  an  adequate 
statistical  analysis.  Fishman  [  9]  proposed  a  method  for 
finding  the  number  of  equivalent  independent  observations 
using  the  spectral  density  function.  A  detailed  evaluation 
of  the  mathematics  involved,  length  of  computer  program 
required,  and  run  time  of  the  program  led  to  the  conclusion 
that  for  the  length  of  the  simulation  lequired  for  this 
study,  the  spectral  density  approach  was  not  appropriate. 

The  determination  of  a  sample  size  would  involve  calculations 
for  lags  up  to  500  or  more  which  alone  would  be  several 
times  the  total  run  time  of  the  entire  simulation  for  all 
sets  of  parameters.  For  simulation  models  where  each 
replication  is  very  lengthy  and  costly,  the  Fishman  technique 
would  be  valuable. 
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Figure  29.  Sample  Printou+  '/or  Centrally  Located  Repai: 
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Figure  30.  Sample  Printout  For  Patrolling  Repai 


As  a  first  step  in  designing  a  method  of  sample  size 
selection,  the  values  of  the  autocorrelation  function  were 
examined  for  lags  of  100  to  600  for  the  patrolling  strategy 
with  Y"  600,yw«30,  N«6,  and  T-10,  Some  values  art;  displayed 
in  Table  IS  below. 


TABLE  19.  VALUES  OF  AUTOCORRELATION  FUNCTION 


LAG 

AUTOCORRELATION 

100 

.868 

150 

-.062 

200 

-.047 

250 

-.098 

300 

-.024 

350 

-.017 

400 

.057 

450 

.055 

500 

.048 

550 

-.028 

600 

-.045 

It  is  clear  that  the  autocorrelation  function  does  not 

asymptotically  approach  zero  but  exhibits  some  cyclic 

behavior.  The  autocorrelation  closest  to  zero  was  -.0002 

which  occurred  for  a  lag  of  533.  Since  in  the  simulation 

model,  we  desire  independence  between  replications  (blocks 

of  observationa  of  some  fixed  size)  or  at  least  a  small 

covariance,  the  ACOVAR  subroutine  was  used  to  determine 

covariance  between  consecutive  blocks.  From  a  preliminary 

check  of  various  values,  the  covariance  between  nonadjacent 

blocks  was  always  several  orders  of  magnitude  less  than  the 

covariance  becween  adjacent  blocks.  Obviously,  if  the 

covariance  between  adjacent  blocks  was  made  sufficiently 

small,  the  covariance  for  nonadjacent  blocks  could  be 
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considered  zero.  Some  values  of  covariance  for  various 
sample  sizes  for  A  «600,  JA-30,  n»6,  T*10  are  presented 
in  Table  20. 

TABLE  20  COVARIANCE  FOR  VARIOUS  SAMPLE  SIZES 

SO  BLOCKS 


SAMPLE  SIZE  EACH  BLOCK  PATROLLING  CENTRALLY  LOCATED 


100 

-.105 

-.191 

200 

.041 

.197 

300 

-.027 

.062 

400 

.042 

.021 

500 

.042 

-.038 

600 

-.088 

-.015 

700 

-.021 

-.033 

800 

-.005 

-.028 

After  a  comparison  of  various  combinations  of  sample  sizes 

for  all  the  parameter  sets,  a  size  of  600  repairs  was  chosen 

as  producing  the  lowest  covariances  for  all  sets.  By  using 

the  same  sample  size  for  all  simulation  runs,  one  is  able 

to  make  a  more  valid  comparison.  The  covariances  of  adjacent 

blocks  for  a  sample  size  of  600  and  six  machines  for  various 

parameters  are  displayed  in  Table  21. 

TABLE  21  COVARIANCE  ESTIMATES 
N-6,  p*  30,  50  BLOCKS 

X  T  PATROLLING  CENTRALLY  LOCATED 


600 

10 

-.089 

-.015 

300 

10 

-.026 

-.113 

600 

100 

.095 

.255 

300 

100 

.196 

-.028 

After  steady  state  was  reached,  the  simulation  run  was 
broken  into  replications  (blocks).  Each  block  contained 
600  repairs  and  the  entire  run  was  comprised  of  50  blocks. 
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Fifty  was  chosen  so  that  large  sample  statistics  could  be 

later  applied  in  a  comparison  of  the  results  of  each 

strategy.  The  blocks  were  classified  as  odd  or  even  and 

a  value  of  each  variable  stored  for  every  block.  A  mean 

and  variance  was  then  calculated  for  the  odd  blocks  and 

for  the  even  blocks.  A  pooled  nean  and  variance  was 

determined.  The  pooled  variance  was  calculated  from 

Var(X  ♦  Y)  -  1/4  (Var  X  ♦  Var  Y  ♦  Cov(X»Y)  ) 

2 

-  1/4  (Var  X  ♦  Var  Y)  ♦  1/4  Cov(:;,Y) 

-  1/4  (Var  X  ♦  Var  Y)  ♦ 2  £  £  (X-X)  (yY) 

n*  i  4 J 


-  1/4  (Var  X  +  Var  Y)  ♦  2(n-l)  C 


where  C  is  the  covariance  term  between  adjacent  blocks 
previously  determined.  The  covariance  between  nonadjacent 
blocks  was  assumed  to  be  zero. 


The  estimates  of  the  total  system  efficiency  for  each 
set  of  parameters  and  strategy  are  shown  in  Table  22  below. 


TABLE  22  SUMMARY  OF  EFFICIENCY  ESTIMATES 

/*-  30 

N  A  T  PATROLLING  CENTRALLY  LOCATED 


MEAN 

STD  DEV 

MEAN 

STD  DEV 

6 

600 

10 

88,58 

.54 

91.68 

.45 

300 

10 

75.76 

.98 

79.71 

1.12 

5f  ■» 

100 

57,05 

.98 

56.40 

1,71 

3v 

100 

35,74 

.85 

29,34 

,94 

18 

600 

10 

67,08 

1,33 

52,95 

1.71 

300 

10 

40, OS 

1.56 

26,52 

,96 

600 

100 

25.20 

,77 

9,25 

,44 

300 

100 

12,82 

,47 

4,59 

.32 
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Reliability  of  Estimates 


Since  we  have  purposely  chosen  the  number  of  replica¬ 
tions  to  be  sufficiently  large  to  use  large  sample  (normal) 
statistics,  we  invoke  the  central  limit  theorem  and  determine 
confidence  intervals  for  the  estimates  of  total  system 
efficiency.  Thus,  from  the  normal  distribution  we  have 


*Z 


dll 


x  -  )U 

s/r?r 


*/i 


and  P(*'z«/l5  C  /A  <  *  +  ***£; 

'In 


where  S  is  sample  standard  deviation 

X  is  the  estimate  of  total  system  efficiency 
n  is  the  sample  size,  in  our  case,  50 
Assuming  a  significance  level  of  .05,  we  have  the  confidence 
intervals  shown  in  Table  23. 


TAFLE  23.  95%  CONFIDENCE  INTERVALS 

FOR  ESTIMATES  OF  TOTAL  SYSTEM  EFFICIENCY 

f*  «  30 


N 

6 


Patrolling  Centrally  Located 


Lower 

Upper 

Lower 

Upper 

T 

Bound 

Bound 

Bound 

Bound 

600 

10 

88,43 

88,73 

91.56 

91.81 

300 

10 

75.49 

76. 03 

79.39 

80.02 

600 

100 

56.78 

57,32 

55.93 

56.88 

300 

100 

35.50 

35.97 

29.08 

29.60 

600 

10 

66, 71 

67.45 

52.47 

53,42 

300 

10 

39.65 

40.52 

26.25 

26.78 

600 

100 

24.99 

25.42 

9.13 

9.37 

300 

100 

12.69 

12.95 

4.50 

4.68 
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Comparison  of  Strategies 

Since  the  samples  were  large  enough  for  the  central 
limit  theorem  to  be  invoked,  one  can  test  the  difference 
between  the  estimates  of  total  system  efficiency  foe  the 
two  strategies  for  a  given  set  of  parameters  and  number 
of  machines.  Assume  the  samples  are  from  two  normal 
populations  where  the  patrolling  strategy  has  mean  ui  and 
the  central  location,  mean  U2  with  variances  unknown.  The 
statistic 


where  Xj  is  the  estimate  from  the  patrolling  strategy, 

X?  is  the  estimate  from  the  centrally  located 
strategy, 


S,  and  S 
deviations 


ars  the  respective  sample  standard 


■  n£  *  50 


has  a  standard  normal  distribution.  The  null  hypothesis 
that  -  u2  "  0  is  teste  against  the  alternative  that 
ui  -  u 2  i  0.  The  critical  region  is  |  l\iZ </2.  Foro<«  .05, 
one  must  have  |z|  £  1.96.  The  statistic  for  each  case  is 
presented  in  Table  24,  The  ntll  hypothesis  is  rejected 
in  all  cases.  There  is  a  significant  difference  in  the 
means  for  the  two  strategies  in  all  cases, 
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In  the  first  two  cases  with  six  machines,  one  can  see 
that  the  centrally  located  strategy  produced  a  higher  total 
system  efficiency.  For  the  remaining  cases,  the  patrolling 
strategy  had  a  higher  total  system  efficiency.  For  the 
remaining  cases,  the  patrolling  strategy  had  a  higher  total 
system  efficiency, 

TABLE  24.  TEST  OF  DIFFERENCE  OF  MEANS 
TOTAL  SYSTEM  EFFICIENCY 
yU-  30 

N  A  T  Z 


6 

600 

10 

-31.2 

300 

10 

-18.7 

600 

100 

2.3 

300 

100 

30.2 

18 

600 

10 

46.0 

300 

10 

53.4 

600 

100 

402. 

300 

100 

321. 

In  a  qualitative  sense  one  can  observe  from  Table  24 
that:  (1)  As  the  walk  time  between  machines  increased  with 

other  parameters  constant,  the  total  system  efficiency  for  both 
strategies  decreased  and  the  centrally  located  strategy  had 
a  larger  decrease. 

(2)  As  the  failure  rate  increased  with  other  parameters 
constant,  the  total  system  efficiency  decreased, 

(3)  For  given  ^  and  T,  as  the  number  of  machines 
increased,  the  total  system  efficiency  decreased. 
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SECTION  VII 


COMPUTER  SIMULATION  LANGUAGES 

The  discussion  of  simulation  techniques  throughout 
the  major  part  of  this  report  is  essentially  independent 
of  the  programming  language  used,  However,  at  some  point 
the  analyst  must  decide  what  language  to  use:  either  a 
general  purpose  language  (GPL)  or  a  simulation  program¬ 
ming  language  (SPL). 

The  purpose  of  this  section  is  to  highlight  the  dif¬ 
ferences  between  GPL  and  SPL,  and  examine  the  impact  of 
using  a  particular  language  on  the  model.  This  section 
is  not  a  detailed  discussion  of  programming  languages. 

The  reader  is  directed  to  the  appropriate  reference  cited 
in  the  discussion  for  an  in  depth  treatment  of  features 
of  each  language - 

Regardless  of  the  language  used,  certain  basic  func¬ 
tions  are  required: 

•  Pseudorandom  number  generation 

•  Creation  of  random  variates 

•  A  time-flow  mechanism 

•  Data  storage 

•  Statistical  analysis 

•  Data  output 

•  Diagnostic  information  on  logic  and  other  errors. 
In  the  following  discussion,  we  will  comment  on  each  of 


these  features. 
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General  Purpose  Languages 

An  analyst  with  access  to  a  computer  system  would 
have  some  general  purpose  language  available  for  use  in 
simulation.  The  language  may  be  FORTRAN,  ALGOL,  COBOL, 
or  PL/1.  Let's  examine  the  advantages  and  disadvantages 
of  using  a  GPL  for  simulation  rather  than  an  SPI.. 

The  chief  advantage  in  using  a  GPL  is  flexibility. 

The  analyst  is  free  to  write  any  subroutine  he  believes 
is  required  for  the  simulation.  An  SPL  is  usually  rigidly 
structured  so  that  the  analyst  must  use  what  is  built  into 
the  language.  He  may  or  may  not  be  able  to  conveniently 
add  his  own  subroutines.  Also,  an  SPL  may  imply  that  the 
problem  be  expressed  in  a  particular  form,  such  as  flow 
through  a  block  diagram  as  in  GPSS. 

A  second  advantage  of  a  GPL  is  universal  availability. 
FORTRAN  compilers,  for  example,  exist  for  almost  every 
make  and  size  computer.  This  permits  greater  compatibility 
between  machines  for  simulations  which  mr”-  be  transferred 
between  computer  installations. 

In  terms  of  the  basic  requirements,  the  analyst  must 
provide  a  pseudorandom  number  generator  and  a  means  of 
creating  random  variates.  This  is  not  really  a  disad¬ 
vantage  since  by  providing  the  random  number  generator 
the  analyst  knows  exactly  what  he  has.  He  probably  has 
already  evaluated  the  generator's  behavior  and  has  demon- 
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strated  its  acceptability.  By  using  *«  SPL  pseudorandom 
number  generator,  the  analyst  is  probably  assuming  that 
the  generator  has  been  tested  and  demonstrate*  randomness 
in  a  statistical  sense.  This  unfortunately  could  be  an 
invalid  assumption  leading  to  cataclysmic  results.  Often 
the  random  number  generator  in  an  SPL  cannot  be  directly 
accessed  for  statistical  testing  and  the  user  is  left 
unsure  of  its  behavior. 

In  a  GPL,  the  analyst  is  free  to  choose  the  time  flow 
mechanism  which  he  believes  is  most  appropriate  for  that 
particular  simulation.  An  SPL  may  force  the  analyst  into 
using  a  particular  time  flow  mechanism, usually  next  event. 

The  analyst  must  also  provide  a  means  jf  storing  and 
manipulating  data  for  the  output  report.  But  he  is  free 
to  use  graphical  or  tabular  formats  in  whatever  represen¬ 
tation  is  suitable.  This  flexibility  is  not  always  present 
in  SPLs. 

The  chief  disadvantage  of  a  GPL  is  in  diagnostic 
information  on  the  program.  The  GPL  compiler  will  detect 
syntax  errors  in  language  usage  but  does  not  have  the 
elaborate  logic  checks  present  in  SPL. 
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Simulation  Programming  Languages 

The  most  popular  Simulation  Programming  Languages 
are  GASP,  GPSS,  SIMSCRIPT,  and  SIMULA.  Table  25  dis¬ 
plays  the  current  availability  of  these  languages. 


Table  25.  The  Availability  of  Simulation 
Programming  Languages 


Computer  Manufacturer 

GPSS 

SIMSCRIPT  II 

SIMULA 

Burroughs 

X 

X 

X 

CDC 

X 

X 

Honeywell 

X 

X 

IBM 

X 

X 

UNIVAC 

X 

X 

X 

The  discussion  of  each  language  will  include  its  advantages 
and  disadvantages  as  well  as  how  it  views  the  world. 

GASP 

The  General  Activity  Simulation  Program  (GASP)  is 
unlike  other  SPLs  since  it  is  written  in  FORTRAN  and  can 
be  used  on  any  machine  with  a  FORTRAN  compiler .  GASP  is 
a  set  of  FORTRAN  subroutines  and  functions  which  are 
linked  by  a  nuiin  program  called  the  GASP  EXECUTIVE.  The 
main  program  starts  the  simulation,  enables  sequencing  of 
time,  monitors  intermediate  results,  and  initiates  printout 
of  results. 
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The  chief  advantages  of  GASP  are  that  it  is  not 
necessary  to  learn  a  new  language  or  obtain  a  new  compiler, 
it  employs  the  modular  approach  used  in  simulation,  and  it 
has  machine  independence. 

The  disadvantages  relate  to  GASP's  basis  in  FORTRAN. 
Problems  are  encountered  in  input/output  formating  and 
in  debugging  the  program. 

GASP  employs  a  next  et'ent  time  flow  mechanism. 

An  event  selection  function  is  used  to  sequence  the  exe¬ 
cution  oi  event  routines,  The  simulation  switches  from 
event  to  event  according  to  the  times  events  are  scheduled 
to  occur.  The  analyst  must  write  event  routines  to 
specify  system  state  changes  that  occur  at  event  times 
and  the  future  events  that  are  generated  by  event  occur¬ 
rences.  An  event  file  is  used  to  store  attributes,  such 
as  occurrence  times,  and  events  are  inserted  in  the  file 
and  removed  from  it  in  chronological  order. 

Basically  GASP  views  the  world  as  consisting  of 
elements,  attributes,  events,  decision  rules,  processes, 
states,  and  values.  Elements  may  be  permanent  or  tempor¬ 
ary  and  they  interact  through  events.  Events  are  FORTRAN 
subroutines  written  by  the  user  which  change  one  or  more 
elements , 

The  primary  reference  for  GASP  is  37  . 
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GPSS 


The  General  Purpose  Simulation  System  (GPSS)  develops 
the  simulation  model  in  terms  of  a  block  diagram.  The 
structure  and  action  of  a  system  is  described  using 
blocks  which  represent  a  step  in  the  movement  of  entities 
through  the  system. 

One  advantage  of  GPSS  is  that  no  knowledge  of  com¬ 
puter  operation  is  assumed.  The  analyst  need®  only  a 
block  diagram  using  the  standard  GPSS  blocks  to  directly 
transcribe  these  into  the  proper  punch  card  formats.  The 
simulation  is  then  ready  to  be  run.  A  second  advantage 
is  that  GPSS  is  easy  to  learn.  Thirdly,  it  is  a  powerful 
tool  for  problems  that  fit  a  block  structure. 

The  disadvantages  include  limited  arithmetic  state¬ 
ments,  requirement  for  a  separate  compiler,  and  slow  exe¬ 
cution  time. 

The  GPSS  time  flow  mechanism  has  a  next  event  orien¬ 
tation.  Temporary  entities  are  called  transactions  and 
may  encounter  time  delays  in  any  block.  When  that  hap¬ 
pens,  the  transaction  is  placed  in  a  Future  Events  Chain 
in  ascending  order  of  departure  from  the  block.  Other 
transactions  are  in  a  Current  Events  Chain.  The  Current 
Events  Chain  is  scanned  at  each  clock  time  to  attempt  to 
move  transactions  into  the  next  block.  If  movement  is 
possible,  the  transaction  moves  through  as  many  blocks 
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as  possible  until  Cl)  it  encounters  a. time  delay  and 
is  placed  into  the  Future  Events  Chain,  (2)  is  blocked 
from  entering  a  next  block  and  remains  in  the  Current 
Events  Chain,  or  (3)  is  destroyed.  After  the  current 
events  chain  is  completely  scanned  without  being  able 
to  move  airy  transaction,  the  clock  is  updated  to  the 
block  departure  of  the  first  transaction  in  the  Future 
Events  Chain. 

GPSS  views  the  world  as  consisting  of  blocks,  trans¬ 
actions,  and  equipment.  The  simulation  model  is  developed 
in  terms  of  a  block  diagram  of  the  flow  of  transactions 
through  the  system.  A  specific  set  of  block  types  exists 
with  which  the  model  can  be  built.  Program  flow  is  achieved 
by  linking  elemental  blocks  together. 

The  primary  references  for  GPSS  are  39  and  40. 

SIMSCRIPT  II 

SIMSCRIPT  II  is  a  powerful  simulation  language  which 
can  be  used  as  a  GPL  in  its  own  rigV:.  SIMSCRIPT  is  de¬ 
fined  on  several  levels;  this  allows  the  analyst  to  adopt 
as  many  features  as  required  for  the  particular  problem 
at  hand. 

The  advantages  of  SIMSCRIPT  include  more  flexibility 
and  power  for  modifying  system  states  than  other  languages, 
emphasis  on  manipulation  of  entity  attributes,  explicit 


163 


use  of  temporary  entities ,  and  use  of  set  operations. 

In  addition,  the  language  is  English* like.  For  example, 
some  statements  are: 

ADD  E  TO  V 

SUBTRACT  E  FROM  V 

DEFINE  A  AND  B  AS  REAL  VARIABLES 

LET  S  -  X  *  Y 

SCHEDULE  A  NEW.  ARRIVAL  IN  6  MINUTES 
The  diagnostic  capabilities  of  SIMSCRIPT  are  quite  good 
and  some  syntax  errors  are  corrected  by  the  compiler. 

The  chief  disadvantage  is  the  complexity  of  the 
language  and  the  necessity  for  another  compiler. 

The  time  flow  mechanism  is 'hext  event."  The  system 
mechanism  observes  data  cards  and  previously  scheduled 
events ,  orders  them  by  event  ti/ne,  and  causes  them  to 
occur  when  the  time  arrives.  The  events  are  called  as 
subprograms,  as  in  GASP.  The  basic  unit  of  action  is 
an  activity  and  an  event  is  an  instant  in  time  which 
starts  or  stops  an  activity. 

SIMSCRIPT's  view  of  the  world  is  based  on  the  notion 
that  the  state  of  a  s?  stem  can  be  described  in  terms  of 
entities,  attributes,  and  sets  (groups  of  entities). 

The  analyst  must  explicitly  specify  all  permanent  and 
temporary  entities  with  a  complete  list  of  attributes 
and  set  memberships. 
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SIMULA 


SIMULA  is  not  as  widely  used  as  the  other  languages. 
It  has  a  block  structure  based  on  that  of  the  pareftt 
language  ALGOL,  where  both  static  and  dynamic  aspects 
are  found  in  the  block.  Each  block  describes  a  component 
of  the  system.  SIMULA  is  mentioned  here  only  to  acquaint 
the  reader  with  its  existence.  The  primary  reference  is 
38. 

Considerations  in  Selecting  a  Language 

Some  considerations  in  selecting  a  particular  pro¬ 
gramming  language  are: 

•  Availability 

•  Programming  costs 

•  Execution  costs 

•  Validity  of  results 
Programming  costs  are  affected  by 

•  Experience  of  the  programmer 

•  Ease  of  programming  in  the  language 

•  Ease  of  debugging 

•  Transferability  between  machines 

•  Availability  of  standard  routines  and  functions 

•  Flexibility  in  using  time  flow  mechanisms 
Execution  costs  are  determined  by 

•  Efficiency  of  the  compiler 

a  Time  flow  mechanism  suitability  to  the  problem 


•  Storage  requirements 

Validity  of  the  results  are  influenced  by 

•  Modularity  of  the  resultant  prograa 

•  "Naturalness"  of  the  language  with  respect  tc 
the  problem  simulated 

•  Tolerance  by  the  compiler  of  non-standard 
constructs 

•  Diagnostic  features  and  error  checking  routines 

•  Quality  of  the  random  number  generators 

Final  Comments  on  Languages 

The  general  advantage  of  an  SPL  over  a  GPL  is  that 
it  helps  the  analyst  crystallize  his  thinking  and  reduce 
model  formulation  time.  The  SFL  can  provide  a  way  of 
thinking  about  the  problem  which  can  be  a  great  asset 
if  the  world  view  is  appropriate  to  the  problem.  Second 
arily,  an  SPL  eases  the  burden  of  obtaining  data  and 
statistics,  and  provides  extensive  error  checking  to 
prevent  misuse  of  the  language  and  t  xocate  logic 
errors . 
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SECTION  VIII 


MILITARY  APPLICATIONS  OF  COMPUTER  SIMULATION 

The  early  history  of  operations  research  is  deeply 
rooted  in  military  applications.  Likewise,  computer 
simulation  has  been  adopted  as  a  major  evaluation  tool 
by  military  decision  makers.  Although  many  problems  such 
as  logistics,  inventory,  and  maintenance  are  common  to  • 
both  the  military  and  industrial  environment,  several 
uses  are  uniquely  military.  In  this  section  we  will 
examine  several  areas  of  military  application  of  simula¬ 
tion. 

An  Overview 

Let's  briefly  examine  three  areas  where  simulation 
is  used  as  a  significant  analysis  tool.  One  important 
area  is  aircraft  design.  Although  computer  models  are 
used  throughout  all  design  phases,  they  are  gaining  in¬ 
creasing  importance  in  the  preliminary  or  conceptual 
phase.  The  use  of  models  allows  engineers  to  establish 
performance  specifications,  evaluate  design  proposals, 
and  reconfigure  subsystems  to  lessor,  vulnerability  of 
components . 

Another  area  of  significant  application  is  nuclear 
survivability/vulnerability  analysis.  Survivability 
in  a  nuclear  environment  is  extremely  important  for 
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military  aircraft  and  weapon  systems.  Digital  models  are 
used  to  calculate  the  effects  due  to  a  nuclear  burst  and 
evaluate  the  survivability/vulnerability  of  systems  lo¬ 
cated  in  the  vicinity.  Effects  usually  include  those  due 
to  neutrons,  gammas,  induced  gammas,  blast,  thermal,  and 
X-rays . 

A  third  case  is  simulation  of  the  Air  Force  logistics 
system.  In  an  era  of  complex  and  sophisticated  weapons 
systems,  extensive  maintenance  and  inventory  facilities 
are  required.  Digital  models  are  used  to  determine  the 
resource  requirements  for  current  and  projected  Air  Force 
equipment  and  weapon  systems. 

There  are  many  other  fields  which  use  simulation  ex¬ 
tensively  bit  the  intention  here  was  to  merely  highlight 
several  widespread  areas .  References  31  and  35  are 
excellent  bibliographies  of  current  simulation  literature 
and  indicate  the  wideranging  application  of  modeling. 

Computer  Simulation  in  Electronic  Wo* fare 

Another  area  of  frequent  application  of  simulation  to 
military  opeiations  research  is  in  electronic  warfare 
analysis.  Electronic  warfare  (EW)  is  defined  as  military 
action  involving  the  use  of  electromagnetic  energy  to 
determine,  exploit,  reduce,  or  prevent  hostile  use  of 
the  electromagnetic  spectrum  and  action  which  retains 


168 


friendly  use  of  thAt  spectrum,  We  will  examine  simulation 
in  the  division  of  EW  known  as  electronic  countermeasures 
(EcM) ,  i,  e,,  actions  taken  to  prevent  or  reduce  an 
enemy's  effective  use  of  the  electromagnetic  spectrum. 

The  decision  maker  is  faced  with  the  problem  of  deter¬ 
mining  which  basic  concept  the  attacking  force  should 
employ:  ignore  the  defense  and  accept  any  consequent 
losses;  avoid  the  defense;  deny  information  tc  the  defend¬ 
er's  sensors,  such  as  through  the  use  of  radar  absorbing 
materials,  reactive  or  evasive  tactics,  or  exploitation 
of  terrain  masking  and  clutter;  degrade  the  operation  of 
the  defense  system  with  noise  jamming,  deception  techniques, 
or  confusion;  or  destroy  the  defense  sensors  and  weapons. 
Simulation  is  used  to  aid  the  decision  maker  in  selecting 
a  concept  to  employ. 

Three  types  of  simulation  are  used  in  ECM  analysis  — 
digital  models,  simulators,  and  flight  tests.  A  simula¬ 
tor  refers  to  a  man-in- the- loop  simulation  and,  within 
an  EW  context,  is  defined  as  a  man-machine  tool  which 
combines  operators,  analog/digital  computers,  end  actual 
hardware,  such  as  real  radar  displays  or  operational  ECM 
components,  to  simulate  the  real  world  system.  Flight 
tests  are  tests  involving  real  aircraft  and  simulated 
threat  radars  with  associated  operators  which  are  per¬ 
formed  in  an  instrumented  or  controlled  environment. 
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Figure  31.  The  Structure  of  Simulation  in 


Therefore,  in  considering  the  structure  of  simulation 
for  ECM  analysis  we  have  Figure  31,  a  special  case  of 
Figure  2  in  Section  I, 

The  Advantages  and  Limitations  of  Modeling 

In. choosing  among  digital  models,  simulators,  and 
flight  tests,  the  decision  maker  is  faced  with  distinctive 
advantages  and  limitations.  Let's  look  at  some  specific 
uses  of  digital  models  and  simulators.  Both  can  be  used 
to  verify  analytical  results,  obtain  results  not  amenable 
to  purely  analytical  solution,  and  conduct  parametric 
and  sensitivity  studies.  They  also  can  provide  a  base 
for  flight  test  design,  verify  instrumentation  require¬ 
ments  for  flight  tests,  evaluate  equipment  specifications, 
optimize  equipment  usage,  and  evaluate  tactics.  The  ad¬ 
vantages  of  digit  il  models  and  simulators  include  rigid 
test  control,  ease  in  making  configuration  changes,  lower 
cost  than  flight  tests,  and  the  limited  instrumentation 
required.  Simulators  enhance  realism  over  pure  digital 
models  by  including  human  interactions  and  actual  ECM 
hardware  or  radar  equipment.  The  limitations  on  both 
models  include  restrictions  on  numbers  of  aircraft  and 
ECM,  less  realism  than  flight  tests,  and  lack  of  many 
real  world  interactions, 

Flight  tests  may  be  used  in  the  same  fashion  as  a- 
bove  for  the  other  models  but  also  can  evaluate  real  e- 
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quipment,  confirm  tactics,  and  evaluate  operational 
readiness.  The  advantages  of  flight  testing  include 
the  capability  to  use  multiple  aircraft  and  multiple 
penetration  aids,  incorporation  of  real  world  phenomena 
and  interactions,  and  high  realism.  The  limitations 
include  extensive  ground  instrumentation,  limited  on¬ 
board  instrumentation,  limited  test  control,  fixed 
geography  of  test  range,  and  restrictions  to  several 
penetrators  against  only  one  radar. 

Table  26  is  a  comparison  of  the  three  simulations 
according  to  several  measures.  It  is  clear  that  tradeoffs 
in  time,  cost,  and  realism  are  involved  in  the  selection 
of  a  particular  simulation  model. 

A  Hypothetical  EW  Study 

To  demonstrate  how  simulation  models  are  used  in  an 
EW  analysis,  consider  the  following  hypothetical  study 
task.  The  overall  objective  is  to  identify  areas  of 
highest  payoff  for  employment  of  ECM  onboard  manned 
aircraft  penetrating  an  air  defense  system  at  high,  medium, 
and  low  altitudes.  In  particular,  the  analyst  would 
examine  the  effectiveness  of  surface-to-air  missiles  and 
antiaircraft  artillery  employed  against  the  penetrators 
and  those  factors  which  are  significant  in  limiting  the 
performance  of  the  defense  system. 
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Several  Important  faatures  should  be  incorporated 
into  the  analysis  study.  The  analysis  should  include 
the  penetration  force's  use  of  tactics  and  ECM,  penetrator 
acquisition  by  enemy  weapons  systems,  enemy  tracking 
system  performance  in  both  clear  and  jamming  environments , 
weapon  trajectory  and  fuzing  errors,  and  vehicle  vulnera¬ 
bility  due  to  warhead  fragmentation. 

A  composite  modeling  approach  could  be  employed  in 
the  study.  Target  acquisition  can  be  modeled  by  digital 
computer  techniques  as  can  simulation  of  missile  guidance, 
missile  and  shell  trajectories,  missile  fuzing,  and 
vehicle  vulnerability.  However,  human  tracking  performance 
should  be  based  on  simulator  and  flight  test  data.  Figure 
32  dep  icts  the  interrelationship  of  various  subtasks  and 
modeling  techniques  utilized  within  the  study. 

In  order  to  make  accurate  judgments  on  the  basis  of 
aircraft  survivability,  the  effects  of  target  tracking 
errors  must  be  transformed  to  weapon  miss  distances. 

Then  by  modeling  the  vulnerable  components  of  the  aircraft 
and  blast  fragmentation  of  the  warhead,  the  probability 
of  killing  the  penetrator  can  be  scribed  as  a  function 
of  weapon  miss  distance.  Figure  33  shows  a  graph  of 
illustrative  output  from  such  digital  model.  By  syste¬ 
matically  evaluating  the  same  defense  configuration  with 
and  without  various  types  of  electronic  countermeasures, 
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Figure  33  .  Typical  Model  Output  Illustrating  Penetrator 
Kill  as  a  Function  of  Missile  Miss  Distance 


a  family  of  probability  curves  could  be  generated  de¬ 
picting  the  effectiveness  of  ECM  against  terminal  weapons 
systems.  The  probability  of  penetrator  kill  can  also 
be  described  as  a  function  of  missile  miss  distance  and 
downrange  distance  from  the  surface-to-air  missile  site. 
Figure  34  illustrates  the  variations  in  vehicle  kill 
probability  as  a  function  of  warhead  detonation  position. 

The  data  resulting  from  an  evaluation  of  tracking 
system  performance,  say  for  antiaircraft  artillery  could 
be  summarized  by  a  cumulative  probability  that  the  shell 
missed  by  less  than  a  given  distance.  Figure  35  shows 
such  a  graph  as  a  function  of  tracking  mode. 

The  three  outputs  above  typify  the  probabilistic 
relationships  which  can  be  derived  from  data  generated 
by  digital  models.  A  systematic  variation  of  various 
parameters,  ECM  techniques,  force  sizes,  or  tactics  would 
result  in  a  thorough  investigation  for  the  proposed  task. 
As  a  result,  an  assessment  can  be  made  of  the  vulnerable 
areas  of  the  defense  system  and  those  functions  against 
which  ECM  can  be  most  effectively  employed. 


b 


177 


TRACKING  MODS 


Figure  3*.  Effectiveness  of  AAA  Tracking  System  Performance 


1.0 


35 .  Penetrator  Kill  as  a  Function  of  Distance 
Pros  Missile  Site 


Computer  Simulation  in  Logistics 

The  effects  on  operations  of  supply  distribution 
policies,  spares  stockage  levels,  and  maintenance  capabil¬ 
ities  is  usually  difficult  or  impossible  to  determine 
using  analytical  techniques.  Demands  for  supplies,  spare 
parts,  and  repair  work  occur  at  random  intervals  and 
frequently  the  demand  rates  cannot  be  correlated  to  other 
activity  levels.  For  example,  the  removal  rates  of  air¬ 
craft  engines  may  decrease  for  a  period  of  time  when 
flying  activity  increases.  In  these  cases,  the  inter¬ 
actions  of  various  logistics  functions  can  best  be  studied 
using  simulation. 

A  goal  of  the  supply  function  is  to  preposition  the 
limited  available  assets  where  they  can  fill  the  greatest 
need  with  the  least  delay.  A  large  portion  of  the  military's 
resources  is  invested  in  a  class  of  items  known  as  recover¬ 
able  spares.  These  items  may  be  removed  from  an  aircraft 
when  they  fail,  replaced  by  a  like  serviceable  unit,  and 
the  failed  unit  is  then  sent  to  a  repair  facility.  Some 
of  the  units  may  cost  up  to  $500,000  each.  For  example, 
if  an  Inertial  Navigation  Unit  (INU)  fails  on  an  aircraft, 
it  must  be  sent  to  a  depot  for  repair.  When  the  depot 
completes  the  repair  of  the  INU,  a  decision  must  be  made 
as  to  where  to  locate  the  unit  until  it  is  needed.  It 
may  require  several  days  to  ship  the  unit  to  a  base. 
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The  depot  may  decide  to: 

a.  Ship  the  unit  to  any  one  of  several  bases  where 
tnc  unit  may  be  needed,  or 

b.  Retain  the  unit  at  the  depot  until  it  is  known 
which  base  needs  it  most. 

At  the  time  the  decision  is  made,  none  of  the  bases  may 
have  a  "desperate"  need  for  one  more  unit.  I£  base  A  has 
the  lowest  stock  on  hand  and  the  unit  is  shipped  to 
base  A,  it  may  develop  that  base  B  uses  up  its  on-hand 
stock  first  and  becomes  "desperate"  for  a  replacement. 

In  this  case,  a  "bad  decision"  has  been  made.  If,  on 
the  other  hand,  the  depot  decides  to  retain  the  unit  and 
base  A  actually  develops  a  need  for  the  unit,  then  the 
delay  of  shipping  time  occurs  before  base  A  receives  the 
unit  and  it  may  appear  that  the  depot  made  a  "bad  decision" 
in  retaining  the  unit.  Clearly  the  only  way  to  make  the 
"right  decision"  every  time  is  to  have  complete  knowledge 
of  the  future.  Lacking  this,  one  must  consider  the  known 
probabilities  of  error  for  various  alternatives  aid  se¬ 
lect  the  cours  of  highest  probable  payoff. 

Several  decision  making  rules  for  this  redistribution 
problem  are  possible.  The  current  policy  in  the  Air  Force 
is  to  establish  an  authorized  stock  level  for  each  type 
of  unit  for  each  base  and  when  a  base  sends  a  unit  to 
the  depot  for  repair,  the  depot  sends  the  base  a  like  unit 
as  replacement. 
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An  alternative  policy  proposed  by  RAND  Corporation  [42] 
is  called  Real  Time  METRIC.  Under  this  policy  a  "need" 
based  on  on-hand  stock  and  expected  demand  rate  is  Com* 
puted  for  each  base.  A  depot  "reluctance"  is  also  computed. 
If  the  "need"  at  one  or  more  bases  exceeds  the  "reluctance" 
then  the  unit  is  shipped  to  the  base  with  the  greatest 
"need".  Otherwise,  the  depot  retains  the  unit  until  the 
next  decision  time. 

Would  this  policy  give  better  results?  Simulation 
techniques  were  used  to  test  ♦he  policy  and  results  indi¬ 
cated  that  a  301  reduction  in  base  backorders  could  occur. 
This  decision  policy  is  being  included  in  the  Advanced 
Logistics  System.  Other  redistribution  policies  are 
also  being  developed  and  tested  through  simulation. 
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TRICI-J  OF  THE  TRADE 

Certain  computer  programming  techniques  are  used 
in  many  different  types  of  simulation  models.  They  pro* 
vide  fundamental  capabilities  and  are,  for  the  most  part, 
simple  but  enlightening  uses  of  elementary  concepts. 

For  the  uninitiated  the  method  of  implementation  may  not 
be  obvious;  therefore,  some  "tricks  of  the  trade"  are 
documented  in  this  section. 

Sorts 

The  programmer  may  be  required  to  sort  an  array  of 
values,  for  example,  into  ascending  order  of  magnitude. 
There  are  at  least  two  techniques—  a  bubble  sort  and  a 
virtual  sort. 

Bubble  Sort 

In  a  bubble  sort  an  array  of  values  is  checked  ele¬ 
ment  by  element  and  compared  to  the  succeeding  value. 

If  the  objective  is  to  order  the  elements  by  increasing 
magnitude,  then  when  a  succeeding  value  is  smaller  than 
its  predecessor  the  elements  are  reversed.  The  FORTRAN 
subroutine  is  shown  in  Figure  36  .  Note  that  if  the 
list  of  values  is  already  completely  sorted,  only  one 
pass  will  be  required. 
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c  *  *  * 


*  *  * 


C  *  *  *  BUBBLE  SORT  *  *  * 

C 

C 

SUBROUTINE  SORT  (A,  ISTART,  ISTOP) 
DIMENSION  A  (100) 

M  -  ISTART  ♦  1 
N  -  ISTOP 

1  ISW  -  1 

DO  2  I  •  M,  N 

•  IF  (A(I-l).LE.A(I))  GO  TO  2 
TEMP  -  A(I-l) 

A(I-l)  -  A(I) 

A(I)  -  TEMP 
ISW  -  0 

2  CONTINUE 
N  -  N-l 

I?  (ISW, EQ. 0)  GO  TO  1 

RETURN 

END 

Figure  36,  Bubble  Sort  FORTRAN  Subroutine 
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Virtual  Sort 


In  a  bubble  sort,  the  original  ordering  of  the  values 
is  lost  since  the  ultimate  effect  is  to  reorder  the  elements 
according  to  increasing  size.  In  a  virtual  sort,  the  end 
result  is  the  same  but  the  original  array  is  not  changed. 

A  new  indexing  array  is  used  as  a  pointer  to  establish  the 
new  ordering.  The  FORTRAN  subroutine  is  shown  in  Figure  37  • 


C  *  *  * 

C  *  * 

c  *  *  * 

c 


10 


2u 


*  *  * 

VIRTUAL  SORT  *  *  * 

*  *  * 


SU BROUT I NE  VIRSORT (A ,M, I START , I STOP ) 
DIMENSION  A(100),  M(100) 

DO  10  I  -  I  -  ISTART,  ISTOP 

M(J)  -  I 

IMAX  ■  ISTOP  -  1 

DO  20  I  -  ISTART,  IMAX 

K  -  I  ♦  1 

DO  10  J  -  X,  ISTOP 
IF  (A(M(I).LE.A(M(J))  GO  TO  20 
ITEMP  -  M(I) 

M(I)  -  M(J) 

M(J)  *  ITEMF 
CONTINUE 
RETURN 
END 


Figure  37.  Virtual  Sort  FORTRAN  Subroutine 
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Search  for  Extrema 


One  of  the  most  elementary  exercises  in  computer 
programming  is  to  determine  the  maximum  or  minimum  from 
among  a  set  of  values.  A  simple  FORTRAN  DO  Loop  and 
a  greater  or  less  than  test  will  suffice.  Figure  38 
is  the  FORTRAN  subroutine  for  finding  the  maximum  value 
in  an  array  and  returning  the  value  of  the  array  index. 

The  only  point  of  significance  is  the  test  IF  (A(I) .LT.XMAX) . 
Note  that  if  two  elements  in  the  array  have  the  same 
value  XMAX,  the  index  of  the  latter  one  will  be  chosen 
as  IMAX.  Hid  the  test  been  IF(A(I) . LE.XMAX) f  then  the 
first  value  encountered  would  be  retained  as  IMAX.  The 
programmer  should  be  aware  of  this,  especially  if  the  in¬ 
dex  has  some  effect  on  results,  such  as  alternative  solu¬ 
tions  . 
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C  *  *  *  SEARCH  FOR  MAXIMUM  * 

C 

SUBROUTINE  FINDMAX  (A,IMAX) 
DIMENSION  A(100) 

XMAX  -  *1.E  10 
DO  10  I  -  1,  100 
IF(A(I) .LT.XMAX)  GO  TO  10 
XMAX  -  A(I) 

IMAX  -  I 
10  CONTINUE 
RETURN 
END 


Figure  38,  Maximua  Search  FORTRAN  Subroutine 


TABLE  LOOK-UP 

If  a  random  variate  is  to  be  selected  from  on$  of 
the  standard  distributions,  it  is  usually  best  to  use  the 
traijformation  methods  given  in  Section  III.  However, 
random  variates  sometimes  are  needed  from  a  distribution 
which  is  difficult  or  impossible  to  describe  using  the 
known  probability  density  functions.  For  example,  sup¬ 
pose  one  desires  to  generate  the  length  (in- hours)  of 
sn  aircraft  sortie  where  5%  of  all  sorties  are  1  hour 
long,  45%  are  7  hours  long,  40%  are  9  hours  long,  and 
10%  are  24  hours  long. 

See  Figure  39. 

.40 
.30 
prob . . 20 
.10 
0 

0123456789  10  -  --  24 

Sortie  Length  (hrs) 

Figure  39.  Distribution  of  sortie  lengths 

In  this  case  a  table  look-up  procedure  should  be  used. 
First  a  table  must  be  established  thus: 

F  (X)  X 

t5£  t 

.50  7 

.90  9 

1.00  24 
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To  use  this  table  to  generate  .mdom  sortie  lengths,  the 

FORTRAN  function  in  Figure  40  is  called. 

FUNCTION  SORLEN  (DUMMY) 

DIMENSION  F  (4),  HOURS  (4) 

DATA  F/. 05,. 50,. 90, 1.00/, HOURS/1. 0,7. 0,9. 0,24.0/ 

R  -  KANG(RNUM) 

DO  .  1-1,4 

IF  (F(I).GE.R)  GO  TO  2 

1  CONTINUE 

2  SORLEN  -  HOURS  (I) 

RETURN 

END 

Figure  40.  A  FORTRAN  function  for  table  look-up. 

Thus,  if  the  random  number  generated  is  0.62,  the  random 
variate  9.'  will  be  returned. 

A  variation  of  this  technique  involves  interpolation. 
If  the  length  of  sorties  is  described  by  the  graph  in 


Figure  41,  the  function  in  Figure  42  would  be  used. 


Figure  41.  Cumulative  distribution  of  sortie  lengths 
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FUNCTION  SORLET  (DUMMY) 

DIMENSION  F(9),H(P) 

DATA  F/.0,.l,.l,.3,.6,,75,,75,.95,1.0/ 

DATA  H/0. ,2. ,6. ,7. ,8. ,9. ,24. ,25. ,26./ 

R  »  RANG  (RANUM) 

DO  1  I  -  1,9 

IF (F (I) .GE.R)  GO  TO  2 

1  .  CONTINUE 

2  SORLET  -  H(I-l)  +  (H(I)  -  H(I-1))*(R-F(I-1)) 

1  /  (F(I)-F(I-l)) 

RETURN 

END 

Figure  42.  Table  look-up  with  interpolation. 


Thus,  if  the  generated  random  number  is  0.20,  the  returned 
sortie  length  will  he  6.5  hours. 
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Tablo  27.  Continuous  Probability  Distributions 
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Table  28.  Discrete  Probability  Distributions 
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Table  29.  ASCII  Character  Set 
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GLOSSARY 


The  fol lowing  glossary  Is  a  selection  of  frequently  occurring 
simulation  terms  whose  definitions  are  essential  to  understand  or 
even  converse  In  the  military  operations  research  world. 

The  computer  terminology  Is  based  on  the  A)P  Glossary. 

NAVSO  P-3097,  Department  of  the  Navy,  December  1970. 
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absolute  address 

(1)  An  address  that  Is  permanently  assigned  by  the  machine 
deslqner  to  a  storage  location. 

(2)  A  pattern  of  characters  that  Identifies  a  unique  storage 
location  without  further  modification. 

access  time 

(1)  The  time  Interval  between  the  Instant  at  which  data  are 
called  fc»*  from  a  storage  device  and  the  Instant  delivery 
begins. 

(2)  The  time  Interval  between  the  instant  at  which  data  are 
requested  to  be  stereo  and  the  instant  at  which  storage 
Is  started. 

accumulator 

A  register  In  which  the  result  of  an  arithmetic  or  loqle 
operation  Is  formed. 

accuracy 

The  degree  of  freedom  from  error,  l.e.,  the  degree  of  conformity 
to  truth  or  to  a  rule.  Accuracy  Is  contrasted  with  precision, 
for  example,  four  place  numerals  are  less  precise  than  six 
place  numerals,  nevertheless  a  properly  computed  four  place 
numeral  might  be  more  accurate  than  an  Improperly  computed 
six  place  numeral. 

address 

An  Identification  as  represented  by  a  name,  label,  or  number 
For  a  register,  location  In  storage,  or  any  other  dat.i  source. 


ADP 

Automatic  Data  Processing 


air  defense 

All  defensive  measures  designed  to  destroy  attacking  enemy 
aircraft  or  missiles  In  the  earth's  atmosphere. 


ALGOL 

ALGOrlthmlc  Languaqe.  A  language  primarily  used  to  express 
computer  programs. 

alqorltlim 

A  prescribed  set  of  all  the  letters  In  a  language.  Including 
letters  with  diacritical  signs  where  appropriate;  punctuation 
marks  are  not  part  of  an  alphabet. 


alphameric 

same  as  alphanumeric 
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alphanumeric 

Pertalnlni  to  a  character  set  that  contains  letters ,  digits, 
and  usually  other  characters  such  as  punctuation  marks. 

Synonymous  with  alphameric. 

American  Standard  Code  for  Information  Interchange 

A  standard  seven-olt  coded  character  set  developed  by  the 
American  National  Standards  Institute  (ANSI)  to  be  used  for 
Information  Interchange  among  Information  processing  systems, 
communications  systems,  and  associated  equipment.  Abbreviated 
ASCII.  Same  as  OSASCII. 

analog 

Pertaining  to  representation  by  means  of  continuously  variable 
physical  quantities.  Contrast  with  digital. 

analoq  comouter 

( i }  A  computer  In  which  analog  representation  of  data  Is  mainly 
used. 

(2)  A  computer  that  operates  on  analog  data  by  performing  physical 
processes  on  these  data.  Contrast  with  digital  computer. 


analysis 

The  methodical  Investigation  of  a  problem,  and  the  seporttlcn 
of  the  problem  Into  smaller  related  units  for  further  detailed 
study. 

analyst 

A  person  who  defines  problems  and  develops  ilgorlthms  and  procedures 
for  their  solution. 

arithmetic  unit 

The  unit  of  a  computing  system  that  contains  the  circuits  that 
perform  arithmetic  operations. 


array 

An  arrangement  of  elements  In  one  or  mcr«  dimensions, 
artificial  intelligence 

The  capability  of  a  device  to  perform  functions  that  arc*  normally 
associated  with  human  Intelligence,  such  as  reasoning,  learning, 
and  self-improvement.  Related  to  machine  learning. 

artlflclil  languane 

A  language  based  on  a  set  of  prescribed  rules  that  are 
established  prior  to  Its  usage.  Contrast  with  natural  language. 

artillery 

Complete  projectile-firing  weapons  consisting  of  cannon  or 
missile  launchers  on  suitable  carriages  or  mounts.  Field 
artillery  cannons  are  classified  according  to  caliber  as: 
11ght-120mm  and  less  med1um-121-160mm,  heavy-1 61 -21 Omm,  very 
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heavy-greater  than  210mm. 
assemble  ' 

To  prepare  a  machine  language  propram  from  a  symbolic  lanquane 
proqram  by  substituting  absolute  operation  codes  for  symbolic 
operation  codes  and  absolute  or  relocatable  addresses  for 
symbolic  addresses. 

assembler 

A  computer  program  that  assembles, 
assembly  language. 

(1)  A  machine-oriented  programming  language  which  belongs  to  an 
assembly  program  or  system. 

(2)  In  writing  Instructions  using  an  assembly  language,  the  pro¬ 
grammer  Is  primarily  concerned  with  a  label  flild,  *n 

oper:  Jon  field,  and  an  operand  field.  It  Is  possible 
to  relate  the  symbolic  coding  to  Its  associated  flowchart 
If  desired,  by  appending  comments  to  each  Instruction  line 
or  program  segment. 

attenuatl on 

Decrease  In  Intensity  of  a  signal,  beam,  or  wave  as  a  result 
of  absorption  of  enerqy  and  o'  scattering  out  of  the  pa*h  o* 
a  detector,  but  not  Including  the  reduction  due  to  geometric 
spreading,  l.e.,  the  inverse  square  of  distance  effect. 

attrlbuter 

Descriptive  parameters  associated  with  entitles, 
attrition 

The  reduction  of  the  effectiveness  of  a  force  caused  by  loss  of 
personnel  and  material. 

batch  processing 

Is  processing,  without  unscheduled  Interruption,  of  a  group  of 
Items  prepared  or  required  for  one  or  more  related  operations. 

BCD 

Binary-coded  decimal  notation. 

bias 

The  amount  by  which  the  average  of  a  sot  of  values  departs  from 
a  reference  value. 

binary 

(1)  Pertaining  to  a  characteristic  or  property  Involving  a 
selection,  choice,  or  condition  In  whtbh  there  are  two 
possibilities. 
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(2)  Pertaining  to  the  number  representation  system  with  a 
radix  of  two. 

binary  code 

A  code  that  makes  use  of  exactly  two  distinct  characters,  usually 
0  and  1. 

binary  coded  decimal  notation 

Positional  notation  In  which  the  Individual  decimal  digits 
expressing  a  number  In  decimal  notation  are  each  represented  by 
a  binary  numeral*  e.g.,  the  number  twenty»three  Is  represented 
0010  0011  In  the  8-4-2-1  type  of  binary-coded  decimal  notation 
and  by  1011  In  binary  notation.  Abbreviated  BCD. 

binary  search 

A  dichotomizing  search  In  which  the  number  of  Items  of  the  set 
Is  divided  Into  two  equal  parts  at  each  step  of  the  process. 
Appropriate  adjustments  ai'e  usually  made  for  dividing  an  odd 
number  of  Items. 


A  binary  digit.  Same  as  shannon.  See  check  bit.  Information 
bits,  parity  bit,  sign  bit. 


block 

(1)  A  set  of  things,  such  as  words,  characters,  or  digits 
handled  as  a  unit. 

(2)  A  collection  of  contiguous  records  recorded  as  a  unit. 

Blocks  are  separated  by  block  gaps  and  each  block  may 
contain  one  or  more  records. 

block  diagram 

A  diagram  of  a  system,  Instrument,  or  computer  In  which  the 
principal  parts  are  represented  by  suitably  associated 
geometrical  flqures  to  show  both  the  basic  functions  and  the 
functional  relatloshlos  among  the  parts.  Centrist  with 
flowcharts. 

buffer 

A  routine  or  storage  used  to  compensate  for  a  difference  In  rate 
of  flow  of  data,  or  time  of  occurrence  o *  events  when 
transmitting  data,  from  oee  device  to  another. 


bug 

byte 


A  mistake  or  malfunction. 

A  sequence  of  adjacent  binary  digits  used  to  represent  a 
character  (usually  6  or  8  bits  long). 
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card  column 

(1)  A  single  line  of  punch  positions  parallel  to  the  short 
edge  of  a  3  1/4  by  7  3/8  Inch  punched  card. 

(2)  Positions  In  a  line  parallel  to  the  vertical  edqe  of  the 
card. 

card  punch 

A  device  to  recprd  Information  In  cards  by  punching  holes  In 
the  cards  to  represen  letters,  digits,  and  special  charatters. 

card  reader 

A  device  which  senses  and  translates  Into  Internal  form  the  holes 
In  punched  cards. 

cathode  ray  tube  display 

A  device  that  presents  data  In  visual  form  by  means  of  controlled 
electron  beams.  Abbreviated  CRT  display. 


chad 

The  piece  of  material  removed  #hen  formlnq  a  hole  or  notch  1r«  a 
storage  medium  such  as  punched  tape  or  punched  cards.  Synonymous 
with  chip. 

character 

A  letter,  dlqlt,  or  other  symbol  that  Is  used  as  part  of  the 
organisation,  control,  or  representation  of  data.  A  character 
Is  often  In  the  form  of  a  spatial  arrangement  of  adjacent  or 
connected  strokes, 

characteristic 

The  Integral  part  of  a  logarithm.  For  example.  In  the  expression, 
log  643-2.808,  the  .808  Is  the  mantissa  and  the  2  Is  the 
characteristic. 

character  recognition 

The  Identification  of  graphic,  phonic,  or  ether  characters  by 
automatic  means. 

closed  loop 

A  loop  from  which  there  Is  no  exit  other  than  by  Intervention 
from  outside  the  program. 

closed  ship 

Pertaining  to  the  operation  of  a  computer  facility  In  which  most 
productive  problem  programming  Is  performed  by  a  group  of  programming 
specialists  rather  than  the  problem  orlqlnators.  The  use  of 
the  computer  Itself  may  also  be  described  as  cleded=shoo  If  full 
time  trained  operators,  rather  than  user/programmers  serve  as  the 
operators.  Contrast  with  open  shop. 


code 


A  set  of  unambiguous  rules  specifying  the  way  In  which  data 
may  be  represented,  e.g.,  the  set  of  correspondences  In  the 
standard  Synonymous  with  coding  scheme. 

Command  and  control  system 

The  facilities,  equipment,  commumt cations,  procedures,  and 
personnel  essential  to  a  comnander  for  planning,  directing, 
and  controlling  operations  of  assigned  forces  pursuant  to  the 
missions  assigned. 

command  language 

A  source  language  consisting  primarily  of  procedural  operators, 
each  capable  of  Invoking  a  function  to  be  executed. 

common  business  oriented  langua*v 

A  specific  language  by  which  business  data  processing  procedures 
may  be  precisely  described  In  a  standard  form.  The  language  Is 
Intended  not  only  as  a  means  for  directly  presenting  any 
business  program  to  any  suitable  computer,  for  which  a  compiler 
exists,  but  also  as  a  means  of  communicating  such  procedures 
among  Individuals.  Commonly  referred  to  as  COBOL. 


corocn  fluid 

A  field  that  can  be  accessed  by  two  or  more  Independent  routines, 
compile 

To  prepare  a  machine  language  program  from  a  computer  program 
written  In  another  programming  language  by  making  use  of  the 
overall  logic  structure  of  the  program,  or  generating  more  than 
one  machine  Inttructlon  for  each  symbolic  statement,  or  both,  as 
well  as"pe#€erm1ng  the  function  of  an  assembler. 


compiler 

A  program  that  compiles, 
complement 

A  number  that  can  be  derived  from  a  specified  number  by  stbtractlng 
It  from  a  second  specified  number.  For  example.  In  radix  notation, 
the  second  specified  number  may  be  a  given  power  of  the  radix  or 
one  less  than  a  given  power  of  the  radix.  The  negative  of  a  number 
Is  often  represented  by  Its  complement. 
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computer 

A  data  processor  that  can  perform  substantial  computation. 

Including  numerous  arithmetic  or  Ionic  operations,  without 
Intervention  by  a  human  operator  durlno  the  run. 

computer  program 

A  series  of  Instructions  or  statements,  In  a  form  acceptable 
to  a  computer,  prepared  In  order  to  echleve  a  certain  result. 

computer  word 

A  sequence  of  bits  or  characters  treated  as  a  unit  and  capable 
of  being  stored  In  one  computer  location.  Synonymous  with 
machine  word. 

conditioned  »vent 

Events  ui-nd  to  trigger  events  caused  by  a  new  system  state, 
congruence 

The  congruence  x«y  (mod  m)  Is  read  X  Is  oonqruent  to  y  modulo 
ni .and  means  that  (x-y)  is  divisible  by  m. 

console 

That  part  of  a  computer  used  for  communication  between  the 
operator  or  maintenance  enqlneer  and  the  compu 

control  card 

A  punched  card  containing  Input  data  or  parameters  for  Initializing 
or  modifying  a  program, 

control  proqram 

A  collective  or  nenerai  term  for  all  routines  In  the  ooer»t\ng 
system  that  contribute  to  the  management  of  resources.  Implement 
the  data  organization  or  communications  conventions  of  thfc 
operating  system,  or  contain  privileged  operations. 


copy 

To  reKoeduce  data  In  a  new  location  or  other  destination,  leaving 
the  source  data  unchanged,  although  the  physical  form  of  tne 
result  may  differ  from  that  of  the  source.  For  example,  to  copy 
a  deck  of  cards  onto  a  magnetic  tape.  Contrast  with  duplicate. 

counter 

A  device  such  as  a  register  or  storage  location  used  to  represent 
the  number  of  occurrences  of  an  e*ent. 


CPU 

Central  Processing  Unit. 
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data 

A  representation  of  facts,  concepts,  or  Instructions  in  a 
formalized  manner  suitable  for  communication.  Interpretation, 
or  processing  by  humans  or  automatic  Hans. 


debug 

To  detect,  locate,  and  remove  mistakes  from  a  routine  or 
malfunctions  from  a  computer. 


deck 

A  collection  of  punched  cards.  Synonymous  with  card  deck, 
diagnostic 

Pertaining  to  the  detection  and  Isolation  of  a  malfunction  of 
mistake. 


digit 

(1)  A  symbtl  that  represents  one  or  the  non-negative  Integers 
smaller  than  the  radix.  For  example.  In  decimal  notation, 
a  digit  Is  one  of  the  characters  from  0  ttt  9. 

(2)  Synonymous  with  numeric  character. 

digital  computer 

(1)  A  computer  In  which  discrete  representation  of  data  Is  mainly 
used. 

(2)  A  computer  that  operates  on  discrete  data  by  performing 
arithmetic  and  logic  processes  on  these  data.  Contrast 
with  analog  computer. 

documentation 

(1)  The  creating,  collecting,  organizing,  storing,  citing,  and 
disseminating  of  documents  or  the  Information  recorded  In 
documents. 

(2)  A  collection  of  documents  or  Information  on  a  given  subject, 
double  precision 

Pertaining  to  t.nv.  use  of  two  computer  wordssto  represent  a  number. 


dummy 

Pertaining  to  the  characteristic  of  having  the  appearance  of  a 
specified  thing  but  not  having  the  capacity  to  function  as  such. 
For  example,  a  dummy  character,  dummy  plug,  or  a  dummy  statement. 


dump 

To  copy  the  contents  of  all  or  part  of  a  storage,  usually  from 
an  Internal  storage  Into  an  external  storage. 

dynamic  programming 

In  operations  research,  a  procedure  for  optimization  of  a 
multistage  problem  solution  wherein  a  number  of  decisions  are 
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available  at  each  stage  of  the  process, 
electromagnetic  radiation 

Radiation  made  up  o<'  oscillating  electric  and  (Magnetic  fields  and 
propagated  with  the  speed  of  light.  Includes  gamma  radiation. 
X-rays,  ultraviolet, .visible  and  Infrared  radiation,  and  radar 
and  radio  waves. 

electronic  counter-countermeasures 

That  major  subdivision  of  electronic  warfare  Involving  actions 
taken  to  insure  our  own  effective  use  of  electromagnetic 
radiations  despite  the  enemy's  use  of  countermeasures . 

electronic  countermeasures 

That  major  subdivision  of  electronic  warfare  Involving  actions 
taken  to  prevent  or  reduce  the  effectiveness  of  enemy  equipment 
and  tactics  employing  or  affected  by  electromagnetic  radiations 
and  to  exploit  the  energy's  use  of  such  radiations. 

electronic  deception 

The  deliberate  radiation,  ('eradiation,  alteration,  absorption, 
or  reflection  of  electromagnetic  radiations  In  a  manner  Intended 
to  mislead  an  enemy  In  the  Interpreatlon  of  data  received  by  his 
electronic  equipment  or  to  present  falser  Indications  to  electronic 
systems. 

electronic  jamming 

The  deliberate  radiation,  reradiation,  or  reflection  of 
electromagnetic  signals  with  the  object  of  Impairin'*  the  use 
of  electronic  devices  by  the  enemy. 

electronic  warfare 

Military  action  Involving  the  use  of  electromagnetic  energy  to 
determine,  exploit,  reduce  or  prevent  hostile  use  of  the 
electromagnetic  spectrum  and  action  which  retains  friendly 
use  of  the  electromagnetic  spectrum.  Alsu  called  EW.  There 
are  three  dlvlslors  within  electronic  warfare: 

1.  electronic  warfare  support  measures— That  division 
of  electronic  warfare  Involving  actions  taken  to 
search  for.  Intercept,  locate,  and  Identify 
Immediately  radiated  electromagnetic  enerqy  for  the 
purpose  of  Inrvnedl ate  threat  recognition.  Thus, 
electronic  warfare  support  measures  provide  a  source 
of  Information  required  for  Immediate  action  Involving 
electronic  countermeasures,  electronic  counter*counter- 
measures,  avoidance,  targeting,  and  other  tactical 
employment  of  forces.  Also  called  ESM. 
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2.  electronic  countermeasures— That  division  of  electronic 
warfare  Involving  actions  taken  to  prevent  or  reduce  an 
enemy's  effective  use  of  the  electromagnetic  spectrum. 

Also  called  ECM.  Electronic  countermeasures  include: 

a.  electronic  Jamming— The  deliberate  radiation,  re- 
radiation,  or  reflection  of  ilectromagnetic  energy 
devices,  equipment,  or  systems  being  used  by  an 
enemy. 

b.  electronic  deception— The  deliberate  radiation, re- 
radiation,  alteration,  absorption,  or  reflection  of 
electromagnetic  energy  it'  a  manner  Intended  to  mislead 
an  enemy  -n  the  interpretation  or  use  of  Information 
received  by  Ms  electronic  systems.  There  are  two 
categories  of  electronic  deception:  (1)  manipulative 
deception— The  alteration  or  simulation  of  friend’ y 
electromagnetic  radiations  to  accomplish  deception. 

12)  Imitative  deception- -The  Introduction  of  radiations 
Into  enemy  chanrels  which  imftite  his  own  emissions. 

3.  electronic  counter-countermeasures— That  division  of  electronic 
warfare  Involving  actions  taken  to  Insure  friendly  effective 
use  of  the  electromagnetic  spectrum  despite  the  enemy' f  use 

of  electronic  warfare.  Also  called  ECCM. 

emulate 

To  Imitate  one  system  with  another  such  that  the  Imitating 
system  accepts  the  same  data,  executes  the  same  programs,  and 
achieves  the  same  results  as  the  Imitated  system.  Contrast 
with  simulate. 

entity 

An  object  by  which  the  system  can  be  defined,  l.e.,  a  component 
of  the  system  represented  In  the  model. 

error  correcting  code 

A  code  In  which  each  acceptable  expression  conforms  to  specific 
rules  of  construction  that  also  define  one  or  more  equivalent 
nOnacceptable  expressions,  so  that  ff  certain  errors  occur  In 
an  acceptable  expression  the  result  will  be  one  of  Its 
equivalents  and  thus  the  error  can  be  corrected. 

error  detecting  code 

A  code  In  which  each  expression  conforms  to  specific  rules  of 
construction,  so  that  If  certain  errors  occur  In  an  expression, 
the  resulting  expression  will  not  conform  to  the  rules  of 
construction  and,  thus,  the  presence  of  the  errors  Is  detected. 
Synonymous  with  self-checking  codet 

Euler's  function 

In  modular  arithmetic  tha  number  of  positive  Integers  less  than 
m  and  relatively  prime  to  m.  Denoted  ${m)  .  For  m  prime,  d(m)» 
m-1 ,  r 
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event 

Points  of  initiation,  alteration,  or  conclusion  of  activities. 


Fibonacci  series 

A  series  of  Integers  In  which  each  inteqer  is  equal  to  the  sum 
of  the  two  preceding  Integers  in  the  series.  The  series  Is 
formulated  mathematically  by  xf-x^i+x*..?.  where  Xn»0,  xi»l, 
i.e.,  0,1,1,2,3,5,8,13,21...  -  1 


field 

In  a  record,  a  specified  area  used  for  a  particular  category 
of  data,  e.g.,  a  group  of  card  columns  used  to  represent  a  wage 
rate  set  of  bit  locations  In  a  computer  word  used  to  express  the 
address  of  the  opened. 


RSFD 

First-in-First-out  priority  test. 

first  geneeation  computer 

A  computer  utilizing  vacuum  tebe  components. 

fixed  point  arithmetic 

A  method  of  calculation  In  which  operations  take  place  In  an 
Invariant  manner,  and  in  which  the  computer  does  not  consider 
the  location  of  the  radix  point.  This  Is  illustrated  by  desk 
calculators  or  slide  rules,  with  which  the  operator  must  keep 
track  of  the  decimal  point.  Similarly  with  many  automatic 
computers,  in  which  the  location  of  the  radix.point  is  the 
programmers  responsibility. 


flgg 

(1)  Any  of  various  types  of  indicators  used  for  Identification, 
e.g.,  a  wordmark. 

(2)  A  character  that  signals  the  occurrence  o^  some  oondltion, 
such  as  the  end  of  a  word. 

(3)  Synonymous  with  mark,  sentinel,  tag. 

Flight  Similator 

An  electronic  device  to  simulate  the  entire  characteristics  of  the 
flight  and  operation  of  military  aircraft.  It  is  used  to  test 
and  check  out  pilots  in  the  operation  of  aircraft  and  the  use 
of  electronic  equipment  employed  in  actual  flight  tactics  and 
combat  operation. 

floating  point  arithmetic 

A  method  of  calculation  which  automatically  accounts  f e  •  the 
location  of  the  radix  point.  This  usually  is  accomplished  by 
handling  the  number  as  a  signed  mantissa  tfcpes  the  radix  raised 
to  a  Integral  exponent,  I.e.,  the  decimal  number  to  an  Integral 
•xnonent,  i.e.f  the  decimal  number  tc  an  Integral  exponent,  i.e., 
;he  decimal  number  *88.3  might  be  written  .083x102;  the  binary 
lumber  +.0011  as  +.11x22. 
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flowchart 

A  graphical  representation  for  th?  definition,  analysis,  or 
solution  of  a  problem.  In  which  symbols  are  used  to  represent 
operations,  data,  flow,  equipment,  etc.  Contrast  with 
block  diagram. 

FORTRAN 

FjRihila  TRANs latlng  system.  A  language  primarily  used  to  express 
computer  programs  by  arithmetic  formulas. 

general  purpose  computer 

A  computer  that  Is  designed  to  handle  a  wide  variety  of  problems, 
ground  controlled  Interception 

A  radar  technique  which  permits  control  of  friendly  aircraft  or 
guided  missiles  for  the  purpose  of  effecting  Interception. 

guided  missile 

An  unmanned  vehicle  moving  above  the  surface  if  the  earth,  whose 
trajectory  of  flight  path  Is  eapablemof  being  altered  by  an  external 
or  Internal  mechanism. 

guided  missile  (air-to-air) 

An  air-launched  gelded  missile  for  use  against  air  targets, 
gv.ided  missile  (alr-to-surface) 

An  air-launched  guided  missile  for  use  against  surface  targets. 

guided  missile  (surface-to-surface) 

A  surface-launched  guided  missile  for  use  against  surface  targets. 

hal f-word 

A  contiguous  sequence  of  bits  or  characters  which  comprises 
half  a  computer  word  and  Is  capable  of  being  addressed  as  a  unit. 

hard  copy 

A  prttfced  copy  of  machine  output  In  a  visually  readable  form,  for 
example,  printed  reports,  listings,  documents,  etc. 

hardware 

Physical  equipment,  as  opposed  to  the  computer  program  or  method 
or  use,  e.g., mechanical,  magnetic,  electrical,  or  electronic 
devices.  Contrast  with  software. 

header  card 

(1)  A  card  that  contllns  Information  related  to  the  data  In 
cards  that  follow. 

(2)  A  card  that  contains  supplemental  Information  related  to  the 
data  on  the  succeeding  card(s).  Contrast  with  trailer  card. 
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heuristic 

Pertaining  to  exploratory  methods  of  problem  solving  1h  which 
solutions  are  discovered  by  ^valuation  of  the  progress  made 
toward  the  final  result.  Contrast  with  algorithm. 


hexadecimal 

Same  as  sexadecimal. 


high  order  position 

Pertaining  to  the  weight  or  significance  assigned  to  the  digits 
of  a  number ,  l.e.,  In  the  number  267349  the  highest  order  digit 
Is  2,  the  lowest  order  digit  it  9.* 

Hollerith 

Pertalnlrg  to  a  particular  type  of  code  or  punched  card 
utilizing  12  rows  per  column  and  usually  80  columns  per  card. 

Initialize 

To  set  counters,  switches,  and  addresses  to  zero  or  other  starting 
valees  at  the  beginning  of,  or  at  prescribed  points  In,  a 
con.outer  routine. 


Input 

Pertaining  to  a  device,  process,  or  channel  Involved  In  the 
Insertion  of  data  or  states,  or  tp  the  data  or  states  Involved. 

Instruction 

A  statement  that  specifies  an  operation  and  the  values  or 
locations  of  Its  operands.  In  this  context,  the  Perm  Instruction 
Is  preferable  to  the  terms  command  or  order  which  are  sometimes 
used  synonymously. 

Integer  programing 

In  operations  research,  a  class  of  procedures  for  locating  the 
maximum  o**  minimum  of  a  function  subject  to  constraints,  where 
some  or  all  variables  must  have  Integer  values. 

Interceptor 

A  manned  aircraft  utilized  for  Identification  and/or  engagement 
of  airborne  objects. 

Interface 

A  shared  boundary.  An  Interface  might  be  a  hardware  component 
to  link  two  devices  or  It  might  be  a  portion  of  storage  or 
registers  accessed  by  two  or  aoee  computer  programs. 

Interpreter 

(1)  A  computer  program  that  translates  and  executes  each  source 
language  statement  oefore  translating  and  executing  the 
next  one. 
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(2)  A  device  that  prints  on  a  punched  card  the  data  already 
punched  In  khe  card. 


An  abbreviation  for  input/output, 
job  control  statement 

A  statement  in  a  job  that  Is  used  In  identifying  the  job  or 
describing  Its  requirements  to  the  operating  system. 

keypunch 

A  keyboard  actuated  device  that  punches  holes  in  a  card  to 
represent  data. 

language 

A  set  of  representations,  conventions,  arid  rules  used  to  convey 
Information. 

left- justify 

(1)  To  adjust  the  printing  positions  of  characters  on  a  page 
so  that  the  left  margin  of  the  page  Is  regular. 

(2)  By  extension,  to  shift  the  contents  of  a  register  so  that 
the  most  significant  digit  is  at  some  specified  position 
of  the  register.  Contrast  with  normal  size. 

(3)  To  align  characters  horizontally  so  thattthe  left-most 
character  of  a  string  Is  In  a  specified  position. 

Life  cycle  cost 

The  total  cost  of  an  Itam  or  system  over  its  full  life.  It 
Includes  the  cost  of  development,  acquisition,  ownership 
(operation,  maintenance,  support,  etc.)  and,  where  applicable, 
disposal. 

LIFO 

Last- In-FI rst-out  priority  basis, 
line  printer 

A  device  that  prints  all  characters  of  a  line  as  a  unit.  Contrast 
with  character  printer. 

linear  programming 

In  operations  research,  a  procedure  for  locating  tne  maximum  or 
minimum  of  a  linear  function  of  variables  that  are  subject 
to  linear  constraints.  Synonymous  with  linear  optimization. 
Abbreviated  LP. 


load 

In  programming,  to  enter  data  Into  storage  or  working  registers. 
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load- and -go 

An  operatinq  technique  In  which  there  are  no  stops  between 
the  loddlng  arid  execution  phases  of  a  program,  and  which  may 
include  assembling  or  compiling, 

logical  record 

A  collection  of  Items  Independent  of  their  physical  environment. 
Portions  of  the  same  logical  record  may  be  located  In  different 
physical  records. 


loop 

A  sequence  of  Instructions  that  Is  executed  repeatedly  until  a 
terminal  condition  prevails. 

low  order' position 

Pertaining  to  the  weight  of  significance  assigned  to  the  digits 
of  a  number,  l.e.,  In  the  number  396148,  the  low  order  digit 
Is  8. 

machine  Instruction 

An  Instruction  that  a  machine  can  recognize  and  execute, 
machine  language 

(1)  A  language  that  Is  used  directly  by  a  machine. 

(2)  The  set  o*  Instructions  expressed  In  the  number  system 
basic  to  a  computer,  toqether  with  symbolic  operation 
codes  with  absolute  addresses,  relative  addresses,  or 
symbolic  addresses. 

macro  Instruction 

An  Instruction  In  a  source  language  that  Is  equivalent  to  a 
specified  sequence  of  machine  Instructions. 

magnetic  core 

A  configuration  of  magnetic  material  that  1r,  or  Is  Intended 
to  be,  placed  In  a  spatial  relationship  to  current-carrying 
conductors  and  whose  magnetic  properties  are  essential  to  its 
use.  It  may  be  used  to  concentrate  an  Induced  magnetic 
field  as  in  a  transformer  Induction  coll,  or  armature,  to  . 
retain  a  magnetic  polarization  for  the  purpose  of  storing 
data,  or  for  Its  nonlinear  properties  as  in  a  logic  element. 

It  may  be  made  of  such  ««ter1al  as  Iron,  Iron  oxide,  or  ferrite 
and  In  such  shapes  as  wires,  tapes,  toroids,  rods,  or  thin  film. 

magnetic  disc 

A  flat  circular  plate  with  a  magnetic  surface  on  which  data  can 
be  stored  by  selective  magnetization  of  portions  of  khe  flat 
surface. 
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Magnetic  drum 

A  right  circular  cylinder  with  a  magnetic  surfiee  on  which  data 
can  be  stored  by  selective  magnetization  of  portions  of  the 
curved  surface. 

magnetic  storage 

A  storage  device  that  utilizes  the  magnetic  properties  of 
materials  to  store  data,  e.g.,  magnetic  cores,  tapes,  and  films. 

magnetic  tape 

(1)  A  tape  with  a  magnetic  surface  on  which  data  can  be  stored 
by  selective  polarization  of  portions  of  the  surface, 

(2)  A  tape  of  magnetic  material  used  as  th?  constituent  In  some 
forms  of  magnetic  cores. 

main  frame 

Same  as  central  processing  unit, 
management 

A  process  of  establishing  and  attalnfig  objectives  to  carry  out 
responsibilities.  Management  coaiiits  of  those  continuing  actions 
of  planning,  organizing,  directing,  coordinating,  controlling,  and 
evaluating  the  use  of  men,  money,  materials,  and  facilities  to 
accomplish  missions  and  tasks. 

mantissa 

The  fractional  part  of  a  logarithm.  In  the  expression,  log  643-2.808, 
the  .808  Is  the  mantissa  and  the  2  Is  the  characteristic. 

Markov  chain 

A  probabilistic  model  of  events  In  whlbh  the  probability  of  an 
event  Is  dependent  only  on  the  event  that  precedes  It. 


mask 

(1)  A  pattern  of  characters  that  is  used  to  control  the 
retention  or  elimination  of  portions  of  another  pattern 
of  characters. 

(2)  A  filter. 

mass  storage  device 

A  device  havlnq  a  large  storage  capacity,  e.g.,  magnetic  disc, 

magnetic  drum. 

mathematical  model 

A  mathematical  representation  of  a  process,  device,  or  concept, 
matrix 

(1)  In  mathematics,  a  two-dimensional  rectangular  array  of 
quantities.  Matrices  are  manipulated  In  accordance  with 
the  rules  of  matrix  algebra. 
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(2)  In  computers,  an  array  of  any  number  of  dimensions. 


median 

An  average  of  a  series  of  quantities  or  values;  specifically, 
the  quantity  or  value  of  that  Item  which  Is  so  positioned  In 
the  series,  when  arranged  In  order  of  numerical  quantity  or  value, 
thfct  there  are  an  equal  number  of  items  of  greater  magnitude  and 
lesser  maqnltude. 


memory 

Same  as  storage. 


mini  computer 

Generally  It  Is  a  small,  general-purpose  digital  computer  with 
a  central  processor  and  core  memory  (approximately  4096  vsrds), 
and  weighs  about  85  pounds.  The  cabinet  holding  these  two 
components  would  be  about  the  size  of  an  electric  typewriter. 

It  has  a  small  word-size,  12  bit  or  less,  and  Is  economically 
priced,  usually  from  $4,000  to  $13,000.  Mostly  they  are  used 
In  the  on-line,  real-time  environment  and  are  bull*  Into  larger 
systems  as  special-purpose  data  reducers  and  controllers.  Current 
machines  vary  In  word  length  ,  Input/output  facilities,  Instruction 
sets,  software,  and  performance. 


mnemonic  symbol 

A  symbol  chosen  to  assist  the  human  memory,  e.q.,  an  abbreviation 
such  as  "mpy:  for  multiply. 


model 

Any  representation  of  a  real  world  system. 


Monte  Carlo  method 

A  method  of  obtaining  an  approximate  sol-.  ...on  to  a  numerical 
problem  by  the  use  of  random  . c ,  „  the  random  walk 
method. 


nanosecond 

A  billionth  of  a  second. 


natural  language 

A  lanquaqe  whose  rule?,  refs  vet  and  describe  current  usage  rather 
than  prescribe  asaga.  Contrast  with  artificial  language. 


noise 

(1)  Random  variations  of  one  or  more  characteristics  of  any  entity 
such  as  voltage,  currant,  or  data. 
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(2$  A  random  signal  of  known  statistical  properties  of  amplitude, 
distribution,  and  spectral  density* 

(3)  Loosely,  any  disturbance  tending  to  interfere  with  the  normal 
operation  of  a  device  or  system. 

object  code 

Output  from  a  compiler  dr  assembt>w»wh1ch  Is  Itself  executable 
machine  code  or  Is  suitable  for  processing  to  produce  executable 
machine  code. 

object  program 

A  fully  compiled  or  assembled  program  that  Is  ready  to  be  loaded 
Into  the  computer.  Synonymous  with  target  program.  Contrast 
with  source  program. 

OCR 

Optical  character  recognition. 

octal 

(1)  Pettalnlng  to  a  characteristic  or  property  Involving  a 
selection,  choice  or  condition  in  which  there  are  eight 
possibilities. 

(2)  Pertaining  to  the  number  representation  system  with  a  radix 
of  eight, 

offline 

(1)  Pertaining  to  equipment  or  devices  not  under  control  of  the 
centtal  processing  unit. 

(2)  Descriptive  of  a  system  and  of  the  peripheral  equipment 

of  devices  In  a  system  In  which  the  operation  of  peripheral 
equipment  Is  not  under  the  control  of  the  central  processing 
unit. 

online 

(1)  Pertaining  tt  equipment  or  devices  under  control  of  the 
central  processing  unit. 

(2)  Pertaining  to  a  user's  ability  tu  Interest  with  a  computer . 
open  shop 

Pertaining  to  the  operation  of  a  computer  facility  In  which  most 
productive  problem  programming  Is  performed  by  the  problem 
originator  rather  than  by  a  group  of  programmlgg  specialists. 

The  use  of  the  computer  Itself  may  also  be  described  as  open 
shop  If  the  program  user/programmer  also  serves  as  the  operator, 
rather  than  a  full  time  trained  operator.  Contrast  wtth  closed 
shop. 

opttatlng  system 

Software  which  controls  the  execution  of  computer  programs  end 
which  may  provide  scheduling,  debugging,  Input/output  control, 
accounting,  compilation,  storage  assignment,  data  management,  and 
related  services. 
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operations  research 

The  use  of  the  scientific  method  to  provide  criteria  for 
decisions  concerning  the  actions  of  peopM,  machines,  and 
other  resources  In  a  system  Involving  repeatable  operations. 
Synonymous  with  operations  analysis. 

optical  character  recognition 

The  machine  Identification  of  printed  characters  through  use  of 
light-sensitive  devices.  Contrast  with  magnetic  ink  character 
recognition.  Abbreviated  OCR. 


order 

In  modular  arithmetic,  the  order  of  X  modulo  m  Is  the  least 
positive  exponent  n  with  xn«*1  (mod  M)  where  x  and  m  are 
relatively  prime. 

output 

Pertaining  to  a  device,  process,  or  channel  Involved  In  an 
output  process,  or  to  the  data  or  states  Involved. 

overflow 

(1)  That  portion  of  the  resdlt  of  an  operation  that  exceeds  the 
capacity  of  the  Intended  unit  of  storage. 

(2)  Pertaining  to  the  generation  of  overflow  *;  In  (1). 

(3)  Contrast  with  underflow. 

overlay 

The  technique  of  repeatedly  using  the  same  blocks  of  Internal 
storage  diring  different  stages  6f  a  program.  When  one  routine 
Is  no  longer  needed  In  storage,  another  routine  can  replace  all 
or  part  of  It. 


pack 

To  compress  data  In  a  storage  medium  by  taklnq  advantage  of  known 
characteristics  of  the  data,  in  sucu  a  way  that  the  original 
data  can  bto  recovered,  e.g.,,to  compress  eata  In  a  Storage 
medium  by  making  use  of  bit  or  byte  locations  that  would  otherwise 
go  unused. 

page 

A  segment  of  a  computer  program  which  has  a  virtual  address  and 
can  be  located  In  main  storage  or  In  auxiliary  storage. 

paging 

The  scheme  used  to  locate  pages  to  move  them  between  main  storage 
and  auxiliary  storage  or  to  exchange  them  with  pages  of  the  same 
or  other  computer  programs. 
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parameter 

A  variable  that  Is  given  a  constant  value  for  a  specific  purpose 
or  process. 

pebetration 

A  form  of  offensive  maneuver  which  seeks  to  break  through  the 
enemy's  defensive  position,  widen  the  gap  created,  and  destroy 
the  continuity  of  his  postions. 

penetration  aids 

Techniques  and/or  devices  employed  by  aerospace  systems  to  Increase 
the  probability  of  weapon  system  penetration  of  an  enemy  defence. 
Examples  are:  low  altitude  flight  profiles,  trajectory  adjustmeits, 
reduceo  radar  cross-sections  of  attack  vehicles,  Improved  vehicle 
hardness  to  effects  of  defense  engagements,  terrain  avoidance 
radar,  bomber  defense  missiles,  decoys,  chaff,  electronic 
countermeasures ,  etc.  Penetration  aids  are  used  by  an  offensive 
system  to  penetrate  more  effectively  enemy  defenses. 

periphery  and  auxiliary  equipment 

Card  and  tape  Input-output  devices,  printers  for  hard  copy  output 
drums  and  discs  for  auxiliary  memory  storage,  magnetic  tape  to 
microfilm  dvvlces, “optical  character  recognition  equipment 
(which  may  also  be  In  the  computer  category),  cathode  ray  tube 
displays,  transaction  recorders,  magnetic  ink  character 
recognition  equipment,  sophisticated  plotter-display  output 
In  a  photographic  or  pveudophotograph'c  form,  and  photographic 
and  magnetic  sheet,  strip,  or  chip  r,emory  storage  devices. 

Automatic  data  processing  data  transmission  facilities  also 
come  within  this  buoader'deflnltlon  when  connected  to  an  on-base 
computer. 

Polish  notation 

Same  as  preflc  notation. 

possible 

A  term  used  to  qualify  a  statement  made  under  conditions  wherein 
*  some  evidence  eilsts  to  support  the  statement.  This  evident®  Is 
sufficient  to  warrant  mention,  but  insufficient  to  warrant 
assumption  as  true. 

prefix  notation 

A  method  of  forming  mathematical  expressions  In  which  each 
operator  precedes  Its  operands.  For  example,  In  prefii 
notation,  the  expression  "(a  plus  b)  multiplied  by  c"  could 
be  represented  by  +abc.  Synonymous  with  Lukasiewicz  notation, 
parenthesis-free  notation,  Polish  notation. 
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primary  event 

Events  which  may  cause  the  creation  of  a  new  future  event  regardless 
of  the  system  state. 

prlmatlve  root 

In  modular  arithmetic,  the  prlmatlve  root  of  m  Is  a  nuhfcer 
whose  order  Is  equal  to  0  (m),  [Euler's  function] 

probability 

A  mapping  of  the  real  line  Into  the  Interval  [0,1]  on  the  real 
line. 


probable 

A  term  used  to  qualify  a  statement  made  under  conditions'  wherein 
the  available  evidence  Indicates  that  the  statement  Is  factual 
until  there  Is  further  evidence  In  confirmation  or  denial. 

problem  oriented  language 

A  programming  language  designed  for  the  convenient  expression 
of  procedures  used  In  the  solution  6f  a  wide  class  of  problems. 

process 

An  activity  that  requires  time, 
program 

(1)  A  series  of  actions  proposed  In  order  to  achieve  a  certain 
result. 

(2)  Loosely,  a  routine. 

(3)  To  design,  write,  and  test  a  program  as  In  (1). 

(4)  Loosely,  to  write  a  routine. 

programmer 

A  person  mainly  Involved  In  designing,  writing,  arid  testing 
corpfuter  programs. 

programming  flowchart 

A  flowchart  representing  the  sequence  of  operations  1r  a 
program. 

programming  language 

A  language  used  to  prepare  computer  programs, 
pseudovndom  number  sequence 

A  sequence  of  numbers,  determined  by  some  defined  arithmetic 
process,  that  Is  satisfactorily  random  for  a  given  purpose, 
such  as  by  satisfying  one  or  more  of  the  standard  statistical 
tests  for  randomness.  Such  a  sequence  may  approximate  any 
one  of  several  statistical  distributions,  such  as  uniform, 
normal,  or  gaus<1in  distribution. 
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punch  card 

A  card  suitable  for  punching  In  a  pattern  of  holes  to  represent 
data. 

radar 

Radio  detection  and  ranging  equipment  that  determines  the  distance 
and  usually  the  direction  of  objects  by  transmission  and 
return  6f  electromagnetic  energy. 

radar  clutter 

Unwanted  signals,  echoes,  or  Images  on  thd  face  of  the  display 
tube  which  Interfere  with  observation  of  desired  signals. 

radar  coverage 

The  limits  within  which  objects  can  be  detectdd  by  one  or  more 
radar  stations. 

radar  echo 

The  signal  Indication  of  an  object  which  has  reflected  enerfy 
transmitted  by  a  radar. 

radix 

In  positional  representation,  the  Integer,  If  It  exists,  by 
which  the  significance  of  the  digit  place  must  be  multiplied  to  give 
Ihe  significance  of  the  next  higher  digit  place.  For  example, 

In  decimal  notation,  the  radix  of  each  place  1$  ten;  of  the 
fives  place  Is  two*  Synonymous  with  base. 

random  access 

Same  as  direct  access. 

random  variable 

A  function  raapplgg  the  sample  space  Into  the  real  line. 

range 

The  difference  between  the  highest  and  lowest  value  that  a 
quantity  cr  function  may  assume. 

real  time 

Pertaining  to  the  performance  of  a  computation  during  the  actual 
time  that  the  related  physical  process  transpires,  In  order  that 
results  of  the  computation  dan  be  used  In  guiding  the  physical 
process. 

register 

A  device  capable  of  storing  a  specified  amount  of  data  such  as 
one  word. 

replication 

repeating  a  run  of  a  program  for  the  same  combination  of  parameters 
but  with  different  random  variations. 
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roundoff 

To  delete  the  least  s1gnif1caik*41g1t  or  digits  of  a  numeral, 
ana  to  adjust  the  part  retained  In  accordance  with  so m  rule. 

routine 

An  ordered  set  of  Instructions  that  may  hav*  some  general  or  frequent 
use. 


run 

A  single,  continuous  performance  of  a  computer  program  for  a 
specified  combination  of  parameters. 


scan 

To  examine  sequentially,  part  by  part, 
search 

To  examine  a  set  of  Items  for  one  or  more  having  a  desired  property. 

second  generation  computer 

A  computer  us<ng  solid  state  components. 

semantl cs 

The  relationships  between  symbols  and  their  meanings, 
simulate 

To  represent  certain  features  of  the  behavior  of  a  physical 
or  abstract  system  by  the  behavior  of  anothhr  system. 

simulation 

The  representation  of  certain  features  of  the  behavior  of  a 
physical  or  abstract  system  by  the  behavior  6f  another  system, 
e.g.v  the  representation  of  physical  phenomena  by  means  of 
operations  performed  by  a  computer  or  the  represettatlon  pf 
operations  of  a  computer  or  the  representation  of  operations 
of  a  computer  by  those  of  another  computer. 

simulator 

In  electronic  warfare,  a  mar-machine  tool  which  combines  operators, 
analog/digital  computers,  and  actual  hardware  to  simulate  a  real 
world  system. 


skew 

The  angular  displacement  of  a  symbol  or  data  medium  from  the 
Intedded  or  Ideal  placement. 

small  arms 

All  arms,  Including  automatic  weapons,  up  to  and  Including 
.60  caliber  and  shotguns. 
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smooth 

To  apply  procedures  that  decrease  or  eliminate  rapid 
fluctuations  In  data. 

software 

The  program*  and  routines  used  tc  extend  the  capability  of  automatic 
data  processing  equipment.  The  types  of  software  are  as  follows: 

a.  Basic  Software  comprises  those  routines  and  programs 
designed  to  extend  or  facilitate  the  use  of  particular 
automatic  data  processing  equipment,  the  requlrvnent 
for  which  tafcls  into  account  the  design  characteristics 
of  such  equipment.  This  software  Is  usually  provided  by 
the  original  equipment  manufacturer  and  Is  normally 
essential  to  and  a  part  of  the  system  configuration 
furnished  by  h’m.  Examples  of  basic  software  are 
executive  and  operating  programs;  diagnostic  programs; 
compilers;  assemblers,  utility  routines,  such  as  sort- 
merge  and  l.iput/output  oonversion  routines;  file 
management  programs  and  data  management  programs. 

Data  management  programs  are  commonly  linked  to,  and/or 
under  the  control  of,  the  executive  or  operating  programs. 

b.  Application  Software  consists  of  those  routines  and  programs 
designed  by  or  for  automatic  data  processing  equipment 
users  to  accomplish  specific,  mission-oriented  tasks, 

jobs  or  functions  using  the  automatic  data  processing 
equipment  and  basic  software  available.  Applications 
software  may  be  either  general  purpose  packages,  such 
as  demand^deposlt  accounting,  payroll,  machine  tool 
control,  etc.;  or  specific  application  programs  tailored 
to  accomplish  a  single  or  limited  numter  of  users' 
functions,  such  as  base  level  personnel,  depot 
maintenance,  missile  or  satellite  tracking,  etc.  Except 
for  general  purpose  packages  which  are  acquired  directly 
from  software  vendors  or  from  the  orlgWal  eq’Jpment 
manufacturers,  this  type  of  software  Is  normally  developed 
by  the  user,  either  with  In-house  resources  or  through 
contract  servl ces . 


sort 

To  segregate  items  Into  groups  according  to  some  definite  rules. 
Same  as  order. 


source  language 

The  language  from  which  a  statement  is  translated. 


source  progttm 

A  computer  program  written  In  a  soiree  language.  Contrast  with 
objective  program. 
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special  puroose  computer 

A  computer  that  Is  designed  to  handle  a  restricted  class  of  problems, 
statistic 

A  random  variable  whose  values  are  determined  by  sample  data, 
storage 

Pertaining  to  a  device  Into  which  data  can  be  entered,  In  which 
they  can  be  held,  and  from  which  they  can  be  retrieved  at  a 
later  time.  Synonymous  with  memory. 

string 

A  linear  sequence  of  entitles  such  as  characters  of  entitles 
such  as  characters  of  physical  element. 

syntax 

(1)  The  structure  of  expressions  In  a  language. 

(2)  The  rules  governing  the  structure  6f  a  language. 

system 

A  collection  of  Identifiable  parts  capable  of  Interacting  In 
such  a  way  that  the  entire  collection  functions  together. 

table  look-up 

A  procedure  for  obtaining  the  function  value  corresponding  to 
an  argument  from  a  table  of  function  values. 

third  generation  computer 

A  computer  utilizing  solid  loqlc  technology  components, 
time  sbarlnq 

(1)  Pertaining  to  the  Interland  use  of  the  time  of  a  device. 

(2)  A  computer  operation,  under  control  of  an  executive 
routine  Incorporating  a  scheduling  priority  algorithm, 
which  effectively  enables  computer  availability  to  a 
multitude  of  users  virtually  simultaneously. 


translate 

lo  transform  statements  from  one  language  to  another  without 
significantly  changlnq  the  meaning. 

translate* 

A  program  which  translates  from  one  programming  language  Into 
another  programming  language. 

transliterate 

To  convert  the  characters  of  one  alphabet  to  the  corresponding 
characters  of  another  alphabet. 
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truncate 

To  terminate  a  computational  process  in  accordance  w1th*some 
rule,  e.g.,  to  end  the  evaluation  of  a  power  series  at  a 
specified  term. 

underflow 

Pertaining  to  the  condition  that  arises  when  a  machine  computation 
yields  ?  nonzero  result  smaller  than  the  smallest  nonzero 
quantity  that  the  Intended  unit  of  storage  Is  capable  of  storing. 
Contrast  with  overflow. 


user 

Anyone  who  requires  the  use  of  services  of  a  computing  system 

or  Its  products. 

validation 

Testing  the  agreement  between  the  behavior  of  a  simulation  model 

and  the  real  world  system. 

variable 

A  quantity  that  can  assume  any  of  a  given  set  of  values. 

verification 

Ensuring  that  a  simulation  model  behaves  as  the  anilyst  Intends. 

verify 

(1)  To  determine  whether  a  transcription  of  data  or  other 
operation  has  been  accomplished  accurately. 

(2)  To  check  the  results  of  keypunching. 

virtual  address 

*  A  symbol  that  can  be  used  as  a  valid  address  part  but  does  not 

necessarily  designate  an  actual  location. 

vulnerability 

(1)  The  susceptibility  of  a  nation  or  military  force  to  any 
action  by  any  means  through  which  Its  war  potential  or 
combat  effectiveness  may  be  reduced  or  Its  will  to  fight 
diminished. 

(2)  The  characteristics  of  a  system  which  causes  It  to  suffer  a 
definite  degradation  (Incapability  to  perform  the  designated 
mission)  as  a  result  of  having  been  subjected  to  a  certain 
level  of  effects  In  unnatural  (manmade)  hostile  environment. 


war  game 

A  simulation,  by  whatever  means,  of  a  military  operation  Involving 
two  or  more  opposing  forces,  using  rules,  data,  and  procedures 
designated  to  deptit  an  actual  or  assumed  real  life  situation. 
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word 


A  character  string  or  s  bjt  string  considered  as  an  entity: 
normally  the  smallest  entity  which  can  be  fetched  from,  'tored 
In,  or  addressed  In  memory. 


word  length 


A  measure  of  the  size  of  a  word,  usually  specified  In  units 
such  as  characters  or  binary  digits. 


write 

To  record  data  In  a  storage  device  or  a  data  medium.  The 
recording  need  not  be  permanent,  such  as  the  wrltlnq  on  a 
cathode  ray  tube  display  device. 
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Discrete,  196 
Pseudorandom  numbers: 
Definition,  39 
Generation,  39 


,  131 


it*. 


— -  - 


■  i-imwa 


Randomness  tests: 

Chi  square,  50 
Correlation  coefficient,  55 
Coupon  collector,  58 
Gap,  61 

Kolmogorov-Smlrnov,  53 
Permutation,  58 
Poker,  57 

Runs  above  S  below  the  mean,  54 

Runs  up,  57 

Serial ,  60 

Sign,  51,  123 

Spectral,  59 

Random  variates: 

Binomial,  73 
Exponential,  70 
Gamma,  75 
Geometric,  70 
Negative  binomial,  71 
Normal,  76 
Poisson,  74 
Uniform,  72 
Rejection  method,  73 
Replication,  97 

Runs  above  &  below  the  mean,  54 
Runs  up,  57 

Sample  size  determination,  110 
Sampling  procedures: 

Permanent  entitles,  103 
Temporary  entities,  110 
'"carch  for  maximum,  190 
Serial  test,  60 
Sign  test,  51,  123 
Simulation: 

Advantages,  29 
.  Closed  loop,  15 
In  EW,  168 
In  Logistics,  180 
Military  applications,  167 
Model  structure,  11,  169 
Study  guidelines,  30 
Simulation  languages,  157 
GASP,  160 

General  Purpose,  158 
GPSS,  160 
SIMSCRIPT  II,  163 
SIMULA,  165 


Software,  digital  computer,  22 
Sorts: 

Bubble,  187 
Virtual,  189 
Spectral  test,  59 
Steady  state,  99,  102 
Stratified  sampling,  105 
Synthesis,  16 
System  representation,  79 
Systems  analysis,  2 
Systems  engineering,  2 

Table  look-up,  192 
Time: 

Generation  times  for  PRN,  63 
State  of  the  model,  82 
Time  flow  mechanisms,  84 

Fixed  time  increment,  84,  87 
Next  event,  87 
Transformation  effect,  35 
Truncation,  101 

Uniform  variates,  72 

Validation,  127 
Verification,  127 

War  game,  11 
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